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al  solutions,  receiver  and  translator  configurations,  and  effects  of  vehicle 
dynamics  on  receiver  operation  and  results. 
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1.0 


INTRODUCTION 


Use  of  the  Global  Positioning  System  (GPS)  for  Amy,  Navy  and  Air  Force  test 
and  training  ranges'  Time,  Space  and  Position  Information  (TSPI)  requirements 
provides  potential  improvements  in  accuracy  and  mobility  and  less  cost  to  the 
nation  In  comparison  to  current  range  instrumentation  systems.  This  document 
was  prepared  to  provide  members  of  the  Triservices  ranges'  test  and  evaluation 
community  with  an  orientation  to  GPS,  examination  of  the  error  sources  affect¬ 
ing  GPS  performance,  discussion  of  accuracy  measures  developed  to  assess  GPS 
performance,  and  a  variety  of  GPS  applications  and  related  issues  of  interest 
to  the  range  analyst  and  planner.  Means  of  evaluating  accuracy  of  results 
using  the  GPS  capability  and  complementing  it  to  improve  local  accuracy  or  ex¬ 
tend  test  periods  for  particular  ranges'  scenarios  are  discussed. 

GPS  will  be  fully  operational  in  the  late  1980's  with  an  array  of  18  satellites 
that  create,  in  effect,  an  accurate  RF  three  dimensional  grid  all  over  the 
world  and  extending  thousands  of  miles  above  the  surface.  GPS  overcomes  pre¬ 
vious  navigation  systems  limitations  of  insufficient  accuracy,  limited  area, 
and  two  dimensional  results.  The  use  of  GPS  is  not  restricted  to  number  of 
participants,  weather  conditions,  visibility,  or  location  (except  for  range 
requirements  such  as  data  link  line-of-slght) .  Common  use  of  GPS  technology 
will  allow  standardization  of  data  processing  and  interfaces. 


1.1 


GPS  Development  and  Operation 


The  GPS  program  has  been  in  existence  for  approximately  fourteen  years.  It 
was  begun  when  the  United  States  Air  Force,  Army,  Marine  Corps,  Navy  and 
Defense  Mapping  Agency  combined  their  technical  resources  to  develop  a  high¬ 
ly  accurate  space-based  navigation  system.  The  development  of  GPS  is  being 
managed  by  a  Joint  Program  Office  (JPO)  at  the  Air  Force  Systems  Command, 

Space  Division,  in  Los  Angeles.  The  management  team  consists  of  DoD,  DoT 
and  NATO  personnel . 

The  GPS  program  is  being  conducted  through  three  successive  phases.  During 
Phase  I,  GPS  performance  and  feasibility  were  validated.  Currently  in  Phase 
II,  GPS  operational  effectiveness  is  being  validated.  In  Phase  III,  produc¬ 
tion  and  deployment  of  GPS  will  be  accomplished. 

The  GPS  comprises  three  major  segments:  Space,  Control  and  User  as  shown  in 
Figure  1 .1 . 

The  Space  Segment  will  consist  of  a  number  of  NAVSTAR  satellites.  Each  satel¬ 
lite  transmits  an  L1  (1575.42  MHz)  and  an  L 2  (1227.6  MHz)  signal.  L1  is  mod¬ 
ulated  by  a  Precise  (?)  signal  and  a  Clear/Acquisition  (C/A)  signal.  L2  is 
modulated  by  either  a  P  or  C/A  signal .  Superimposed  on  these  spread  spectrum 
signals  is  a  navigation  message  containing  ephemeris,  clock  and  atmospheric 
propagation  correction  data.  Table  1.1  summarizes  satellite  signal  character¬ 
istics.  Each  satellite  is  positioned  In  a  10,900  nautical  mile  circular  orbit 
(with  period  equal  to  12  sidereal  hours)  at  an  inclination  of  55  degrees.  In 
the  fully  operational  configuration,  there  are  18  satellites  in  six  orbital 
planes  with  three  satellites  In  each  plane.  This  configuration  is  designed 
to  provide  a  minimum  of  four  satellites  in  view  to  any  user,  thereby  ensuring 
worldwide  coverage. 

The  Control  Segment  Includes  a  number  of  Monitor  Stations  and  Ground  Antennas 
located  throughout  the  world,  together  with  a  Master  Control  Station.  The 
Monitor  Stations  each  contain  a  GPS  receiver  which  oassively  tracks  all 
NAVSTARs  within  view  and  collects  orbital  data.  Trie  Master  Control  Station 
processes  the  data  collected  by  the  Monitor  stations  and  estimates  orbital 
elements,  clock  parameters  and  atmospheric  propagation  corrections  for  each 
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PSEUDO-RANGE  RATE  DATA. 
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NAVSTAR  satellite.  These  estimates  are  used  to  update  the  navigation  message  - 
of  each  satellite,  and  the  updated  messages  are  transmitted  to  the  satellites 
via  the  Ground  Antennas. 

The  User  Segment  consists  of  Individual  User  Equipment  (UE)  sets  and  support 
equipment.  Each  UE  set  passively  receives  data  transmitted  by  NAVSTAR  satel¬ 
lites  and  provides  navigation  and  time  Information  to  a  host  vehicle.  Receiver 
operation  Is  Illustrated  In  Figure  1.2. 

The  concept  of  GPS  rests  on  the  principle  of  trllateratlon.  That  Is,  when 
the  positions  of  three  points  (suitably  located)  are  known  together  with  the 
distance  (or  range)  from  each  to  a  fourth  point,  then  the  position  of  the 
fourth  point  can  be  determined.  In  GPS,  the  positions  of  the  NAVSTAR  satel¬ 
lites  are  contained  In  their  respective  navigation  messages.  The  User  Infers 
range  to  a  given  satellite  by  measuring  the  propagation  time  of  the  encoded 
radio  signal  transmitted  by  the  satellite.  The  user  Inferred  range  obtained 
by  conversion  of  measured  propagation  time  Is  called  pseudo-range.  Because 
of  the  clock  offset,  the  user  requires  pseudo-ranges  from  at  least  four 
(vice  three)  satellites  to  estimate  his  three  dimensional  position  and  clock 
error  relative  to  GPS  time. 

In  a  completely  analogous  manner,  the  user  obtains  pseudo-range  rates  and 
extracts  his  velocity  and  the  drift  rate  of  his  clock  relative  to  GPS  time  by 
measuring  the  received  (Doppler  shifted)  frequencies  of  the  carrier  signals 
from  four  (or  more)  NAVSTAR  satellites. 


BPS  RECEIVER  SIMPLIFIER  RLOCK  DIAGRAM 


N 

IU 

oc 


SITE  CLOCK 
AND  FREQ 
REFERENCE 


1.2 


Document  Organization  and  Summary 


Sections  2.0  and  3.0  develop  accuracy  evaluation  concepts  and  Section  4.0  con¬ 
tains  several  applications  considerations.  In  Section  2.0,  the  effects  of 
various  GPS  error  sources  are  discussed  with  the  objective  of  establishing  an 
error  budget  for  GPS  satellite  to  receiver  range  measurements.  Two  error 
budgets  (one  for  real-time  and  one  for  post- test  uses)  are  developed  and  the 
GPS  error  specification  is  presented.  The  contributors  to  the  error  budgets 
are  allocated  by  responsibility  to  the  Space  Segment,  the  Control  Segment,  or 
the  User  Segment.  Although  the  error  sources  In  each  segment  are  subject  to 
some  manipulation,  attention  is  focused  primarily  on  error  sources  which  are 
most  directly  affected  by  the  user. 

The  User  Segment  errors  are  strongly  dependent  on  link  margin  (Slgnal-to-Nolse 
Ratio),  code  selection,  method  of  Ionospheric  delay  compensation,  and  whether 
effective  use  Is  made  of  the  available  Doppler  measurement  data.  Because  of 
this  strong  dependence,  the  main  thrust  of  Section  2.0  deals  with  these  Issues. 

The  definition,  relationship  and  Interpretation  of  various  accuracy  measures 
that  can  be  used  to  assess  GPS  navigation  and  time-transfer  performance  are 
presented  In  Section  3.0.  There  exists  a  variety  of  accuracy  measures  cur¬ 
rently  used  to  assess  the  performance  of  navigation  and  weapons  systems. 

These  measures  Include  such  entitles  as  Root-Sum-Square  (RSS)  quantities, 
spherical  and  circular  standard  error,  and  spherical  and  circular  error  prob¬ 
able.  The  proliferation  of  such  measures  raises  questions  regarding  their 
respective  Interpretation  and  mutual  relationships.  In  Section  3.0,  various 
commonly  used  accuracy  measures  are  defined.  The  assumption  of  Gaussian  prob¬ 
ability  distribution  Is  Invoked  to  enable  quantitative  interpretation  of  each 
such  measure  and  Its  relation  to  other  measures.  The  main  thrust  of  Section 
3.0  Is  to  provide  the  analyst  with  some  guidelines  which  will  enable  proper 
selection  and  calculation  of  the  appropriate  accuracy  measure  for  his  particu¬ 
lar  application. 

Section  4.0  contains  eleven  subsections  dealing  with  a  variety  of  consider¬ 
ations  for  GPS  Range  applications.  An  attempt  Is  made  to  relate  the  Inter¬ 
acting  effects  of  the  various  subsections  and  of  the  previous  sections. 
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Subsection  4.1  amplifies  the  limitations  to  real-time  data  processing  and 
possible  refinements  for  post-test  processing.  Subsection  4.2  addresses  re¬ 
lative  and  absolute  accuracy  considerations  and  applications.  "Relative" 
Implies  that  the  position,  velocity  and  time  data  of  test  participant s) 
with  respect  to  each  other  Is  the  significant  Information  as  compared  to 
"absolute"  which  requires  data  In  a  global  coordinate  system.  Relative  Is 
primarily  applicable  to  small  arena  or  local  scenarios.  Global  scenarios 
might  use  relative  techniques  for  some  requirements  such  as  propulsion  and 
guidance  system  performance  evaluation  and  for  scoring. 

Subsection  4.3  provides  Information  about  the  current  and  future  number  of 
satellites  available.  Two  effects  of  major  Importance  for  Range  tests  are 
the  periods  when  a  sufficient  number  of  satellites  are  above  the  horizon 
and  the  accuracy  of  the  results  due  to  their  positions  relative  to  the  test 
participants.  Determination  of  test  participants'  TSPI  using  the  current 
GPS  satellites  Is  restricted  to  a  few  periods  of  a  few  hours  per  day.  Length 
of  the  periods  depends  to  same  extent  on  the  number  of  parameters  to  be  esti¬ 
mated  and  user  latitude  and  longitude.  Subsection  4.4  continues  the  discus¬ 
sion  of  accuracy  considerations  with  the  emphasis  on  two  and  three  dimen¬ 
sional  solutions. 

Complementing  the  current  GPS  satellite  array  with  compatible  ground  stations 
can  extend  test  periods,  make  testing  Independent  of  the  satellites,  Improve 
accuracy  and  llne-of-slght  considerations,  and/or  provide  range  C3  capabili¬ 
ties.  Subsection  4.5  provides  considerations  of  this  concept  referred  to  as 
"pseudo-satellites"  or  "Inverted  GPS". 

Test  vehicle  dynamics  Induce  tracking  errors  and  can  Inhibit  acquisition  of 
data  with  traditional  range  Instrumentation.  GPS  technology  Is  subject  to 
the  same  problems  but  In  different  manners  than  the  range  analyst  Is  accus¬ 
tomed  to  evaluating.  Subsection  4.6  Identifies  these  problems  and  presents 
Information  about  their  effects.  Methods  of  "aiding"  signal  acquisition  and 
track  are  reviewed. 

Subsection  4.7  discusses  uses  of  GPS  receiver  and  translator  configurations. 
WSMC  experience  Includes  use  of  GPS  receivers  during  two  Mlnutaman  IC8M 
launches  and  development  of  a  GPS  translator  equipped  sonobuoy  array  for 


reentry  vehicle  scoring.  A  missile-borne  -translator  Is  also  being  evaluated. 
Section  4.7  discusses  receiver  aiding  and  various  approaches  to  translator 
designs. 

Raw  measurement  data  rather  than  navigation  solution  results  are  preferable 
for  test  and  training  range  data  processing  In  order  to  evaluate  the  output 
and/or  apply  more  sophisticated  processing  techniques.  Section  4.8  discusses 
these  needs  from  the  standpoint  of  making  provision  In  user  equipment  under 
development  for  these  data  requirements.  Another  consideration  the  range 
analyst  must  take  Into  account  is  the  ability  to  record  and/or  transmit  the 
data  since  operational  equipment  modifications  are  frequently  forbidden  or 
restricted. 

Subsection  4.9  continues  discussions  of  GPS  receivers  being  developed  with 
descriptions  of  their  characteristics  and  specifications,  wherever  possible. 

Section  4.10  discusses  applications  of  Kalman  Filter /Smoother  (KFS)  tech¬ 
niques  to  GPS  measurements  and  the  complementary  benefits  of  GPS  data  and 
Inertial  navigation  system  data. 

Subsection  4.11  continues  the  discussion  of  GPS  data  processing  and  develops 
subtle  concepts  In  Doppler  data  processing. 
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ERROR  8U0GETS 


For  each  operational  combination  of  a  code  (C/A  or  P),  a  link  (L-j,  L2,  ...) 
and- a  Space  Vehicle  (SV)  within  view  of  a  properly  equipped  GPS  user,  there 
exist  two  types  of  measurements  available  to  the  user.  These  measurement 
types  are,  respectively,  pseudo-time  delay  and  pseudo-Ooppler  shift  of  a 
signal  propagating  from  the  SV  to  the  user  receiver. 

When  the  pseudo-.time  delay  is  scaled  by  the  speed  of  light,  a  pseudo-range 
measurement  is  obtained.  Denoting  this  measurement  by  R,  it  follows  that 

R  ■  R  +  CTy  +  eR 

where  R  Is  the  range  (from  the  Indicated  location  of  the  SV  to  the  true  loca¬ 
tion  of  the  user  receiver),  C  Is  the  speed  of  light,  Ty  is  the  user  clock  off¬ 
set  relative  to  GPS  time,  and  eR  Is  the  combined  error  In  the  pseudo-range 
measurement  from  all  sources. 

In  a  somewhat  similar  fashion,  the  pseudo-Ooppler  shift  cycles  can  be  counted 
over  an  arbitrary  time  interval  and  converted  to  an  Incremental  pseudo-Ooppler 
range  measurement. 

Introducing  the  Increment  notation 

to 

exit*  *  x  (V  - x  <v 

The  Incremental  pseudo-Ooppler  range  measurement  over  time  Interval  [^ ,  t2] 

Is  expressed  by 

.  t2  *2  t2  *2 

[r]^  ■  DOt<|  +  cCTu^t1  +  ^er^t1  ’ 

where  R,  Ty  and  C  have  the  same  meaning  as  before,  and  ep  is  the  combined  error 
In  the  pseudo-Ooppler  range  measurement  from  all  sources. 

There  exist  numerous  methods  of  processing  pseudo-range  and  pseudo-Ooppler 
range  measurements  to  perform  various  navigation  functions.  Some  of  these 
methods  will  be  explored  In  detail  in  a  later  section.  Here  the  more  Important 


methods  are  examined  in  brief. 


The  most  direct  method  of  GPS  navigation  is  to  utilize  four  or  more  pseudo¬ 
range  measurements  (from  suitable  SVs)  and  solve  for  user  position  and  time 
offset.  And  in  a  similar  fashion,  it  is  possible  to  solve  for  user  velocity 
and  frequency  offset  by  utilizing  four  or  more  incremental  pseudo-Doppler 
range  measurements  which  have  been  approximately  differentiated  to  obtain 
pseudo-range  rates. 

Another  method  (applicable  primarily  to  post-test  processing)  consists  of  two 
steps.  First  the  pseudo-range  and  incremental  pseudo-Doppler  range  measure¬ 
ments  are  differenced  and  averaged  to  obtain  a  more  accurate  initial  value  of 
pseudo-range  to  each  SV.  Then  the  smoothed  initial  value  of  pseudo-range  is 
added  to  the  incremental  pseudo-Doppler  range  measurements  for  each  SV  to  ob¬ 
tain  high  precision  pseudo-range  measurements  at  a  high  data  rate.  This  tech¬ 
nique  (developed  in  Section  4.11)  is  illustrated  in  Figure  2.1. 

Furthermore,  in  addition  to  the  above  navigation  methods  (which  are  absolute 
or  global  in  nature),  there  exist  two  differential  navigation  methods. 

Relative  navigation  is  defined  as  the  determination  of  the  coordinate  dif¬ 
ferences  (In  a  suitable  frame  of  reference)  between  two  different  users 
at  a  comnon  set  of  times.  Delta  navigation  is  defined  as  the  determination 
of  the  coordinate  differences  between  two  or  more  locations  occupied  by  a 
single  user  at  different  times. 

In  any  method  of  navigation,  the  consideration  of  accuracy  performance  is 
generally  a  matter  of  major  concern.  The  assessment  of  navigation  perform¬ 
ance  is  a  two  step  process.  The  first  step,  which  is  the  main  subject  of 
this  section,  consists  of  categorization,  characterization,  and  quantification 
of  the  various  error  contributors  in  the  user  measurement  equations.  The 
second  step,  which  Is  the  subject  of  Section  3.0,  consists  of  quantification 
of  the  process  of  error  propagation  through  the  user  algorithms  together 
with  the  definition  and  calculation  of  appropriate  accuracy  measures  for  the 
purpose  of  performance  analysis. 

The  end  result  of  the  first  of  the  two  steps  in  navigation  performance  analysis 
is  a  budget  for  user  measurement  errors.  This  error  budget  usually  provides 
estimates  of  the  RMS  values  of  the  various  error  contributors ,  as  well  as  the 
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RSS  value  of  all  errors  combined.  In  the  range  domain,  the  actual  value  of 
combined  measurement  error  is  called  the  User  Equivalent  Ranging  Error  (UERE), 
and  the  RSS  value  of  all  errors  combined  is  called  the  one-Sigma  UERE  (1<j  UERE). 

The  GPS  system  specifications  document  (Ref  [4])  contains  an  error  budget 
specification  which  is  presented  in  Table  2.1.  The  errors  have  been  placed 
in  either  Space,  Control,  or  User  Segment  categories.  Each  error  in  the  Space 
and  Control  categories  is  specified  to  have  an  RMS  value  less  than  the  amount 
shown  for  each  such  error  source  for  up  to  24  hours  after  SV  upload.  The 
amounts  shown  for  the  User  category  error  sources  are  based  on  use  of  the 
P-code  for  both  navigation  and  dual  frequency  (L-| ,  L2)  ionospheric  delay  com¬ 
pensation  and  on  use  of  local  measurements  of  temperature,  pressure  and  rela¬ 
tive  humidity  for  tropospheric  delay  compensation.  With  no  compensation, 
ionospheric  and  tropospheric  delay  errors  can  each  vary  approximately  2-20 
meters,  depending  on  elevation  angle  (SV  relative  to  user)  and  atmospheric 
conditions. 

The  error  budget  specification  is  a  useful  guide  for  certain  real-time  applica¬ 
tions  of  GPS.  The  specification  does  not,  however,  indicate  the  potential 
improvement  in  post-test  applications  utilizing  better  ephemeris  and  clock 
data  in  conjunction  with  better  processing  techniques.  Neither  does  the  error 
budget  specification  give  an  indication  of  the  temporal  variabilities  of  the 
constituent  errors  which  significantly  influence  the  performance  of  navigation 
algorithms  which  rely  on  filtering  enhancements.  Since  the  specification 
assumes  use  of  the  P-code,  there  is  no  indication  of  error  magnitudes  for 
C/A-code  users.  Finally,  the  specified  error  budget  does  not  include  noise 
terms  for  Doppler  measurements,  nor  the  variation  of  noise  error  levels  with 
link  margin. 

In  the  remainder  of  this  section,  the  issues  which  have  been  raised  herein, 
and  which  are  not  addressed  in  the  GPS  error  budget  specification,  will  be 
explored  further.  Some  of  the  details  with  respect  to  these  issues,  however, 
will  be  given  in  subsequent  sections  of  the  report.  The  culmination  of  this 
exploration  will  be  the  presentation  of  two  basic  error  budgets.  The  first  of 
these  will  apply  nominally  to  various  real-time  GPS  applications,  while  the 
second  will  be  representative  of  post-test  applications  in  which  ideal  con¬ 
ditions  hold  and  special  signal  processing  techniques  are  employed.  Each  of 
these  budgets  will  be  subject  to  some  manipulation  when  specific  applications 
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SYSTEM  UERE 


are  considered  in  later  sections. 


2.1  Pseudo-Range  Noise  Error 

The  GPS  system  specificationsA(Ref  [4])  states  that  whenever  the  SV  is  above 
a  5°  elevation  angle,  atmospheric  losses  are  less  than  2.0  dB,  and  the  receiv 
ing  antenna  has  gain  equivalent  to  0  DBI  for  a  RH  circular  polarization  (3 
081  for  a  linear  polarization) ,  the  RF  signal  levels  at  the  receiving  antenna 
output  shall  not  be  less  than  the  values  shown  in  Table  2.2. 


TABLE  2.2  RECEIVED  SIGNAL  LEVELS  (SPECIFICATION) 


CHANNEL 

C/A-CODE 

P-CODE 

h 

-160  DBW 

-163  DBW 

L2 

-166  DBW* 

-166  DBW* 

*P  code  only  will  normally  be  transmitted  on  \.^. 


The  available  receiver  noise  power  density,  referred  to  the  output  of  the 
receiver  antenna,  is  expressed  by 


-23 

where  k  *  1.38  X  10  Joule  /  Kelvin  is  the  Boltzmann  constant  and  TR  is 
the  effective  receiver  system  noise  temperature  referred  to  the  antenna  output 
With  an  assumed  noise  temperature  of  58QK,  the  available  noise  power  dens  ity 
is  equal  to  -201  OBW/Hz.  Thus  the  estimated  carrier  power- to-noise  power 
density  ratio  on  each  link  is  given  by  the  following  table. 


TABLE  2.3  CARRIER-TO-NGISE  DENSITY  RATIOS 


The  variance  of  the  pseudo-ranging  error  due  to  receiver  noise  is  expressed  by 


where  W  is  the  code  chip  width,  BR  is  the  bandwidth  of  the  code  receiver/ 
tracker,  and  C/Nn  is  the  carrler-to-noise  density  ratio.  8.  is  assumed  to  be 
1  HZ,  and 


293.3  meters  for  C/A  code, 

29.33  meters  for  P  code. 

Based  on  the  above  stated  link  assumptions,  the  RMS  receiver  noise  errors  for 
the  four  link  and  code  combinations  considered  are  given  in  the  following 
table. 


TABLE  2.4  PSEUDO-RANGING  NOISE  ERRORS  (RMS) 


C/A-CODE 

P-CODE 

m 

m 

aR  (METERS) 

2.614  5.216 

0.369  0.522 

2.2  Pseudo-Range  Ionospheric  Delay  Compensation 


In  order  to  perform  dual  frequency  ionospheric  delay  compensation  (Ref  [18]) 
with  pseudo-range  measurements,  the  measurements  on  the  two  carrier  links  are 
first  differenced.  Thus 


4R  a  R2  *  R^ 


and  if  the  correction  is  to  be  applied  to  link  1,  the  difference  is  multiplied 


where  f.]  and  are  the  respective  link  carrier  frequencies*,  to  obtain  a 
measurement  of  ionospheric  delay  on  link  1. 

Since  the  receiver  noise  terms  in  and  L2  are  uncorrelated,  the  variance 
of  the  noise  in  the  measurement  of  ionospheric  delay  on  is  simply 


V1  ,  2  .  2, 

K  (<tr  +  ) 

R1  R2 


After  compensation  for  ionospheric  delay  on  ,  the  variance  of  the  residual 
measurement  noise  is  given  by 


(1+K)2  c  2  +  K2  a  2 

R,  r2 


A  2  . 
fl  «R,  + 


(f?- 


f|)2 


Note  that 


2  2  2 
0R+O  *  "R,  +  "D, 

because  of  correlation  of  the  noise  terms  in  R-j  and  aR  a  (R2  -  R-|).  Also 
note  that  the  residual  noise  variance  has  the  same  value  whether  L1  or  L2  is 
selected  for  compensation. 

Based  on  the  assumed  link  margins,  previously  stated,  the  RMS  noise  errors 
before  and  after  dual  frequency  ionospheric  delay  compensation  are  summarized 
in  the  following  table. 


For  and  l>2>  f^  =  1575.42  MHz  and  f2  *  1227.6  MHz,  respectively. 
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TABLE  2.5  RMS  NOISE  ERRORS  BEFORE  AND  AFTER  DUAL  FREQUENCY  IONOSPHERIC  DELAY 
COMPENSATION 


C/A-CODE 

h 

<jR  (METERS) 

2.614 

0.369  (BEFORE) 

<jQ  (METERS) 

9.018 

0.988 

aR+0  (METERS) 

10.454 

1.238  (AFTER) 

It  should  be  noted  that  If  variations  In  link  margins  occur  which  preserve 
the  link  ratios  (such  as  variations  due  to  receiver  gain  or  losses  and  those 
due  to  atmospheric  effects),  the  RMS  noise  errors  In  the  above  table  can  be 
compensated  for  these  variations  simply  by  scaling  by  the  Inverse  of  the 
square  root  of  the  carrier-to-noise  density  ratio. 

2.3  Nominal  Real-Time  Error  Budget 

A  nominal  real-time  error  budget  for  GPS  is  now  presented  In  Table  2.6.  In 
this  table,  both  C/A  code  and  P  code  users  have  been  considered,  and  slowly 
and  rapidly  fluctuating  errors  have  been  separately  categorized. 

Multipath  errors  have  been  arbitrarily,  estimated  to  have  an  RMS  value  equal  to 
W/25,  and  have  been  assumed  to  be  correlated  (at  least  insofar  as  group  or 
code  delay)  on  L-|  and  The  estimated  value  corresponds  to  the  system 
specification  (Table  2.1),  but  It  is  recognized  that  multipath  error  Is 
heavily  dependent  on  design  and  application  factors.  The  correlation  assump¬ 
tion  on  the  two  links  explains  why  multipath  error  does  not  apply  to  dual 
frequency  ionospheric  delay  compensation. 

The  nominal  real-time  error  budget  In  Table  2.6  assumes  dual  frequency 
Ionospheric  delay  compensation  on  whichever  code  is  used  for  navigation. 

This  assumption  may  not  hold  In  all  applications,  however.  First,  in  some 
differential  navigation  applications,  it  may  be  possible  to  ignore  Ionos¬ 
pheric  delay  by  simply  treating  It  as  a  common  bias  which  cancels  when  mea¬ 
surement  or  navigation  differences  are  formed.  Second,  and  very  significant 
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la  ERRORS 
(METERS) 


data  rate  and  noise  bandwidth  are  each  assumed  to  be  1  Hz. 


to  C/A  code  users,  the  C/A  code  on  may  be  inhibited  most  or  all  of  the 
time  in  the  operational  GPS.  In  this  event,  C/A  code  users  must  either  fore¬ 
go  Ionospheric  delay  compensation  or  use  alternatives  such  as  modeled  pred¬ 
iction  or  translocation  for  this  purpose. 


Thus  a  modified  real-time  error  budget  in  which  ionospheric  delay  compensa¬ 
tion  Is  considered  Independently  Is  of  Interest  In  certain  applications. 

This  modification  to  the  budget  presented  in  Table  2.6  is  accomplished  by 
replacing  the  RSS  value  In  the  fast  column  with  12.3/1.3  and  by  replacing  the 
lo  system  UERE  with  13.0/4.6.  With  this  modification.  Ionospheric  delay 
compensation  errors.  If  applicable,  must  be  separately  taken  Into  account. 

2.4  Pseudo-Ooppler  Range  Noise  Error  and  Ionospheric  Delay  Compensation 

When  pseudo-range  measurements  are  extracted  from  Doppler  phase,  the  variance 
of  error  due  to  receiver  noise  Is  given  by 


(2Hr  C/Nq 


where  x  Is  the  carrier  wavelength  and  BQ  is  the  bandwidth  of  the  carrier 
phaselock  loop.  Is  assumed  to  be  10  Hz. 

By  analogy  with  code  processing,  the  variance  of  the  noise  In  dual  frequency 
Ionospheric  delay  compensation  of  the  Doppler  phase  measurements  on  is 
given  by 

\  ^  ^l  +  ^  ‘ 

where  K  has  the  same  meaning  as  before.  After  compensation  for  ionospheric 
delay  on  ,  the  variance  of  the  residual  noise  in  the  pseudo-Doppler  range 
measurement  is  given  by 

ar+d  »  ^K)2  'r*  +  ^  % 

4  2.4  2 
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Numerical  values  for  the  above  quantities,  based  on  the  link  margins  in 
Table  2.3,  are  presented  in  the  following  table. 

TABLE  2.7  PSEU00-00PPLER  RANGE  RMS  NOISE  ERRORS  BEFORE  AND  AFTER  IONOSPHERIC 
OELAY  COMPENSATION 


C/A-CODE 

P-C0E 

)E 

h 

L2 

h 

L2 

cr  (METERS) 
r 

0.00085 

0.00219 

0.00121 

0.00219  (BEFORE) 

ad  (METERS) 

0.00363 

0.00387 

(METERS) 

0.00402 

0.00458 

(AFTER) 

From  the  values  in  the  above  table,  estimated  RMS  noise  errors  for  pseudo- 
Delta  range,  ar  ■  [r]£+at  ,  and  pseudo-range  rate,  r  •  ar/at,  can  be  obtained 
by  means  of 


G 


r 


and 


JT"or,  before  Ionospheric  delay  compensation  , 
'/^_<7r+d’  a^ter  ionospheric  delay  compensation, 

JT"cjr  Bg,  before  ionospheric  delay  '■ompensation, 
'/^_ar+d  B0’  after  ionospheric  delay  compensation. 


Numerical  values  for  the  above  quantitites  are  presented  in  the  following 
table  for  L-|  only. 
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TABLE  2.8  RMS  NOISE  ERRORS  FOR  PSEUOO-OELTA  RANGE  ANO  PSEUDO-RANGE  RATE 


C/A  (L1 ) 

P  (L^ 

aAr  (METERS) 

0.00121 

0.00171 

Before  Ionospheric 
Delay  Compensation 

(METERS) 

0.00569 

0.00648 

After  Ionospheric 
Delay  Compensation 

aj.  (METERS/SEC0N0S) 

0.0121 

0.0171 

Before  Ionospheric 
Delay  Compensation 

(METERS/SEC0N0S) 

0.0569 

0.0648 

After  Ionospheric 
Delay  Compensation 

2.5  Ideal  Post-Test  Error  Budget 

In  post-test  applications.  It  Is  natural  to  expect  some  Improvement  In  the 
error  budget  relative  to  that  which  applies  in  real-time.  First,  the  clock  and 
ephemeris  errors  should  be  reduced  by  post-test  filtering  and  smoothing  when 
compared  to  predicted  values  of  these  errors.  Second,  certain  processing  options 
are  available  post-test  which  may  not  be  available  In  real-time.  The  most  im¬ 
portant  of  these  Involves  processing  Doppler  phase  measurements  to  obtain  ultra 
low  noise  pseudo-range  data  at  a  relatively  high  rate  (10  Hz).  To  achieve 
maximum  utility  from  this  Doppler  processing  technique,  continuous  measurements 
(with  no  dropouts  or  Ooppler  cycle  slips)  must  be  simultaneously  available  on 
both  L-|  and  Furthermore,  prior  to  processing  the  pseudo-range  measurements 
for  navigation,  one  pass  must  be  made  through  the  data  to  obtain  the  ionospheric 
delay  correction,  and  a  second  pass  must  be  made  to  smooth  the  Initial  value  of 
pseudo-range. 

An  error  budget  for  post-test  GPS  applications  is  presented  in  Table  2.9. 

This  budget  Is  based  on  the  Ideal  conditions  outlined  above  in  which  the 
ephemeris  and  clock  data  have  undergone  post-test  refinements  by  the  control 
segment;  and  the  user  employs  dual  frequency  ionospheric  delay  compensation 
and  Doppler  smoothing  to  reduce  Initialization  error  in  pseudo-range  measure¬ 
ments  constructed  from  continuous  Doppler  data. 


SEGMENT 


Data  rate  and  noise  bandwidth  are  each  assumed  to  be  10  Hz 


It  will  be  observed  that  Ionospheric  delay  compensation,  receiver  noise,  and 
multipath  errors  are  Indicated  In  both  the  slow  and  fast  columns.  Except 
for  multi  th,  the  numbers  In  the  fast  column  are  from  Table  2.7.  (The  RMS 
multipath  noise  error  (fast  column)  Is  assumed  to  be  x^/25.)  The  values  In 
the  slow  column  are  actually  RMS  bias  errors  induced  by  Doppler  smoothing  of 
the  initial  condition.  These  errors  are  due  almost  totally  to  noise  in  the 
code  measurements  of  pseudo-range,  and  their  magnitude  depends  on  the  length 
of  the  Doppler  smoothing  interval.  Because  the  bandwidth  of  the  code  loop  is 
taken  to  be  1  Hz,  the  RMS  values  of  initital ization  errors  induced  by  Doppler 
smoothing  are  obtained  by  dividing  the  respective  values  for  noise  given  in 
Table  2.6  by  where  T  is  the  duration  (expressed  in  seconds)  of  the  Doppler 
smoothing.  The  values  in  Table  2.9  are  based  on  a  Doppler  smoothing  interval 
of  400  seconds  duration. 

In  those  applications  where  ionospheric  delay  compensation  is  either  not  per¬ 
formed,  or  is  performed  in  some  alternative  fashion  which  is  independent  of 
receiver  noise  error,  the  post- test  budget  in  Table  2.9  must  be  modified  to 
delete  terms  due  to  dual  frequency  ionospheric  delay  compensation.  The 
result  of  this  modification  is  that  the  RSS  values  in  the  slow  column  become 
1.7/1. 6,  the  RSS  values  in  the  fast  column  become  0.0077/0.0077,  and  the  Ict 
system  UEREs  become  1.7/1. 6.  With  this  modification,  ionospheric  delay  com¬ 
pensation  errors,  if  applicable,  must  be  separately  taken  into  account. 

The  values  of  <rAr  and  ct*  computed  in  connection  with  the  real-time  budget 
apply  in  the  post-test  budget  without  change. 

2.6  Summary  of  Error  Budget  Properties 

A  top-level  summary  matrix  is  presented  in  Table  2.10.  This  matrix  shows,  in 
a  qualitative  manner,  the  extent  to  which  various  error  sources  enter  dif¬ 
ferent  navigation  modes  for  both  real-time  and  post-test  applications. 
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3.0 


ACCURACY  MEASURES 


The  theoretical  problem  of  assessing  accuracy  in  any  GPS  application  is  solved 
whenever  the  joint  probability  distribution  is  known  for  the  navigation  errors 
of  interest.  In  the  cases  where  the  UEREs  from  each  SV  are  jointly  Gaussian 
and  the  processing  is  essentially  linear,  the  joint  probability  distribution 
of  the  navigation  errors  is  Gaussian,  and  is  therefore  completely  specified 
by  its  first  and  second  moments  (i.e.,  mean  vector  and  covariance  matrix).  A 
simple  example  here  will  serve  the  purpose  of  illustrating  the  development  of 
first  and  second  moments  of  navigation  errors,  as  well  as  providing  a  convenient 
vehicle  to  introduce  the  terminology  of  accuracy  measures. 

Consider  that  a  user  has  made  pseudo-range  measurements  to  four  suitable  SVs 
and  wishes  to  compute  his  position  and  clock  error.  The  four  pseudo-range 
measurements  can  be  expressed  as  follows 

Rf  -  ^(X-X.)2  +  (Y-Y.)2  +  (Z-Z.)2  +  Ct  +  e.  ,  1-1,2,  3,  4. 

In  this  expression,  i  denotes  the  SV,  R.  is  measured  pseudo-range,  is  UERE, 
(X,  Y,  Z)  are  user  coordinates,  (X^,  Y^,  Z.)  are  indicated  SV  coordinates,  t  is 
user  clock  offset,  and  C  is  the  speed  of  light. 

The  user  solution  for  (X,  Y,  Z,  t),  denoted  by  (X',  Y',  Z',  t')  is  given  by 

R1  -  J  (X'-X.)2  +  (Y'-Y.)2  +  (Z'-Z.)2  +  Ct'  ,  1-1,2,  3,  4. 

Eliminating  R^  from  the  two  sets  of  equations  and  expanding  to  first  order  in 
5X  -  X'-X,  $Y  =  Y'-Y,  <5Z  *  Z'-Z,  and  fit  =  t'-t  leads  to 


(X'-X.)fiX  +  (Y'-Y, )6Y  +  (Z’-Z.)fiZ 

+  Cfit  -  ef  ,  1  ■  1 ,  2,  3,  4. 

The  latter  four  equations  can  be  written  in  matrix  form  as 

H?  *  7  , 


f(X'-Xf)2  +  (Y'-Y.)2  +  (Z'-Z.)2 


where 


is  the  unit  vector  from  the  i-th  SV  to  the  user,  consistent  with  the  solution. 
It  follows  that 
s’  a  , 

and  the  first  and  second  moments  of  ?  are  given  by 
E{?}  *  H“VECF>  *  0  , 

C{8>  a  E{S8T}  »  , 

where  aR  is  the  RSS  UERE  common  to  all  4  measurements,  and  the  4  UEREs  are 
assumed  to  be  statistically  independent. 
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Now 


ox2 

JXY 

®XZ 

Cjxt 

Q 

-< 

X 

<Jy2 

5yZ 

COyt 

ffzx 

aZY 

oz2 

Cazt 

Cfftx 

CctY 

totz 

C2at2 

is  the  symmetric  covariance  matrix  of  the  error  vector  [SX,  $Y,  <SZ,  Cdt]1*. 

In  more  general  applications,  the  processing  may  be  more  involved,  and  the 
error  vector  may  contain  additional  terms  such  as  velocities,  frequency,  etc. 
In  such  cases  the  error  covariance  matrix  is  extended  to  contain  the  variances 
and  cross  covariances  of  each  error  term. 


In  terms  of  the  covariance  matrix  for  the  simple  problem  considered  above, 
several  accuracy  measures  can  be  defined  immediately.  The  trace  of  the 
covariance  matrix  is  given  by 

2 

ag2  *  o^2  +  <jy2  +  OjZ  +  C  at2  , 

and  <7q  is  the  Figure  Of  Merit  (FOM)  for  the  solution.  The  FOM  combines  UERE 
statistics  and  geometry  to  assess  accuracy.  The  RSS  Position  Solution  Error 
is  given  by 


'o^2  +  <jy2  +  a^2 


It  is  important  to  realize  that  by  the  properties  of  orthogonal  matrices, 
both  dg  and  op  are  invariant  with  respect  to  orthogonal  transformations.  Hence 
<jg  and  op  are  independent  of  the  Cartesian  coordinate  system  in  which  the 
solution  is  expressed. 


If  a  particular  coordinate  frame  is  selected  in  which  the  X-Y  plane  is  horizontal 
and  the  Z-axis  is  vertical,  then  the  RSS  Horizontal  Position  Solution  Error  is 
given  by 


2  +  Oy  2 
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and  the  RMS  Vertical  Position  Solution  Error  is  given  by 

av  3  Jz  * 

Each  of  the  quantitites  <jg,  ap,  aH,  <jy,  Cot  is  termed  an  accuracy  measure. 
When  an  accuracy  measure  is  normalized  by  ap,  the  "Dilution  Of  Precision" 

(or  OOP)  resulting  from  geometry  and  processing  is  obtained.  For  the  above 
accuracy  measures,  DOPs  are  defined  as  follows: 

GDOP  3  » 

POOP  3  ap/cR 

HOOP  3  aH/aR  , 

VQOP  3  Oy/aR  ■  and 

TOOP  3  Co^/CTp 

Other  measures  of  accuracy  are  often  defined,  such  as  Spherical  Error  Probable 
(SEP),  Circular  Error  Probable  (CEP),  Spherical  Standard  Error  (c^),  and 
Circular  Standard  Error  (oc).  However  use  of  these  accuracy  measures  is 
generally  restricted  to  applications  In  which  the  error  components  are  jointly 
Gaussian.  These  accuracy  measures  will  be  defined  in  the  sequel,  and  for  the 
Gaussian  case,  relationships  between  the  various  accuracy  measures  will  be 
explored. 

3.1  Three  Dimensional  Error  Probability  Distribution 

Let  X,  Y,  Z  be  zero  mean  jointly  Gaussian  random  variables  with  covariance 
matrix 


which  is  assumed  to  be  nonsingular.  It  can  be  shown  that  there  exists  an 
orthogonal  matrix  l  such  that 

lTCl  *  A  , 

where  a  is  a  diagonal  matrix,  and  moreover,  trC  «  tr.\. 

Defining,  £,  n,  c  by 


Z 


n 


I 


c 


2 


it  follows  that  £,  n,  t  are  zero  mean  independent  Gaussian  random  variables 
with  covariance  matrix 


A 


0 


0  0 

o  2  0 

n 

0  a;2 


Now 


K2 


is  the  equation  of  an  ellipsoid.  With  K  >_  0,  the  probability  that  the  point 
(X,  Y,  Z)  lies  within  the  ellipsoid  is  denoted  by  P^(K;3)  and  is  given  by 
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Pg( K; 3)  Is  called  an  Elliptical  Error  Probability  in  3  dimensions,  and  it  is 
a  simple  matter  to  show  that 


pe(K;3)  »  PpCx3  <  <2]  . 

2 

where  x]  Is  a  “chi -squared"  random  variable  with  3  degrees  of  freedom.  Thus 
elliptical  error  probabilities  can  be  obtained  from  readily  available  tables 
of  x2- probabilities. 

In  contrast,  the  probability  that  the  point  (X,  Y,  Z)  lies  within  a  sphere  of 
radius  R,  denoted  by  P-(R),  is  given  by 


,  .2  2  2 

(22)3/Zo,a  cr  "  d5dn<lc  ‘ 
5  1  C 

^  +  n‘  +  c4*  <  R* 


P-(R)  is  called  a  Spherical  Error  Probability,  and  In  general  bears  no  simple 
^  2 

relationship  to  probabilities  associated  with  x  or  other  comnonly  tabulated 
random  variables.  The  value  of  R  which  results  in  P^  «  0.5  Is  called  the 
Spherical  Error  Probable  and  is  denoted  by  SEP.  Thus 

PS(SEP)  »  0.5  . 

2 

By  reference  to  x  tables  it  Is  seen  that  P^(l;3)  *  0.198748  _  That  is, 

the  probability  that  (X,  Y,  Z)  is  within  the  "lo-ell ipsold"  Is  0.198748  ..., 
and  this  is  the  basis  for  the  definition  of  the  spherical  standard  error. 

The  value  of  R  which  results  In  P^  ■  0.198748  ...  is  called  the  Spherical 
Standard  Error  and  Is  denoted  by  a^.  Thus 

P$(cs)  »  0.198748  ...  *  Pe(1;3)  . 

A  special  case  of  particular  Interest  occurs  when  ae  *  •  i.  and  the  "lc- 

i  * 

ellipsoid"  is  actually  a  sphere.  In  this  case  the  Gaussian  random  variables 
X,  Y,  Z  are  said  to  be  Spherically  Distributed,  and  <?_  »  j.  *  j  «  j_. 


3.2  Two-Oimensional  Error  Probability  Distribution 

By  analogy  with  the  above  development,  if  X  and  Y  are  zero  mean  jointly 
Suassian  random  variables  with  covariance  matrix 


there  exists  an  orthogonal  matrix  L  such  that 
LTCL  -  A  , 

where  A  is  diagonal  and  trC  *  trA  .  With  (s,n)  defined  by 


it  follows  that  £»  n  are  zero  mean  independent  Gaussian  random  variables  with 
covariance  matrix 


a^2  0 

0  a  2 

n 

Assuming  C  is  nonsingular,  it  follows  that 

2  , 

+  »  K2 

a  2 

n 

Is  the  equation  of  an  ellipse.  With  K  >_  0,  the  probability  that  the  point 
(X,  Y)  lies  within  the  ellipse,  denoted  by  P^(K;2),  is  given  by 
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Pe(K;2) 


1 


expT{& +  r?} 


,2  2 

y  +  <  K£ 

V  V  “ 


2!lora 
C  n 


d$dn  . 


P_( K; 2)  is  called  an  Elliptical  Error  Probability  in  2  dimensions,  and 


where  Is  a  "chi-squared"  random  variable  with  2  degrees  of  freedom. 

The  probability  that  the  point  (X,  Y)  lies  within  a  circle  of  radius  R, 
denoted  by  Pq(R)  and  given  by 


PC(R) 


■  // 


lZ  +  n2  <  R2 


2  2 


"P^T?  +  7T} 


d5dn  , 


is  called  a  Circular  Error  Probability.  The  value  of  R  which  results  in 
Pr  *  0.5  is  called  the  Circular  Error  Probable  and  is  denoted  by  CEP.  T 


Now  the  probability  that  (X,  Y)  is  within  the  "la-ellipse"  is  P^( 1 ; 2)  * 
0.393469  ....  and  this  is  the  basis  for  the  Circular  Standard  Error,  denoted 
by  Oq  and  defined  implicitly  by 


Pc(oc)  »  0.393469 


PE(l ;2)  . 


Thus  aj.  is  the  radius  of  the  circle  which  contains  (X,  Y)  with  probability 
0.393469  .... 


In  the  special  case  when  a,  3  a  ,  the  "la-ellipse”  is  a  circle  and  ar  *  a.  * 
In  this  case  the  Gaussian  random  variables  X,  Y  are  said  to  be  Circularly 
Distributed. 


3. 
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One-Dimensional  Error  Probability  Distribution 


I 


In  analogy  with  the  above  developments  for  2  and  3  dimensional  cases,  if  x  is 
zero  mean  Gaussian  random  variable  with  variance  cx2,  define 


Pe(K;1)  * 


/ 


x2  <  k2 

°7- 


JzHc 


J-  '&*2  <x 


and 


Pl(R) 


/ 


x2  «  R2 


v/2nc 


e-x  /2ox2  dx 


Pe(K;1)  is  an  Elliptical  Error  Probability  in  1  dimension  and  PL(R)  is  called  a 
Linear  Error  Probability.  But  it  is  easily  seen  that 

Pe(K;1)  *  PL(Kax)  *  PrCxf  <  k23  , 

where  x-|  is  a  "chi-squared"  random  variable  with  1  degree  of  freedom. 


The  Linear  Error  Probable  (LEP)  is  defined  by 


Pl(LEP)  »  0.5  , 

and  it  follows  that 

LEP  *  (0.675)aL  , 

where  s  ax  is  the  Linear  Standard  Error. 
3.4  Dilution  of  Precision 


The  "Dilution  Of  Precision"  or  OOP  for  each  of  the  spherical,  circular,  and 
linear  accuracy  measures  can  be  defined  in  analogy  with  the  DOPs  which  were 
defined  earlier  for  the  RSS  accuracy  measures  oG,  <jp,  oH,  and  ay.  Thus  for  GPS 
applications,  where  cR  denotes  la  UERE, 
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n 


SEPDOP  *  SEP/ s R  , 

SOOP  • a  as/aR  , 

CEPOOP  *  CEP/cr  , 

COOP  *  aC^aR  * 

LEPOOP  »  LEP/ cR  , 

LOOP  *  a^/aR 

It  is  clear  now  that  a  plethora  of  accuracy  measures  can  be  defined.  Further¬ 
more,  there  are  no  general  algebraic  relationships  between  the  relatively  few 
accuracy  measures  defined  above.  One  thing  is  certain  however.  All  accuracy 
measures  defined  above  are  implicitly  contained  In  the  error  covariance  matrix. 

In  fact,  the  error  covariance  matrix  is  itself  the  ultimate  accuracy  measure  of 
the  types  defined  above.  Moreover,  when  the  errors  are  jointly  Gaussian,  the 
mean  vector  and  covariance  matrix  provide  a  complete  statistical  characterization 
of  the  errors.  Since  the  RSS  quantities  Og,  op,  oH,  oy  can  be  obtained  virtually 
by  inspection  from  the  covariance  matrix,  RSS  errors  are  particularly  convenient 
accuracy  measures. 

3.5  Accuracy  Measure  Relationships 

In  certain  special  cases,  it  is  possible  to  determine  numerical  relationships 
between  variously  defined  accuracy  measures.  These  relationships  will  be 
examined  first  for  the  3  dimensional  case  and  then  for  the  2  dimensional  case. 

For  both  cases  it  is  assumed  that  the  errors  are  jointly  Gaussian. 

Consider  first  the  case  in  which  the  1c  ellipsoid  is  a  sphere,  i.e., 

*  c  »  a  *  j>0  and  the  errors  are  spherically  distributed.  In  this  case, 

5  n  c  , — 


o,  Cp  •  J3  eg,  and  for  K  >_  0, 

Ps(Kcs)  »  PrCx3lK2] 
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For  the  50th  percentile  of  a  X3  random  variable,  it  follows  that  SEP  =  /x|  a 
(1.538  ...)a$. 

When  a  3  a  3  a  >  0  and  a-  3  0,  the  3-dimensional  distribution  is  said  to  be 
K  n  ^  r— 

Circularly  Degenerate.  In  this  case  ap  *y2  a,  and  for  K  <  0, 

PsUa)  3  Pr[X2  i  K2]. 

From  this  relation  and  the  tables,  it  follows  that  3  (0.6657  .  ..)a  and 
SEP  *  (1.1774  ...)a,  respectively. 

Finally  when  a  3  o>  0  and  a  3  a,  3  0,  the  3-dimensional  distribution  is 
s  n  o 

said  to  be  Linearly  Degenerate.  In  this  case  crp  3  a,  and  for  K  >  0, 

P^(Ka)  =  Pp[xf  <  K2]. From  this  relation  it  follows  that  a s  3  (0.2517  .  ,.)a 
and  SEP  3  (0.67448  ...)a. 

For  the  2  dimensional  case,  consider  first  the  case  in  which  the  errors  are 
circularly  distributed,  i.e.,  the  la  ellipse  is  a  circle  with  3  a^  3  a  >  0. 
In  this  case,  ac  3  a,  aH  ac,  and  for  K  >  0  , 

Pc(Kcxc)  3  Pr[X2  <  K2]  . 

2 

From  this  expression  and  the  x£  tables,  it  follows  that  CEP  3  (1.1774  ...)ac- 

In  the  2  dimensional  case  where  a,  3  a  >  0  and  a  3  0,  the  distribution  is 

s  n 

Linearly  Degenerate.  In  this  case  aH  3  a,  and  for  K  >  0, 

Pc(Ko)  3  PrCx2  £  K2]  • 

From  this  relation  it  follows  that  ac  3  (0.515  ...}a  and  CEP  3  (0.67448  ...)a. 

The  accuracy  measure  relationships  for  the  above  special  cases  are  summarized 
in  Tables  3.1  and  3.2. 

As  can  be  seen  immediately  from  these  tables,  the  selected  accuracy  measure 
relationships  vary  significantly  with  the  shape  (i.e.,  eccentricity  and/or 
degeneracy)  of  the  la  ellipsoid  (or  ellipse).  This  is  indicative  that  none 


TABLE  3.1  3  OIMENSIQNAL  ACCURACY  RELATIONS 


CASE 

ffP 

aS 

CP 

SEP 

SEP 

°S 

SPHERICAL  DISTRIBUTION 

1 .732 

1 .126 

1 .538 

CIRCULAR  OEGENERACY 

2.1 24- 

1.201 

1.769 

LINEAR  DEGENERACY 

3.973 

1 .483 

2.680 

TABLE  3.2  2  OIMENSIQNAL  ACCURACY  RELATIONS 


CASE 

ffC 

CEP 

ffC 

CIRCULAR  DISTRIBUTION 

1  .414 

1.201 

1 .177 

LINEAR  OEGENERACY 

1 .942 

1 .483 

1 .310 

of  the  accuracy  measures  discussed  above  (except  the  covariance  matrix  itself) 
can  give  a  complete  assessment  of  all  aspects  of  accuracy.  Consequently,  in 
order  to  gain  clearer  insight  into  the  meaning  and  utility  of  each  scalar 
accuracy  measure,  it  is  necessary  to  relate  each  individually  to  an  appropriate 
probability  distribution.  The  most  useful  distributions  for  this  purpose  are 
the  spherical  and  circular  error  probability  distributions  defined  earlier. 

(Note  that  Pg ( R )  is  the  probability  distribution  of  radial  error  in  3  dimensions, 
and  P(.(R)  is  the  probability  distribution  of  radial  error  in  2  dimensions.) 

3.o  Error  Probability  Relationships  to  Accuracy  Measures 

In  Figures  3.1  -  3.3,  plots  of  P^  versus  R  are  given,  where  R  has  been 
normalized  by  c^,  SEP,  and  crp,  respectively.  In  each  figure,  graphs  of  the 
special  cases  of  spherical  distribution  of  errors,  circular  degeneracy,  and 
linear  degeneracy  have  been  included.  In  Figures  3.4  -  3.6,  plots  of  P^ 
versus  R  are  given,  where  R  has  been  normalized  by  CEP,  and  aH,  respectively. 
In  each  figure,  graphs  of  the  special  cases  of  circular  distribution  of  errors 
and  linear  degeneracy  have  been  included. 

Observe  in  Figure  3.1  that  the  three  graphs  pass  through  the  point  P$  *  0.19875, 

R  ■  (as  they  must,  by  definition)  but  as  R  increases  the  three  graphs  diverge. 
Consequently,  is  not  a  very  good  indicator  of  the  value  of  R  at,  say, 

P  *  0.5  or  higher.  In  Figure  3.2  the  three  graphs  pass  through  Ps  »  0.5, 

R  *  SEP  (as  required)  but  diverge  considerably  at  P  =  0.75  or  higher.  In 
Figure  3.3,  there  is  no  single  point  (other  than  the  origin)  through  which 
the  three  graphs  all  pass.  However,  the  three  graphs  are  very  near  each  other 
in  the  vicinity  of  P$  *  0.77,  R  ■  1.2ap;  and  for  higher  values  of  P^,  the 
graphs  do  not  diverge  as  severely  as  those  in  Figures  3.1  and  3.2. 

An  obvious  conclusion  to  be  drawn  from  Figures  3.1  -  3.3  is  that  while  is 
a  "perfect"  indicator  of  R  for  P$  *  0.19875  and  SEP  is  a  "perfect"  indicator 
of  R  for  Ps  *  0.5,  1.2op  Is  a  "good"  Indicator  of  R  for  P^  4  0.77;  and 
furthermore,  ap  is  a  more  "robust"  Indicator  (i.e.,  less  sensitive  to  the 
shape  of  the  lc  ellipsoid)  of  the  value  of  R  for  P^  ■  0.75,  or  higher,  than 
either  SEP  or  o^.  This  conclusion,  coupled  with  the  fact  that  ap  is  easily 
computed  from  the  covariance  matrix,  makes  crp  a  useful  choice  of  accuracy 
measure  whenever  spherical  error  probabilities  of  .75  or  greater  are  of 
interest. 


Examination  of  Figures  3.4  -  3.6  leads  to  conclusions  for  2  dimensions  which 
are  analogous  to  those  for  3  dimensions.  In  particular,  aH  is  a  useful  choice 
of  accuracy  measure  whenever  circular  error  probabilities  of  0.75  or  greater 
are  of  interest. 

There  are  occasions  when  elliptical  error  probabilities  are  of  interest.  These 

probabilities  have  the  property  of  being  independent  of  the  shape  of  the  la 

ellipsoid  (or  ellipse)  as  long  as  the  distribution  is  non-degenerate.  Elliptical 

2 

error  probabilities  can  be  read  directly  from  tables  for  x  random  variables. 

(The  number  of  degrees  of  freedom  of  the  x  random  variable  equals  the  effective 
number  of  dimensions  of  the  error  vector  including  any  reduction  due  to  degen¬ 
eracy.)  Plots  of  elliptical  error  probabilities  for  cases  where  the  effective 
dimension  of  the  error  vector  Is  three,  two,  or  one,  respectively,  are  presented 
in  Figures  3.7  -  3.9. 


FIGURE  3.1  DISTRIBUTION  OF  3  DIMENSIONAL  RADIAL  ERROR  NORMALIZED  BY  SIGMA  S 


FIGURE  3.2  DISTRIBUTION  OF  3  DIMENSIONAL  RADIAL  ERROR  NORMALIZED  BY  SEP 


FIGURE  3.3  DISTRIBUTION  OF  1  DIMENSIONAL  RADIAL  ERROR  NORMALIZED  BY  SIGMA  P 


FIGURE  3.4  DISTRIBUTION  OF  2  DIMENSIONAL  RADIAL  ERROR  NORMALIZED  BY  SIGMA  C 


FIGURE  3.5  DISTRIBUTION  OF  2  DIMENSIONAL  RADIAL  ERROR  NORMALIZED  BY  CEP 


DISTRIBUTION  OF  2  DIMENSIONAL  RADIAL  ERROR  NORMALIZED  BY  SIGMA  II 
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FIGURE  3.8  EEP  FOR  2  DIMENSIONS  (EFFECTIVE) 


E  3.9  EEP  FOR  1  DIMENSION  (EFFECTIVE) 
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APPLICATION  CONSIDERATIONS 


The  purpose  of  this  section  is  to  provide  the  range  analyst  information  for 
developing  GPS  compatible  range  instrumentation  and  data  processing  systems. 
Range  requirements  vary  from  simple  coarse  position  data  to  sophisticated  Post 
test  data  processing  estimates  and  cannot  be  adequately  covered  in  this  docu¬ 
ment.  Subjects  discussed  in  this  section  are  major  concerns  that  became  appar 
ent  during  the  initial  phase  of  a  Tri service  GPS  Range  Applications  Study 
currently  in  progress.  The  data  processing  subjects  are  refinements  developed 
during  GPS  projects  at  the  Western  Space  and  Missile  Center,  Vandenberg  AFB, 
California.  These  subjects  were  researched  and  developed  by  several  members 
of  the  Systems  Performance  Analysis  Department;  consequently,  individualistic 
styles  are  apparent.  Recommendations  for  further  work  are  provided  in  some 
sections. 


4.1  Real-Time  versus  Post-Test 

As  in  virtually  all  test  and  training  range  applications  there  are  decidedly 
different  methods  and  uses  of  processing  data  in  real-time  ana  post-test  and 
the  same  is  true  for  GPS  data.  As  already  discussed  in  Section  2,  the  accur¬ 
acies  one  can  expect  from  GPS  differ  from  real-time  to  post-test.  Several 
other  sections  in  this  report  address  various  aspects  of  this  issue  (particu¬ 
larly  Sections  4.10  and  4.11);  the  purpose  of  this  section  is  to  give  a  gen¬ 
eral  discussion  of  the  overall  real-time  versus  post- test  considerations.  The 
slowly  varying  errors  referred  to  in  Section  2  are,  for  simplicity,  loosely 
called  biases  and  the  rapidly  varying  errors  we  refer  to  as  noise. 

Many  of  the  Triservices  GPS  test  and  training  applications  will  be  recognized 
as  real-time  scenarios.  Typical  examples  are: 

1.  Ballistic  missile  tracking  for  flight  safety. 

2.  Fleet,  unit  and  formation  relative  positioning. 

3.  Aircraft  or  missile  absolute  positioning  using  a  GPS 
aided  navigation  solution. 

4.  Ship,  vehicle,  or  troop  absolute  positioning  and  deployment. 

5.  Targeted  small  missiles  control  and  tracking. 

Post-test  processing  might  serve  to  merely  improve  upon  real  time  results, 
but  post-test  scenarios  also  form  a  separate  class  requiring  higher  accuracy 
and  larger  data  sets  than  real-time  scenarios.  Thus  post-test  reexamination 
and  analysis  of  GPS  data  can  apply  to  the  above  real-time  list. 

Typical  applications  more  strictly  classified  as  post-test  include: 

1.  Trajectory  reconstruction. 

2.  Instrumentation  accuracy  analysis. 

3.  Surveying. 

4.  Point  positioning  such  as  certain  impact  scoring  scenarios  which 
have  no  real  time  processing  requirements. 

In  addition,  test  scenarios  might  specify  a  "quick  look"  GPS  data  processing 
capability  which  may  generally  be  classified  as  post-test.  This  category 
may  allow  data  processing  at  sites  conveniently  removed  from  the  test  loca¬ 
tion  and  test  conditions. 


This  section  discusses  real-time  versus  post-test  accuracies.  The  analysis 
examines  the  individual  terms  of  the  error  budgets  presented  in  Section  2. 
Estimated  measurement  uncertainties  propagate  directly  into  expected  user 
state  uncertainties.  Two  factors  explain  the  increased  accuracies  attain¬ 
able  post-test:  (1)  improved  removal  of  bias  from  the  GPS  pseudo-range 
measurements,  and  (2)  optimal  processing  over  larger  data  sets.  These 
factors  depend  on  the  quality  of  the  data  and  the  nature  of  its  use,  as  dis¬ 
cussed  below. 

4.1.1  Bias  Removal 

Many  Triservices  GPS  scenarios  will  not  be  affected  by  Space  Segment  or 
atmospheric  delay  biases.  User  state  solutions  for  participants  in  small 
arena  differential  navigation  (discussed  in  Section  4.2)  will  be  offset  by 
a  common  bias.  Therefore,  they  are  not  included  in  either  the  real-time  or 
post-test  processing  error  budgets  for  these  scenarios.  Thus  the  following 
discussion  does  not  fully  pertain  to  relative  positioning  or  differential 
navigation  scenarios. 

The  real-time  versus  post-test  biases  of  chief  concern  appear  in  (1)  satel¬ 
lite  ephemerides,  (2)  satellite  clock  predictions,  and  (3)  atmospheric  delay 
corrections.  The  first  two  errors  are  closely  related  and  difficult  to  dis¬ 
tinguish  in  practice  (Ref  [16],  p  9).  They  are  the  values  of  satellite  state 
which  are  used  In  functional  equations  relating  user  to  satellite  slant 
ranges.  The  atmospheric  refraction  correction  is  an  estimated  bias  which 
is  subtracted  from  pseudo-range  measurements  to  give  corrected  values. 

The  GPS  Control  segment  tracks  GPS  satellites  over  extended  periods  of  time. 

The  error  model  of  the  Control  segment  determination  of  satellite  ephemerides 
Includes  both  bias  and  noise  terms.  Optimal  processing  techniques  provide 
"best  available"  estimates  of  the  true  satellite  ephemerides  which  are  in  turn 
periodically  uploaded  to  the  satellites  and  placed  into  satellite  data  formats. 
Thus  the  GPS  ephemerls  data  contains  residual  offsets  or  biases. 

The  Aerospace  Corporation  can  provide  corrected  ephemeris  estimates  having 
smaller  bias  errors  approximately  thirty  days  after  the  fact.  The  improved 
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reduction  of  bias  results  in  higher  post- test  user  accuracies.  In  GPS  Phase 
III,  the  real-time  ephemerides  will  eventually  be  known  as  precisely  as  Phase 
II  post-test  ephemerides  available  now.  Furthermore,  late  Phase  III  gains 
in  post-test  processing  due  to  the  removal  of  ephemeride  biases  will  not  be 
as  dramatic  as  in  current  Phase  II  except  for  replacement  satellites  (due  to 
Inexact  knowledge  of  their  ephemerides  for  some  time). 

Satellite  clock  and  timing  Information  is  extracted  from  the  GPS  data  format 
in  the  form  of  clock  drift  prediction  model  error  coefficients  and  satellite 
clock  offset  estimates  previously  uploaded  by  the  Control  Segment.  Thus, 
satellite  timing  uncertainties  increase  with  the  length  of  time  since  the 
last  update  and  can  cause  significant  errors  in  real-time.  Accuracy  gains 
in  post- test  data  pocessing  can  be  important  in  this  area.  Satellite  clock 
errors  do  affect  state  estimates,  but  errors  correlated  between  satellites 
can  result  In  a  compensating  value  of  user  clock  offset  such  that  a  common 
Control  Segment  uploaded  clock  error  may  not  degrade  user  coordinate  estim¬ 
ates  noticeably  (although  users  of  SPS  time  would  be  affected). 

Errors  introduced  by  refraction  corrections  are  due  to  errors  In  (1)  the 
estimated  tropospheric  delay  and  (2}  the  estimated  ionospheric  delay.  These 
corrections  may  be  determined  by  (1)  a  priori  constant  correction  and/or 

(2)  a  priori  atmospheric  modeling  with/without  actual  parametric  data.  The 
latter  method  is  generally  more  accurate.  Dual  frequency  P  code  users  have 
the  advantage  of  Improved  ionospheric  corrections  In  real-time  If  the  -L2 
Ionospheric  correction  is  used.  Furthermore,  pooled  data  concerning  atmos¬ 
pheric  conditions  can  increase  both  real-time  and  post-test  results.  Exten¬ 
sive  off  line  modeling  using  available  atmospheric  data  and  dual  frequency 
Ionospheric  delay  measurements  can  improve  the  final  atmospheric  corrections 
for  ionospheric  and  tropospheric  delay  for  post-test. 

4.1.2  Processing  Considerations 

Other  than  the  removal  of  bias  from  pseudo-range  measurements,  the  major  real¬ 
time  versus  post-test  consideration  is  the  method  of  processing  GPS  and  other 
measurements.  The  categories  of  concern  are  (1)  the  available  GPS  measurement 
history  as  compared  to  the  processed  data  span,  (2)  computer  algorithms,  and 

(3)  the  use  of  supporting  non  GPS  data.  In  this  section  GPS  measurements  are 


considered  to  be  simultaneous,  that  is,  pseudo-range  and  delta  pseudo-range 
measurements  are  taken  from  each  of  four  satellites  at  the  same  time  instant 
(thus  assuming  a  four  or  more  channel  receiver): 

There  are  several  different  methods  of  processing  the  measurements,  two  of 
which  are  (1)  compute  a  single  point  estimate  of  user  state  using  only  the 
measurements  taken  at  a  single  time  Instant,  or  (2)  compute  a  point  estimate 
of  user  state  using  a  subset  of  GPS  measurement  history  and  known  equations 
of  motion.  Both  methods  model  the  measurement  process  and  include  both  the 
slowly  and  rapidly  varying  terms  In  the  error  model. 

The  first  method  can  provide  better  post-test  accuracy  (than  real-time)  be¬ 
cause  some  of  the  GPS  error  budget  terms  are  smaller.  Here  the  emphasis  is 
on  the  better  satellite  state  Information  provided  post-test  as  explained 
above.  When  GPS  measurements  are  corrected  for  "known"  biases,  the  UERE  terms 
from  Section  2  are  lower,  and  these  smaller  expected  measurement  errors  lead 
directly  to  smaller  expected  user  state  errors. 

The  second  method  can  provide  more  accurate  results  post- test  than  real-time 
because  of  both  the  above  reason  and  the  constraints  of  the  equations  of 
motion,  assuming  these  are  adequately  specified.  Here  the  emphasis  is  on  the 
rejection  of  random  noise  in  the  GPS  measurements  by  user  software  processing. 
The  effect  of  user  receiver  noise  can  be  reduced,  and  receiver  biases  can  pos¬ 
sibly  be  Identified. 

A  full  error  analysis  Is  required  to  analytically  bound  the  noise  terms  of 
the  GPS  error  budget  on  an  Individual  basis  for  multi  sample  processing. 
Simplified  methods  exist  which  should  suffice  to  Indicate  a  reasonable  final 
UERE.  First,  an  empirical  approach  is  to  find  test  accuracy  results  for 
scenarios  similar  to  that  of  the  user.  If  the  POOP  terms  are  specified  the 
apparent  UERE  may  be  readily  factored  from  the  position  coordinate  residuals. 
Here  the  user  knows  a  correct  single  sample  UERE  (as  from  Section  2),  and 
previous  tests  indicate  a  noise  reduction  scale  factor  for  projected  scenar¬ 
ios.  A  second  method  is  to  classify  expected  user  dynamics  as  high,  medium 
or  low  in  terms  of  user  velocity,  acceleration  and  jerk.  This  simplified 
error  analysis  was  described  in  Ref  [2]  which  lists  expected  user  coordinate 
uncertainties.  Again  the  POOP  term  must  be  specified  in  order  to  factor  the 


4.1-4 


apparent  UERE,  and  the  UERE  is  needed  because  PDOP  is  widely  variable  (albeit 
calculable)  for  GPS  Phase  II  and  will  be  somewhat  variabl e  for  Phase  III. 

Stationary  users  and  those  with  smooth  dynamics  will  achieve  satisfactory 
noise  reduction  in  real-time  and  post-test.  For  stationary  users  Doppler 
smoothing  of  the  pseudo-range  data  can  be  accomplished  because  of  the  comple¬ 
menting  nature  of  delta  pseudo-range  measurements.  An  important  application 
is  accurate  surveying  in  a  matter  of  minutes  (Ref  [3]).  The  reduction  of 
noise  is  necessarily  limited  by  unmodeled  errors  and  multipath  effects,  but 
in  some  cases  noise  terms  in  the  applicable  error  budget  can  be  reduced  by 
{IT,  n  being  the  number  of  independent  samples.  Doppler  smoothing  over 
reasonable  intervals  of  continuous  simultaneous  tracking  has  been  claimed  to 
achieve  this  reduction,  both  for  receiver  and  atmospheric  correction  random 
noise.  For  these  reasons,  insofar  as  real-time  vs  post-test  processing  is 
concerned,  it  can  generally  be  stated  that  post-test  accuracies  will  be 
better  than  real-time  accuracies. 

The  data  rate  at  which  these  measurements  can  be  taken  and  stored  or  relayed 
to  recording  equipment  may  exceed  the  real-time  processing  capabilities  of 
on-board  navigation  computers  and  real-time  test  support  computers  on  ships, 
aircraft,  mobile  stations  or  even  permanent  land  based  sites.  Similarly,  up¬ 
dated  navigation  solutions  would  probably  not  be  available  after  every  GPS 
type  measurement.  Thus  post-test  processing  will  probably  involve  a  larger 
data  set,  and,  additionally,  optimal  batch  least  squares  or  Kalman  filter/ 
smoothers  will  operate  over  the  entire  data  set,  while  real-time  processing 
will  be  basically  sequential  with  no  or  limited  smoothing.  These  considera¬ 
tions  Involve  both  the  processed  data  span  and  computer  algorithms. 

Some  scenarios  might  include  non  GPS  data  such  as  vehicle  state  estimates  from 
radar  and  telemetry  data,  measurements  from  Inertial  devices,  altimeter  read¬ 
ings,  tide  tables,  or  survey  coordinates  of  participants.'  These  additional 
data  sets  improve  both  real-time  and  post- test  user  state  estimates,  or  con¬ 
versely,  GPS  data  augments  traditional  tracking  and  user  state  estimation 
schemes. 


4.2  Relative  versus  Absolute  and 

Small  Arena  versus  Global  Navigation  Issues 

4.2.1  Relative  vs  Absolute  Accuracies 

The  most  direct  method  of  GPS  navigation  is  to  utilize  four  or  more  pseudo¬ 
range  measurements  (from  suitable  SVs)  and  solve  for  user  position  and  the 
user  time  offset  relative  to  GPS  time.  The  user  position  is  solved  in  an 
earth  centered  coordinate  system.  This  method  is  called  absolute  or  global 
navigation.  The  non  absolute  navigation  means  either  that  one  user  set  is 
navigating  relative  to  another  user  set  via  a  data  link,  or  that  the  user 
has  been  at  a  given  point  at  some  prior  time  and  is  navigating  relative  to 
coordinates  measured  at  that  point.  Both  of  these  methods  are  referred  to 
as  differential  navigation  but  for  the  purpose  of  this  section  we  will  refer 
to  the  former  as  relative  and  the  latter  as  delta. 

The  GPS  system  errors  have  been  presented  in  Table  2.1.  These  system  errors 
will  affect  the  accuracy  of  the  absolute  navigation.  When  the  determination 
of  only  the  coordinate  differences  between  two  locations  is  the  main  concern, 
the  differential  navigation  techniques  can  provide  increased  accuracy.  Since 
the  bias  errors  caused  by  transmitter  imperfections,  ephemeris  and  clock  pre¬ 
diction,  and  ionospheric  delay  are  similar  at  both  the  user  location  and  the 
destination,  they  approximately  cancel.  These  errors  are  usually  the  larg¬ 
est  component  of  the  total  error.  The  receiver  bias,  which  is  user  depen¬ 
dent,  is  another  component  of  the  tota1  error.  The  receiver  bias  is  cancel¬ 
led  out  in  the  delta  navigation,  but  not  in  the  relative  navigation.  The 
gains  (one  sigma)  that  can  be  made  in  accuracy  in  relative  navigation  versus 
absolute  navigation  are  summarized  in  Table  4.2.1. 

An  illustration  of  the  effect  produced  by  the  differential  navigation  tech¬ 
niques  is  shown  in  Figure  4.2.1,  (Ref  [9]).  It  is  a  segment  of  the  11  January 
1980  mission  data,  using  Texas  Instrument's  High  Dynamic  User  Equipment  (HDUE). 
When  operating  in  the  normal  (absolute  navigation)  mode,  the  HDUE  provided 
position  estimates  that  had  root  mean  square  errors  about  20  meters  horizon¬ 
tally  and  near  40  meters  vertically.  When  HDUE  was  placed  in  the  differ¬ 
ential  mode,  these  errors  were  immediately  reduced  to  the  neighborhood  of 
5  meters  in  this  example. 
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TABLE  4.2.1  GPS  ERROR  3UDGET  (A8S0LUTE  AND  RELATIVE  NAVIGATION) 


Ionospheric  refraction  generates  range  errors  which  can  be  approximately  can¬ 
celled  out  in  the  relative  navigation  if  the  two  receivers  are  relatively 
near  one  another.  The  cancellation  will,  however,  be  gracefully  degraded 
over  space  and  time  due  to  the  translocation  of  ionospheric  correction. 

The  ionospheric  error  is  dependent  on  both  the  characteristics  of  the  iono¬ 
sphere  at  zenith  (at  a  given  time)  and  the  elevation  angle  to  the  satellite 
from  the  user.  The  ionospheric  delay  at  a  given  elevation  is  increased, 
relative  to  the  delay  for  a  ray  to  a  satellite  at  zenith,  by  the  obliquity 
factor.  The  obliquity  factor's  dependence  on  elevation  is  illustrated  in 
Figure  4.2.2  (Ref  [16],  p  34),  and  Figure  4.2.3  (Ref  [16],  p  34)  shows  typical 
measurements  of  ionospheric  delay  for  an  L-band  signal  received  at  vertical 
incidence.  The  mean  ionospheric  delay  (vertical  incidence)  is  about  ten  nsec 
at  nighttime,  while  it  can  be  as  large  as  50  nsec  during  daytime.  The  delays 
can  be  significantly  larger  in  regions  near  the  geomagnetic  equator  or  near 
the  poles,  particularly  during  the  magnetic  storm  periods.  As  shown  in  Figure 
4.2.2,  the  ionospheric  delay  can  be  three  times  the  values  given  above,  at  low 
elevation  angle.  Thus  the  degradation  of  the  cancellation  of  ionospheric  cor¬ 
rection  errors  over  space  and  time  must  be  considered  in  relative  navigation 
applications. 

The  degradation  due  to  translocation  of  Space  and  Control  Segment  errors  can 
be  bounded  in  a  rather  straightforward  fashion  for  both  relative  and  delta 
navigation.  The  case  for  relative  navigation  between  two  points  A  and  B  is 
illustrated  in  Figure  4.2.4.  Here  it  is  assumed  that  the  pseudo-range  errors 
are  monitored  at  a  known  location  A  and  applied  to  correct  the  pseudo-ranges 
measured  at  B.  The  relative  navigation  solution  is  degraded  only  by  the 
difference  in  pseudo-range  error  at  A  and  B.  Thus,  SV  clock  error  has  no 
effect  on  the  solution,  and  only  the  SV  ephemeris  error  which  is  orthogonal 
to  the  bisector  of  angle  ASVB  and  in  the  plane  ASVB  contributes  to  pseudo¬ 
range  error  differences  at  A  and  B. 

The  case  for  delta  navigation  is  illustrated  in  Figure  4.2.5.  Here  it  is 
assumed  that  the  pseudo-range  errors  are  monitored  at  a  known  initial  loc¬ 
ation  and  applied  to  correct  the  pseudo-ranges  measured  subsequently  as  a 
single  vehicle  moves  away  from  the  initial  location.  The  pseudo-range  error 
bound  shown  in  the  figure  is  based  on  the  first  harmonic  of  a  slowly  varying 
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ephemeris  error.  The  period  of  the  harmonic  is  where  T  is  the  SV  period. 

To  this  bound  on  |  a&[  must  be  added  the  term  C  Af/fQ,  where  f  is  the  nom¬ 
inal  value  of  SV  clock  frequency,  Af  is  the  magnitude  of  the  SV  clock  fre¬ 
quency  from  nominal,  and  C  is  the  speed  of  light.  Assuming  the  values 

i«Mxi  ■ 10  "eters 

and 

4^  *  5  X  10“12  , 

To 

which  are  larger  than  implied  by  the  GPS  error  budget,  it  follows  that 
[  aA  [  <_  0.007  meters/sec. 

4.2.2  Small  Arena  vs  Global  Navigation  Issues 

Small  Arena  applications  generally  relate  relative  positions  to  a  known  posi¬ 
tion.  These  may  involve  the  use  of  a  large  network  of  receivers  or  translators 
for  multiple  objects  of  interest.  Global  Navigation  applications  generally 
involve  the  use  of  receivers  and  the  availability  of  virtually  a  complete 
NAVSTAR  constellation  to  determine  absolute  position.  Both  Small  Arena  and 
Global  Navigation  applications  may  use  post-test  processing,  pseudo  satellites, 
and  other  methods  to  improve  the  solution's  accuracy. 

4. 2. 2.1  Small  Arena  Issues 

Most  Small  Arena  applications  involve  the  use  of  test  ranges  where  controlled 
conditions  and  some  post-test  processing  can  be  employed.  The  main  issue  for 
Small  Arena  applications  using  translators  or  transdigitizers  (discussed  in 
Section  4.7)  Is  spectral  bandwidth.  With  the  C/A  code  signal  this  involves 
the  use  of  *  2  MHz  per  object  of  interest  (with  dual  frequency  P  code  trans¬ 
lators  this  increases  to  *  40  MHz).  It  is  obvious  that  extreme  care  must  be 
used  to  avoid  jamming  and  intermodulation  problems. 

If  only  relative  position  is  desired  ionospheric  refraction  is  less  of  a 
problem  for  the  Small  Arena  as  all  receivers  can  observe  the  same  NAVSTAR 
satellites  with  similar  viewing  angles  thus  allowing  the  atmospheric  delay 


error  to  virtually  cancel.  It  should  be  noted  that  many  Small  Arena  appli¬ 
cations,  while  involving  the  relative  position,  require  that  a  control  point 
"absolute"  position.be  known.  The  absolute  position  of  a  fixed  central  point 
may  be  determined  from  geodetic  survey  references  with  a  root  mean  square 
accuracy  of  <_  1  meter.  This  level  of  absolute  accuracy  may  also  be 
achieved  with  a  GPS  receiver  with  good  POOP  and  post-test  refinement.  If  a 
mobile  control  point  is  required  (such  as  for  GPS-SMILS)  then  the  need  to 
determine  the  absolute  location  of  this  moving  point  virtually  dictates  the 
use  of  a  dual  frequency  P  code  GPS  receiver  with  post-test  refinement. 

Pseudo  satellites  can  be  used  to  augment  the  current  limited  NAVSTAR  con¬ 
stellation  whenever  practical.  Pseudo  satellites  should  improve  the  POOP, 
interval  of  coverage,  and  reduce  the  level  of  post- test  processing. 

4. 2. 2. 2  Global  Navigation  Issues 

The  use  of  GPS  receivers  for  Global  Navigation  generally  augment  other  sources 
of  navigation  data  such  as  inertial  systems,  LORAN,  and  other  navigation  data 
to  provide  absolute  position.  An  example  of  this  is  the  location  of  a  SLBM 
submarine  for  launch  coordinates.  The  nature  of  these  applications  limit  the 
amount  of  practical  post-test  processing  and  thus  limits  the  potential  accuracy 
to  *  10  meters  (one  sigma).  Most  Global  Navigation  applications  are  in  remote 
locations  where  the  availability  of  absolute  geodetic  fixes  (terrestial  bench¬ 
marks)  are  limited. 


4.3  Satellite  Constellations  Current  and  Plannee 

With  the  reduction  of  the  number  of  GPS  satellites  from  24  to  18,  various 
alternative  18-satellite  arrangements  were  studied.  The  objective  of  these 
studies  was  to  maximize  the  availability  of  the  system  to  users  despite  the 
reduction  in  the  number  of  satellites.  Of  the  many  alternatives  considered, 
the  current  baseline  Phase  III,  orbital  configuration  for  the  operational 
phase  employs  18  satellites  in  55-deg  inclined,  circular,  12  hour  orbits  to 
transmit  navigation  signals.  It  will  provide  continuous  three  dimensional 
global  coverage  by  placing  three  satellites  equally  spaced  (120  degrees  apart) 
in  each  of  six  orbit  planes.  These  six  orbit  planes  are  60  degrees  apart  in 
inertial  space,  as  depicted  in  Figure  4.3.1.  This  arrangement  of  18  satel¬ 
lites  has  been  found  to  be  best  as  far  as  accuracy  (good  GDOP)  and  uniform 
global  coverage  is  concerned.  An  alternate  constellation  employs  six  satel¬ 
lites  in  each  of  three  orbital  planes. 


With  only  four,  five  or  six  satellites  in  the  current  development  phase,  GPS 
provides  potential  users  a  few  hours  of  service  each  day.  Navigation  inter¬ 
vals  for  42  cities  around  the  world  are  tabulated  as  follows.  Table  4.3.1 
shows  the  interval  for  four  satellites.  As  the  fifth  and  sixth  satellites 
are  added,  the  periods  of  navigation  get  considerably  longer,  as  shown  in 
Tables  4.3.2  and  4.3.3  (Ref  [7]).  The  total  GPS  daily  coverage  based  on  the 
criterion:  PD0P<6,  is  shown  for  several  locations  of  interest  in  Table  4.3.4 
as  the  number  of  satellites  varies  from  6  to  18. 

The  primary  location  for  GPS  testing  with  the  current  constellation  is  the 
Army  Proving  Grounds  in  Yuma,  Arizona.  The  accuracy  of  3  Dimensional  naviga¬ 
tion  at  Yuma  with  four  Space  Vehicles  (SV)  is  indicated  by  the  solid  line  on 

Figure  4.3.2  (Ref  [7]).  The  SEP  (Spherical  Error  Probable),  defined  in 

Section  3,  during  the  two-hour  period  is  shown.  Near  the  middle  of  the  navi¬ 
gation  interval  when  the  satellite  geometry  is  more  favorable,  the  SEP  is 
about  10  m.  With  the  addition  of  a  fifth  satellite,  the  available  test  time 
is  extended  about  90  minutes,  as  shown  by  the  dotted  line  in  Figure  4.3.2. 

The  accuracy  during  this  additional  time  period,  however,  is  degraded,  be¬ 
cause  of  the  poor  satellite  geometry  at  Yuma.  Similar  figures  for  Grand 

8ahama  Island  and  Cold  Lake,  Canada  are  shown  in  Figures  4.3.3  and  4.3.4. 
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TABLE  4.3.1  NAVIGATION  AVAILABILITY  WITH  FOUR  SATELLITES 


GMT  GMT  HRS  GMT  GMT  HRS  Total 

User  Location  Start  End  Duration  Start  2nd  Duration  Duration  ihrsi 


Acapulco,  Mexico 

14.3 

13.3 

1.0 

1.0 

Anchor aga,  Alaska 

14.3 

16.0 

1.3 

1.3 

Ankara,  Turkey 

1.17 

3.3 

1.3 

1.3 

Sruaaala,  Salglua 

2.3 

3.3 

1.0 

1.0 

Suanoa  Air as,  Argentina 

19.3 

21.0 

1.3 

1.3 

Caap  Parka,  California 

13.3 

13.3 

2.3 

2.3 

Calcucta,  India 

2.3 

4.0 

1.3 

1.5 

Capa  Kennedy,  Florida 

14.3 

16.2 

1.7 

1.7 

Capa  Town,  South  Africa 

20.3 

22.2 

1.7 

1.7 

Caracas,  Vanazuala 

13.3 

16.1 

0.6 

0.6 

Charleston,  Souch  Carolina 

14.2 

16.2 

2.0 

2.0 

Christchurch,  New  Zealand 

3.0 

10.3 

2.3 

2.3 

Cold  Lake,  Canada 

14.2 

16.2 

2.Q 

2.0 

Diego  Garcia 

Hoc  available 

Eglln  Air  Force  Sasa 

14.2 

16.2 

2.0 

2.0 

Farn borough.  United  Kingdom 

2.3 

3.3 

1.0 

1.0 

Fort  Monaouth,  New  Jersey 

14.3 

16.2 

1.7 

1.7 

Fortune,  North  Dakota 

14.2 

16.3 

2.0 

2.0 

14.2 

16.2 

2.0 

2.0 

Grand  Bahaaa  Island 

14.3 

16.2 

1.7 

1.7 

Cuaa 

Hot  available 

Honolulu,  Hawaii 

13.2 

14.3 

1.1 

l.l 

Komis 

1.2 

3.3 

2.0 

2.0 

Laa  Palaas,  Canary  lalanda 

Hoe  available 

Una,  Peru 

Noe  available 

New  Haapshlra  (SCF) 

14.7 

16.2 

1.3 

1.5 

New  London,  Connecticut 

14.7 

16.2 

1.3 

1.5 

Pago  Pago,  American  Samoa 

9.0 

10.2 

1.2 

1.2 

Panama  Canal 

13.2 

13.3 

0.3 

0.5 

tlyadh,  Saudi  Arabia 

1.2 

3.0 

1.3 

1.3 

Soma,  Italy 

2.0 

3.0 

1.0 

1.0 

Seychallaa  Island 

Not  available 

Sldnay,  Australia 

7.3 

9.3 

2.0 

2.0 

Stockholm,  Swadan 

2.3 

4.0 

1.3 

1.5 

Taipei,  Formosa 

3.3 

4.2 

0.7 

0.7 

Tananarive,  Malagasy  Republic 

Hoc  available 

Tel  Aviv,  Israel 

1.2 

3.0 

1.3 

1.3 

Thule,  Creanlaad 

Noe  available 

Tokyo,  Japan 

Noe  aval la bio 

Tremao,  Norway 

2.7 

3.3 

1.0 

1.0 

Vandanberg  Air  Force  Base 

13.3 

13.3 

2.0 

2.0 

Yuma,  Arizona 

13.3 

13.3 

2.3 

2.3 

Hocai  Navigation  lntarvala  o f  lass  chan  O.J  hr  vara  nee  lacludad  bacauaa 
•  they  ara  not  consldarad  useful. 


TA8LE  4.3.2  NAVIGATION  AVAILABILITY  WITH  FIVE  SATELLITE  CONSTELLATION 


User  Location 

GMT 

Start 

GMT 

End 

HRS 

Ouratlon 

GMT 

Start 

GMT 

End 

HRS 

Durtcli 

Acapulco,  Mexico 

12.0 

13.3 

2.3 

Anchorage,  Alaska 

13.2 

16.0 

2.8 

Ankara,  turkay 

0.3 

3.3 

2.8 

Brussals,  Selgiua 

1.3 

3.6 

2.1 

luanos  Alras,  Argentina 

18.2 

21.1 

2.8 

Caap  Parks,  California 

13.0 

13.9 

2.9 

Calcucca,  India 

1.3 

4.2 

2.7 

Capa  Kennedy,  Florida 

14 .0 

17.2 

3.2 

Capa  Town,  South  Africa 

20.3 

23.3 

3.0 

Caracas,  Vanazuala 

14.9 

18.8 

3.8 

Charleston,  South  Carolina 

14.2 

17.2 

2.9 

Christchurch,  New  Zealand 

7.3 

10.7 

3.2 

Cold  Lake,  Canada 

13.0 

16.7 

3.7 

Olago  Garcia 

1.7 

2.3 

0.5 

Eglm  Air  force  Base 

14.0 

16.7 

2.7 

Farnborough.  United  king don 

1.7 

3.5 

1.3 

Fort  Monaouth,  New  Jersey 

14.2 

17.3 

3.3 

Fortune,  North  Dakota 

13.2 

16.7 

3.5 

Galveston,  Texas 

13.7 

16.5 

2.7 

Grand  Bahama  Island 

14.0 

17.2 

3.2 

Cuaa 

Hoe  available 

Honolulu,  Hawaii 

12.0 

14.4 

2.4 

Horauz 

1.0 

3.3 

2.3 

Laa  Fainas,  Canary  Islands 

Noe  available 

Lina,  Peru 

18.6 

19.1 

0.5 

Maw  Haapshlre  (SCF) 

14.6 

17.5 

2.9 

Maw  London,  Connecticut 

14.6 

17.3 

2.9 

Pago  Pago,  Ana rl can  Saaoa 

7.3 

8.3 

0.8 

8.3 

10.4 

1.6 

Panama  Canal 

13.3 

14.3 

0.8 

14.9 

15.9 

1.0 

Riyadh,  Saudi  Arabia 

0.3 

3.1 

2.6 

Rone,  Italy 

1.0 

3.1 

2.1 

Seychelles  Island 

1.2 

2.0 

0.8 

Sidney,  Australia 

7.2 

9.6 

2.3 

Stockholm,  Sweden 

1.5 

4.0 

2.5 

Taipei,  Foraosa 

3.2 

6.3 

3.0 

Tananarive,  Malagasy  Republic 

22.0 

23.3 

1.3 

Tel  Aviv,  Israel 

0.2 

3.1 

2.8 

Thule,  Greenland 

13.7 

14.4 

0.6 

15.2 

17.3 

2.0 

Tokyo,  Japan 

3.2 

3.0 

1.3 

Tronso,  Norway 

1.5 

2.4 

0.9 

2.6 

3.3 

1.2 

Vandanbarg  Air  Force  Sase 

13.0 

15.6 

2.6 

Tuns,  Arizona 

13.5 

15.9 

2.4 

NOTE:  Navigation  Intervals  of 

lass  chan  0. 3 

hr  ware 

not  Included 

because 

they  are  not  considered  useful. 


TABLE  4.3.3 


NAVIGATION  AVAILABILITY  WITH  SIX  SATELLITE  CONSTELLATION 


User  Location 

GMT 

Start 

GMT 

End 

HRS 

Duration 

GMT 

Start 

GMT 

End 

HRS 

Duration 

Total 

Duration  i  HRS ) 

Acapulco,  Moxlco 

13.0 

15.6 

2.6 

16.1 

18.3 

2.2 

«a  •  » 

Anchorage,  Alaska 

0.2 

3.2 

1.0 

4.5 

5.0 

0.5 

13.2 

16.0 

2.3 

16.6 

17.8 

1.2 

5.» 

Ankara,  Turkey 

0.3 

5.0 

4.5 

*.5 

Brussels,  Belgium 

1.4 

3.9 

2.5 

2.5 

Buenos  Alras,  Argentina 

13.2 

23.0 

4.3 

i.d 

Camp  Parks,  California 

13.0 

17.5 

4.5 

..5 

Calcutta,  India 

1.3 

7.3 

5.3 

5.3 

Caps  Kennedy,  Florida 

14.0 

13.3 

4.3 

Capa  Town,  South  Africa 

20.3 

0.5 

4.0 

..0 

Caracas,  Venezuela 

14.9 

13.3 

3.3 

3.3 

Charleston,  Souch  Carolina 

14.2 

13.3 

4.5 

-.5 

Christchurch,  haw  Zeeland 

7.3 

12.5 

5.0 

5.0 

Cold  Lake,  Canada 

13.0 

18.3 

3.3 

5.3 

Olego  Carcia 

1.7 

2.3 

0.5 

5.0 

5.5 

0.3 

1.0 

Eglin  Air  Force  Base 

14.0 

18.3 

4.3 

4.3 

Famborough,  United  Kingdom 

1.7 

3.3 

2.1 

2.1 

Fort  Monmouth,  New  Jersey 

14.2 

13.3 

4.0 

4.0 

Fortune,  North  Qekoca 

13.2 

18.3 

5.0 

5.0 

Calvescon,  Texas 

13.7 

19.0 

5.3 

5.3 

Crand  Bahama  Island 

14.0 

18.8 

4.8 

4.3 

Cuam 

Not  available 

Honolulu,  Hawaii 

12.0 

16.3 

4.3 

4.3 

Hormuz 

1.0 

3.8 

4.3 

4.3 

Las  Palmas,  Canary  Islands 

rs 

3.3 

1.0 

1.0 

Lima,  Peru 

18.3 

22.1 

3.6 

3.6 

New  Hampshire  (SCF) 

14.6 

18.0 

3.4 

3.4 

New  London,  Connecticut 

14.6 

18.0 

3.4 

3.4 

Pago  Pago,  American  Samoa 

7.3 

8.3 

0.3 

3.9 

10.4 

1.6 

10.3 

13.3 

3.0 

5.3 

Panama  Canal 

13.3 

14.3 

0.8 

14.9 

13.9 

1.0 

17.1 

20.3 

3.2 

5.0 

Riyadh,  Saudi  Arabia 

0.3 

5.8 

5.3 

5.3 

Rome,  Icaly 

1.0 

4.0 

3.0 

3.0 

Seychelles  Island 

1.2 

2.0 

0.3 

4.7 

5.5 

1.8 

1.5 

Sidney,  Australia 

7.2 

11.5 

4.3 

4.3 

Stockholm,  Sweden 

1.5 

4.0 

2.5 

4.6 

5.5 

0.9 

14.5 

15.2 

0.7 

4.1 

Taipei,  Formosa 

3.2 

6.3 

3.0 

3.0 

Tananarive,  Malagasy  Republic 

22.0 

23.3 

1.3 

1.3 

Tel  Aviv,  Israel 

0.2 

5.0 

4.3 

*.3 

Thule,  Greenland 

2.0 

3.3 

1.3 

3.7 

4.3 

1.0 

13.7 

17.3 

3.5 

5.3 

Tokyo,  Japan 

3.2 

5.0 

1.3 

1.3 

T corns  a,  Norway 

1.5 

4.4 

2.9 

4.6 

5.5 

0.9 

14.2 

15.2 

1.0 

15.7 

17.0 

1.3 

6.1 

Vandenberg  Air  Force  Base 

13.0 

17.5 

4.5 

4.5 

Tuna,  Arizona 

13.3 

17.8 

4.3 

4.3 

NOTE;  Navigation  Intervals  of 

lass  :hen  0.5 

hr  were 

noc  Included 

because 

they  are  not  considered  useful. 
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TABLE  4.3.4  TOTAL  GPS  DAILY  COVERAGE 


THREE  DIMENSIONAL  ERROR  IN  METERS.  SEP 


10  11  12  13  14  IS  16  17  18  19  20 

TIME  IN  HOURS,  GMT 

FIGURE  4.3,2  THREE-DIMENSIONAL  NAVIGATION  AT  YUMA 


FIGURE  4.3.3  THREE-DIMENSIONAL  NAVIGATION  AT  GRAND  3AHAMA 


When  the  altitude  at  the  user  location  is  known,  2  Dimensional  (horizontal) 
navigation  can  be  performed  (see  Section  4.4  for  an  analysis  of  this  method). 

It  is  noticed  that  the  time  period  available  for  testing  2  Dimensional  naviga¬ 
tion  at  Yuma  is  extended  considerably  as  shown  in  Figure  4.3.5.  This  is  due 
to  the  fact  that  this  mode  of  navigation  makes  use  of  some  geodetic  height 
information  thus  necessitating  the  use  of  only  three  satellites  (see  Section 
4.4). 

The  user  position  errors  in  GPS  navigation  are  determined  by  two  factors: 
the  magnitude  of  the  ranging  errors  and  the  geometry  relative  to  the  selected 
satellites.  The  ranging  errors  are  statistically  described  by  la  UERE,  while 
the  effect  of  the  geometry  is  expressed  by  geometric  dilution  of  precision  (GDOP) 
parameters.  These  parameters  include  POOP,  HDOP,  VDOP  and  TDOP,  as  discussed 
in  Section  3. 

Small  values  of  POOP  and  HDOP  indicate  good  geometry  of  the  selected  satellites 
and,  correspondingly,  small  errors  in  user  position  estimation.  Geometrical 
dilution  can  be  visualized  by  the  volume  of  a  special  tetrahedron  as  shown  in 
Figure  4.3.6.  Let  P  denote  the  user  position  and  let  S  be  a  unit  sphere  cen¬ 
tered  at  P,  The  vectors  from  the  user  position  to  each  satellite  will  inter¬ 
sect  the  unit  sphere  at  four  points.  The  tetrahedron,  formed  by  connecting 
these  points  has  a  volume  V.  As  the  volume  of  the  tetrahedron  becomes  larger, 
the  geometric  strength  of  the  solution  improves  accordingly.  It  is  noted  that 
the  volume  is  maximized  when  one  satellite  is  at  the  user's  zenith,  and  the  other 
three  are  separated  by  1 20°  and  are  as  low  on  the  horizon  as  permitted  by  the 
user's  antenna  elevation  angle. 

GDOP  parameters  depend  only  on  the  orientation  of  the  selected  satellites  re¬ 
lative  to  the  user.  They  are  functions  of  the  user  location  on  earth,  and 
change  during  the  course  of  a  day,  due  to  earth  rotation  and  satellite  motion. 
There  is  no  analytical  method  to  determine  the  statistical  behavior  of  GDOP 
parameters.  The  GDOP  parameters  are  generally  given  in  terms  of  percentiles, 
such  as  the  fiftieth  and  ninetieth  percentile  by  conducting  a  Monte  Carlo  sim¬ 
ulation  covering  all  latitude,  longitude  and  time.  The  lack  of  a  mathematical 
expression  for  the  geographical  distribution  of  GDOP  parameters  has  prevent  d 
thu  establishment  of  a  mathematical  relationship  between  percentile  of  GDOP 
and  that  of  user  navigation  errors. 


A  computer  program  to  perform  Monte  Carlo  simulation  has  been  developed  (Ref 
[13],  [23]).  It  divides  the  earth  into  regions  bounded  by  increments  of  lati¬ 
tude  and  longitude  and  computes  the  GDOP  parameters  for  a  user  at  the  center 
of  these  regions.  Since  the  GDOP  value  varies  during  the  course  of  a  day, 
the  computations  are  repeated  at  successive  time  increments  for  each  region. 
The  results  are  given  in  Table  4.3.5  (Ref  [16],  p  1])  together  with  corres¬ 
ponding  values  of  the  user  error  based  on  1c  UERE  (11.8  to  20.7  feet).  Note 
that  the  results  are  based  on  the  original  24  satellite  constellations.  Cor¬ 
responding  results  based  on  4,  12  or  18  satellites  are  not  available  at  this 
time. 


TABLE  4.3.5  ANTICIPATED  WORLDWIDE  USER  POSITION  ERROR  DISTRIBUTION 


koroonial 

Vorticai 

Tim* 

HDOP 

Umt  Er¬ 
ror  Pa¬ 
ramour 

Uwr  Error 
(Mount 

Uaor  Er¬ 
ror  (Foot) 

VDOP 

|g! 

Uoor  Error 
(Mount 

L'aor  Er¬ 
ror  i  Froil 

TDOP 

l'aor  Er¬ 
ror  Pa¬ 
rameter 

I'ser  Er¬ 
ror  i  N*. 
noae- 
condai 

50th  percen¬ 
tile 

1.39 

1.15 

4. 1-7.2 

14-24 

1.99 

1.39 

5.0-8.8 

16-29 

1.05 

0.73 

8-15 

rms 

1.44 

1.45 

5.2-9. 1 

17-30 

2.16 

2.21 

8.0-13.9 

26-46 

1.21 

1.22 

14-25 

90th  percen¬ 
tile 

1.71 

2.19 

7.9-13.8 

26-45 

2.80 

3.57 

12.9-22.5 

42-74 

1.76 

i 

1.96 

23-40 
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2  versus  3  Dimensional  Navigation 


The  NAVSTAR/GPS  system  is  designed  to  provide  accurate  3  dimensional  position 
to  a  user  anywhera  on  or  near  the  earth.  When  the  entire  eighteen  satellite 
constellation  is  complete  this  3  dimensional  navigation  will  be  possible  at 
almost  all  times.  However,  until  the  constellation  is  complete  many  areas 
will  have  only  a  fraction  of  each  day  when  the  requisite  four  satellites 
(necessary  for  the  3  dimensional  navigation)  are  suitably  located  for  use. 
When  adequate  altitude  information  is  available  a  user  can  use  only  three 
satellites  to  perform  a  2  dimensional  (horizontal)  navigation.  The  uncer¬ 
tainty  in  the  altitude  information,  of  course,  directly  influences  the  accur¬ 
acy  of  this  2  dimensional  navigation.  In  this  section,  3D  navigation  means 
that  a  priori  altitude  information  is  optimally  combined  with  GPS  measure¬ 
ments  to  construct  a  3D  solution.  On  the  other  hand,  2D  navigation  means 
that  a  priori  altitude  information  is  assumed  perfect  and  combined  with  GPS 
measurements  to  construct  only  the  2D  (horizontal)  part  of  the  solution. 

The  general  GPS  3  dimensional  navigation  has  been  discussed  in  Section  3, 
where  no  a  priori  estimate  uncertainties  are  considered.  The  purpose  of  this 
subsection  is  to  analyze  accuracy  measures  for  both  2  and  3  dimensional 
navigation  when  a  priori  altitude  estimates  are  available.  When  an  a  priori 
estimate  of  the  geodetic  height  is  known,  with  a  given  variance,  both  HOOP 
and  POOP  can  be  improved  if  the  a  priori  estimate  is  properly  weighted  using 
the  a  priori  variance. 

Let  an  orthogonal  coordinate  system  be  selected  such  that  the  x-y  plane  is 
horizontal  and  the  z-axis  is  vertical.  The  basic  equations  for  GPS  naviga¬ 
tion  systems  are 

+  ct  *  Rif  i  »  1,2, 3, 4,  (4.4.1) 

where  x,  y,  z  and  t  are  user  position  and  unknown  clock  bias;  c  is  the  speed 
of  light;  x^ ,  y^ ,  z^  are  the  1th  SV  position,  and  R^  is  the  pseudo-range 
measurement  to  the  ith  SV. 


Since  equation  (4.4.1)  is  nonlinear,  a  linearization  using  first-order  Taylor 
expansion  is  applied.  Let  x,  y,  z  and  t  be  the  nominal  values  of  x,  y,  z  and 
t,  and  let  dx,  <sy,  dz  and  dt  be  the  corrections  to  these  nominal  values,  i.e., 
dx  *  x-x,  dy  3  y-y,  dz  3  z-z  and  dt  3  t-t.  The  nominal  value  of  the  pseudo¬ 
range  measurement  from  the  ith.  SV  to  user  is  given  by 


Ri  -/(x-x^2  +  (y-y.,)2  +  (z-z^2  +  ct 


(4.4.2) 


The  offset  of  the  pseudo-range  measurement  from  its  nominal  value  is 


4Ri  *  Ri  -  Ri 


Expansion  of  (4.4.1)  to  first  order  about  the  nominal  leads  to 


(x-x.)dx  +  (y-yi  )sy  +  (z-z^dz 


-  ct 


+  cdt  3  aR^  ,  i  3  1 ,  2,  3,  4. 


These  four  equations  can  be  written  in  a  matrix  form  as 


He  3  a R 


where 


,  aR 


(4.4.3) 


*  C(x-x.)  ,  (y- y1)  ,  (z-z.)  ]/  (R.  -  ct) 

is  the  unit  vector  from  itli  SV  to  user. 


Let  PQ  denote  the  covariance  matrix  of  the  nominal  values  of  x,  y,  z  and  t. 
Since  the  geodetic  height  of  the  user  position  is  known,  it  is  assumed  that 
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z  is  Gaussian  distributed  with  mean  z  and  variance  a The  variances  of 
x,  y  and  t  are  assumed  infinity  .  Thus 


rO  0  0  0 

0  0  0  « 

where  zero  correlation  is  assumed. 

The  minimum-variance  estimates  of  6x,  <5y,  <5Z  and  c<$t  are  given  by 

3  *  [P"1  +  ht  m"1  h]"1  ht  M"1  aR 

where  M  is  the  covariance  matrix  of  the  pseudo-range  measurements.  It  is 

assumed  that  the  four  UEREs  are  statistically  independent  and  a0  is  the  RSS 

2  K 

UERE  common  to  all  four  measurements.  Thus  M  *  a  I  I.,  where  I  is  a  mxm 

A  k  4  m 

identity  matrix.  The  covariance  matrix  of  3  is  given  by 

P  *  CPQ"1  +  HT  M"1  H]"1 


X 

axy 

°xz 

C5xt 

yx 

2 

ay 

ayz 

cayt 

zx 

ffzy 

2 

az 

cazt 

tx 

caty 

cotz 

C2at' 

In  Section  3,  POOP  is  defined  by  (ax2  +  ay2  +  a22)1//2/ffR  and  HOOP  is  defined 

by  (<jx2  +  a y2)1^2/^^.  Both  POOP  and  HDOP  are  functions  of  H,  M  and  PQ. 

Figures  4.4.1,  (a)-(d),  show  therelationship  of  POOP,  HOOP  and  (az/aR)  with 
various  SV  arrangements .  The  geometry  of  the  SVs  is  arbitrary,  and  is  sel¬ 
ected  for  demonstration  only. 


In  the  2  dimensional  navigation  case,  the  a  priori  altitude  estimate  is 
assumed  perfect  and  only  the  horizontal  coordinates  of  the  user  position  are 
estimated.  Hence,  it  is  of  interest  to  investigate  the  accuracy  of  the  2  dim< 
ensional  solution  as  a  function  of  the  unmodeled  altitude  variance. 

It  follows  from  (4.4.3)  that: 

(x-x.)  <sx  +  (y-y  • )  <5y 

- ; - : -  +  c5t  3  e.  ,  i  3  1,  2,  3,  4,  (4.4.4) 

Ri  -  ct  1 

where  ei  3  aRi  -  (z-z^Mz-zJ/tR.-ct). 

The  four  equations  in  (4.4.4)  can  be  written  in  a  matrix  form  as 

He  3  e  , 

where 


and 

e/  *  [(x-x,).  (y-y, ) 3  /  (R,  -  ct). 

Note  that  E(e)  3  0,  and  E(seT)  s  ^  I  +  ^  ffT  where  f  3  Cf , ,  f2,  f3,  f4]T 
with  f.  3  (z  -  z . )/{ R.  -  ct). 


>  • , 

s 


s 

1 

i 

L*' 

s 

"i 


P 


The  least-square  estimate  of  3  is  3  3  (HT  H)"1  H^e,  and  the  covariance 
matrix  of  ts  is  given  by 

c(3)  3  [HT  H]*1  HT  (cjjl  +  ~a\  ffT)  H  [HT  H]"1 

or 

c(3)  3  a j*  [HT  H]*1  +  a2z  [HT  H]"1  HT  ffT  H  [HT  H]'1  . 

2  “2 

Here,  the  HDOP  is  a  function  of  H,  f,  aR  and  Figures  4.4.2,  (a)-(d), 
show  the  relationship  of  HDOP  and  (<?2/aR)  with  various  SVs  geometry.  It 
should  be  noted  that  three  SVs  are  sufficient  to  perform  2  dimensional  navi¬ 
gation,  for  only  x,  y,  and  t  are  to  be  estimated.  The  relationship  of  HDOP 
and  (<JZ/<JR)  for  3  SVs  is  also  included  in  Figures  4.4.2,  (a)-(d). 

As  an  aid  in  quantifying  the  ideas  and  results  of  this  section,  it  is  useful 
to  introduce  the  notion  of  Constellation  Value  (CV).  For  a  given  navigation 
accuracy  criterion  and  a  given  SV  constellation,  CV  is  defined  as  the  prob¬ 
ability  of  satisfying  the  criterion,  as  determined  by  the  fraction  of  each 
day  for  which  the  criterion  is  satisfied  averaged  over  the  entire  globe. 

It  is  of  interest  to  examine  the  behaviour  of  CV  during  the  GPS  constellation 
buildup,  and  this  is  illustrated  for  three  cases  of  interest  in  Figure  4.4.3. 
The  three  cases  are: 

1.  3D  navigation  with  PD0P<6 

2.  2D  navigation  with  HD0P<4 

3.  2D  navigation  with  HD0P<4  assuming  perfect  altitude. 

In  the  first  two  cases,  no  a  priori  altitude  information  is  available,  while 

in  the  third  case  the  a  priori  (or  independent)  measure  of  altitude  is  assum¬ 
ed  perfect. 

The  figure  shows  that  the  CV  is  not  much  improved  by  specifying  only  a  2D 
navigation  criterion  unless  an  independent  measure  of  altitude  is  available, 
in  which  event,  the  CV  improves  considerably.  For  example,  at  the  time  when 
the  constellation  contains  12  satellites,  with  no  independent  altitude  in¬ 
formation  the  CV  for  3D  navigation  (PD0P<6)  is  *60%  and  the  CV  for  2D  navi- 


t 
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gation  (HD0P<4)  is  *70$,  but  with  perfect  altitude,  the  CV  for  2D  navigation 
(H00P<4)  is  *90%. 

Thus,  during  constellation  buildup,  augmentation  of  GPS  navigation  solutions 
with  independent  altitude  measurements  can  significantly  extend  the  daily 
window  of  operability  with  GPS. 
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FIGURE  4.4.2(d)  RELATIONSHIP  OF  HOOP  WITH  A  PRIORI  ALTITUDE  UNCERTAINTY 
IN  2  DIMENSIONAL  NAVIGATION. 


Satellite  configuration  (AZ.EL),  (0*,30*)f  (30*. 40*),  (45*. 50*) ,  (60*. 60*) 


Courtesy  of  JPO  and  Aerospace  Corporation 
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Pseudo-Satellites 


Ground  stations  that  transmit  signals  compatible  with  GPS  satellite  signals 
are  referred  to  as  pseudo-satellites.  An  example  of  their  use  is  at  the 
Yuma  Proving  Grounds  using  a  configuration  called  Inverted  GPS  (Ref  [64]}. 
Their  purpose  is  to  supplement  or  replace  the  need  for  GPS  satellite  signals 
during  demonstration  tests  of  GPS  concepts  and  User  Equipment.  Their  use  in 
Range  instrumentation  systems  would  have  two  major  advantages.  First,  test 
scheduling  would  be  Independent  of  satellite  coverage  and,  second,  they  could 
provide  Improvements  In  data  results  due  to  geometrical  considerations  and/or 
satellite  signal  reception  problems.  Command,  control  and  comnunlcatlons  (C3) 
functions  can  be  combined  with  the  pseudo- satellite  uplink.  For  some  ranges. 
It  might  be  necessary  to  operate  pseudo- satellites  even  after  GPS  is  opera¬ 
tional  with  the  full  18  satellite  array. 

Two  problems  must  be  considered  for  use  of  pseudo-satellites.  The  first 
problem  concerns  situations  when  test  vehicles  would  operate  very  near  the 
pseudo-satellite  and  Its  signal  level  would  be  sufficiently  higher  than  the 
satellites  or  other  pseudo-satellites  to  prevent  their  acquisition  or  track. 

If  there  is  only  one  test  vehicle,  this  might  be  overcome  by  adjusting  the 
transmitter  power  using  various  control  schemes.  If  all  participants  are 
sufficiently  far  from  the  pseudo- satellite,  one  power  level  might  be  found 
that  would  be  acceptable.  Conversely,  the  test  vehicle  could  control  the 
received  signal  level  with  a  directional  antenna  array. 

In  practice,  the  near-far  problem  has  not  been  a  serious  constraint  at  the 
Yuma  Proving  Grounds  Inverted  GPS  Range.  An  acceptable  transmitter  power 
Is  selected  and  GPS  receivers  on  test  vehicles  operate  satisfactorily  except 
for  acquisition  of  a  low  level  signal  when  near  a  pseudo-satellite.  Better 
definition  of  the  constraints  for  various  receiver  and  translator  configur¬ 
ations  and  signal  level  control  techniques  must  be  developed  before  evaluat¬ 
ing  use  of  GPS  technology  In  Range  Instrumentation.  The  range  of  signal 
levels  acceptable  during  various  receiver  and  translator  configurations  oper¬ 
ation  must  be  determined  analytically  and  verified  with  experiments  for  rep¬ 
resentative  cases.  Then  pseudo- satellite  signal  level  control  require¬ 
ments  can  be  defined  and  suitable  configurations  developed  for  particular 


range  situations.  Figure  4.5-1  shows  a  missile  tracking  configuration 
being  considered  at  WSMC  in  which  the  pseudo-satellite  transmitter  power 
level  could  be  controlled  in  an  open  loop  manner  using  a  priori  trajectory 
data  or  In  a  closed  loop  manner  using  telemetered  GPS  receiver  signal  level 
data  or  it  may  not  be  necessary  to  vary  the  level. 

For  some  situations,  use  of  a  different  frequency  might  be  advantageous.  For 
Instance,  if  the  test  vehicle  receives  L1  signals  from  the  GPS  satellites, 
the  pseudo-satellite  could  be  operated  on  l^,  I3  or  another  frequency.  The 
advantage  Is  that  the  pseudo-satellite  signal  level  problem  is  greatly  sim¬ 
plified.  The  disadvantage  is  that  the  receiver  has  to  be  able  to  receive 
on  two  frequencies.  The  downlink  bandwidth  problem  can  be  minimized  by 
overlaying  the  received  signals  if  the  levels  are  compatible.  Obviously, 
the  test  vehicle  receiver  becomes  more  complicated. 

Use  of  C/A  code  would  be  adequate  and  even  preferable  for  many  Range  appli¬ 
cations.  The  pseudo-satellite  uplink  bandwidth  requirements  using  C/A  code 
are  one  tenth  those  for  P  Code.  The  L-j  C/A  code  signal  level  is  3  dB  higher 
than  the  P  Code  level  transmitted  by  the  GPS  satellites. 

The  second  problem  that  must  be  considered  with  pseudo- satellites  is  that 
their  data  Is  not  used  in  a  solution  the  same  way  that  normal  GPS  satellite 
data  Is  used.  Since  the  pseudo-satellites  are  ground  stations,  ephemerldes 
are  not  applicable  and  their  position  Is  provided  In  static  earth-fixed 
coordinates.  For  situations  using  algorithms  developed  for  the  range 
application,  this  can  be  taken  into  account.  For  situations  using  User 
Equipment  developed  for  operations  with  satellites,  there  are  potential 
Incompatibilities.  Other  considerations  for  pseudo-satellites  are  that  the 
range  of  Doppler  effects  Is  lower  than  satellites,  the  station  location  is 
accurately  known  In  an  earth  fixed  coordinate  system  (not  GPS),  and  a  refer¬ 
ence  receiver  operated  with  the  pseudo-satellite  would  provide  data  for  re¬ 
lating  these  coordinate  systems. 

Pseudo-satellites  could  be  operated  on  moving  platforms  such  as  aircraft 
and  ships.  Positions  of  the  pseudo-satellites  could  be  determined  from  GPS 
signals  If  sufficient  satellites  or  other  pseudo-satellites  and  an  accept¬ 
able  combination  of  codes  and  frequencies  are  available.  Dynamic  pseudo- 
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satellites  might  be  necessary  for  geometrical  considerations  or  for  supple¬ 
menting  satellite  coverage  In  situations  where  survey  locations  are  not 
available.  Use  of  a  pseudo-satellite  on  or  near  a  target  would  optimize 
relative  oata.  For  targets  in  danger  of  being  destroyed  a  less  expensive 
translator  configuration  could  be  evaluated;  for  example,  GPS  satellite  and 
pseudo-satellite  signals  relayed  through  the  target  translator  to  an  attack 
vehicle  translator  signal  receiver. 


4.6 


Effects  of  Vehicle  Dynamics  on  Receiver  Operation 


In  addition  to  the  "steady  state"  errors  (discussed  In  Section  2)  Inherent 
In  the  GPS  system,  the  dynamics  of  the  vehicle  containing  the  GPS  User  Equip* 
ment  Introduces  additional  problems  and  errors  In  the  operation  of  GPS  navi¬ 
gation.  In  fact,  the  ability  to  even  acquire  and/or  maintain  track  of  the 
GPS  signals  depends  on  the  dynamics  (velocity,  acceleration  and  jerk)  of  the 
vehicle  as  well  as  the  user's  equipment  configuration  and  operational  environ¬ 
ment.  This  section  presents  a  general  discussion  of  the  elements  of  the  GPS 
user  equipment  which  are  affected  by  the  dynamics  and  the  nature  In  which  they 
are  affected. 

4,6.1  Receiver  Configuration 

The  basic  design  of  the  GPS  receiver  (See  Fig.  1.1)  consists  of  from  one  to 
five  channels,  each  In  turn  Incorporating: 

1)  A  Delay  Lock  Loop  (OLL)  for  the  purpose  of  determining  the  pseudo¬ 
range  from  the  user  to  a  specific  GPS  satellite.  The  DLL  achieves 
this  by  measuring  the  time  delay  that  must  be  Introduced  to  an  Inter 
nally  generated  Pseudo-random  Noise  (PN)  code  so  as  to  achieve  cor¬ 
relation  with  the  PN  code  (either  the  C/A  or  P)  on  the  received  sig¬ 
nal  (see  Table  1.1).  This  DLL  will  be  referred  to  as  the  code  loop. 

2)  A  second  loop  Involving  a  Costas  (or  similar)  mechanization  to 
lock  the  received,  Doppler  shifted,  carrier  signal  (L^  or  L£)  with 
an  Internally  generated  oscillator  frequency;  the  Doppler  shift 
(thus  measured)  Is  used  to  determine  the  pseudo-range  rate  of  the 
user  relative  to  an  SV.  This  loop  Is  called  the  carrier  loop. 

3)  A  stable  crystal  oscillator  used  for  generating  the  Internal  signals 

The  time  position  rf  the  inte  nally  generated  PN  code  must  be  continually 
shifted  as  the  pseu_„-  ant  changes  In  time,  and  the  frequency  of  the  inter¬ 
nally  generated  "earner"  must  be  constantly  shifted  as  the  Doppler  shift  of 
the  received  signal  changes  In  time.  Only  when  both  these  loops  are  "locked" 
can  the  navigation  data  contained  In  the  GPS  signal  be  used.  The  process  of 


acquiring  lock  will  be  termed  "acquisition"  and  the  process  of  maintaining 
the  lock  once  It  has  been  achieved  will  be  termed  “tracking".  The  problems 
of  acquiring  and  maintaining  tracu  due  to  vehicle  dynamics  will  be  discussed 
first  followed  by  a  brief  look  at  the  errors  In  the  GPS  measurements  due  to 
the  dynamics. 

4.6.2  Acquiring  and  Maintaining  Track 

Once  the  desired  constellation  of  four  GPS  satellites  has  been  selected  the 
problem  of  acquisition  of  the  signals  arises.  If  the  GPS  receiver  has  four 
or  more  channels  the  acquisition  of  signals  from  four  or  more  satellites  can 
be  carried  on  simultaneously  and  reacquisition  will  only  be  necessary  when 
track  cannot  be  maintained  (either  a  failure  of  lock  or  changing  the  selected 
constellation).  For  receivers  with  less  than  four  channels  the  acquisition 
must  be  carried  out  over  and  over  as  the  satellite  signals  are  processed 
sequentially.  The  problems  that  vehicle  dynamics  cause  In  the  acquisition  and 
tracking  will  be  discussed  first  since  they  are  always  present  and  will  simply 
be  compounded  In  receivers  with  less  than  four  channels. 

Basic  elements  used  In  the  tracking  loops  are;  1)  a  Phase  Detector  (PO)  used 
to  detect  the  phase  error  between  the  received  signal  and  internal  signal  gen¬ 
erated  by  2)  a  Voltage  Controlled  Oscillator  (VCO),  and  3)  a  filter  (normally 
of  second  order  or  higher)  whose  output  of  the  filtered  error  signal  regu¬ 
lates  the  VCO.  The  entire  loop  demonstrates  a  natural  frequency  with  a  char¬ 
acteristic  loop  bandwidth.  Lock  Is  typically  indicated  when  a  small  phase 
error  Is  achieved  (high  correlation  between  the  received  and  internal  signal). 
There  are  several  mechanizations  that  are  used  to  achieve  lock  but  all  depend 
on  obtaining  a  frequency  (or  time  delay)  difference  between  the  received  sig¬ 
nal  and  VCO  generated  carrier  (or  PN  code)  that  is  small  compared  to  the 
bandwidth  of  the  loop.  A  search  for  the  appropriate  frequency  (or  time  delay) 
may  be  made,  but  If  acquisition  time  Is  to  be  relatively  short  then  the  search 
must  be  restricted  to  frequencies  (and  delays)  relatively  close  to  the  right 
one.  In  the  case  of  the  code  loop  this  translates  into  the  need  for  knowing 
an  approximate  pseudo-range  so  th  t  the  delay  induced  in  the  internally  gener¬ 
ated  PN  code  is  close  to  the  pseudo-transit  time.  In  the  carrier  loop  it  is 
necessary  to  know  an  approximate  value  of  the  relative  velocity  between  the 
user  and  the  satellite  so  as  to  calculate  the  Doppler  shift  expected  on  the 


received  signal.  For  a  stationary  or  slowly  moving  user  the  acquisition  prob¬ 
lem  can  be  solved  just  by  using  approximate  knowledge  of  his  position  and  the 
almanac  data  stored  in  his  GPS  receiver.  However,  for  a  highly  dynamic  user 
the  acquisition  problem  can  become  acute  (or  take  an  unacceptable  time)  unless 
some  other  source  of  information,  such  as  from  an  inertial  navigation  system, 
is  available. 

The  wider  the  bandwidth  of  the.  loop  the  faster  the  acquisition  but  unfortun¬ 
ately  the  probability  of  slipping  cycles  increases  also  (due  to  the  worsened 
system  signal  to  noise  ratio  caused  by  the  wider  bandwidth).  Once  both  loops 
have  locked  up,  the  tracking  of  the  code  loop  is  aided  by  delta  pseudo-range 
information  from  the  carrier  loop.  The  narrower  the  loop  bandwidths  the 
better  the  system  signal  to  noise  ratio;  however,  the  ability  to  maintain 
lock  in  the  face  of  dynamics  increases  with  wider  bandwidths.  Thus  the  design 
objective  for  a  particular  application  takes  the  user's  dynamics  into  account 
by  choosing  the  minimum  bandwidths  consistent  with  the  maximum  dynamics  ex¬ 
pected.  Reference  [16],  (page  104)  details  this  trade-off  and  shows  that  the 
carrier  loop  loses  lock  before  the  code  loop. 

Only  approximate  formulas  for  the  time  required  to  acquire  track  have  been 
derived  for  any  practical  GPS  receiver  loops.  Thus,  quantitative  statements 
on  acquiring  and  tracking  are  peculiar  to  each  set  of  User  Equipment.  Gener¬ 
ally  the  acquisition  is  accomplished  by  locking  on  the  C/A  code  first;  then, 
as  the  carrier  loop  locks  up  and  begins  aiding  the  code  loop  a  transfer  is 
made  to  lock  on  the  P  code.  The  reason  for  this  method  is  that  the  C/A  code 
is  only  1023  bits  repeating  every  1  msec  and  thus  much  easier  to  lock  onto  than 
the  one  week  period  P  code.  In  virtually  all  applications  both  the  50  bps 
data  stream  as  well  as  the  hand-over-word  must  be  read  to  acquire  the  P  code. 

In  addition  the  C/A  code  track  is  easier  to  maintain  on  lock  than  the  P  code 
in  cases  of  high  dynamics  (due  partially  to  the  fact  that  the  C/A  code  has 
a  3  dB  higher  signal  to  noise  ratio  than  the  P  code). 

In  addition  to  the  problems  of  acquiring  loop  lock  it  would  appear  that  an 
additional  problem  due  to  vehicle  dynamics  would  be  performing  the  naviga¬ 
tional  computations.  Most  formulations  use  an  iterative  approach  to  the 
navigation  solution  requiring  an  initial  approximation  of  the  position  and 
velocity  of  the  user.  Depending  on  the  degree  of  dynamics  and  the  computa¬ 
tional  cycle  period  the  initial  position  and  velocity  used  in  the  iteration 


could  be  greatly  in  error  and  thus  require  far  more  time  to  converge.  No 
discussion  of  this  problem  was  noted  in  the  literature;  although  it  is  possibly 
reflected  in  the  "time-to-first-flx"  In  the  data  sheets  of  the  specific  GPS 
User  Equipment  packages. 

GPS  receivers  with  fewer  than  four  channels  must  sequentially  acquire  each 
of  four  satellite  signals  thus  acquisition  problems  due  to  vehicle  dynamics 
are  greatly  magnified  unless  the  dynamics  are  so  low  (e.g.,  a  man  on  foot) 
that  acquisition  can  repeatedly  be  made  very  quickly. 

4.6.3  Tracking  Errors 

Generally  speaking  being  locked  in  a  loop  (hence  tracking)  Is  defined  as  the 
condition  of  small  phase  difference  between  the  received  signal  and  the  in¬ 
ternally  generated  signal.  Practical  loops  can  be  demonstrated  (Ref  [55], 

Ch.  4)  to  maintain  lock  as  long  as  the  Doppler  shift  of  the  received  signal  is 
constant  In  time  (the  relative  velocity  of  the  receiver  and  GPS  satellite 
must  be  constant);  obviously  this  is  virtually  never  a  realizable  case.  As 
the  Doppler  shift  changes  in  time  the  loop  will  attempt  to  adjust  so  as  to 
make  the  phase  difference  small  again  but  there  will  be  a  time  lag  in  the 
adjustments  being  made.  These  lags  translate  themselves  into  pseudo-range 
measurement  errors  in  the  code  loop  and  delta  pseudo-range  errors  in  the  carrier 
loop.  Ref.  [56],  (Ch.  18)  gives  a  detailed  mathematical  treatment  of  acquisi¬ 
tion  and  tracking  limits  for  a  delay  lock  loop  as  functions  of  acceleration. 

No  mathematical  derivation  of  errors  in  pseudo-range  and  pseudo-range  rate 
due  to  acceleration  and  jerk  was  found  in  the  available  literature.  One 
could  and  should  be  developed  along  the  lines  of  Ref  [69],  (Section  C.3). 

The  effects  of  acceleration  on  the  frequency  of  the  GPS  crystal  oscillator 
can  cause  tracking  errors  which  appear  as  errors  in  pseudo-ranges.  Unless 
crystal  oscillators  with  very  low  sensitivities  to  acceleration  are  used  (or 
the  effects  compensated  for)  the  errors  for  high  dynamics  aircraft  and  missiles 
with  GPS  receivers  can  become  excessive  (Ref  [70]).  Finally,  what  might  be 
considered  a  vehicle  dynamics  problem  is  the  need  to  keep  the  GPS  user  antenna 
oriented  (or  beam  pointed)  so  as  to  receive  the  GPS  signals  throughout  man¬ 
euvers  for  applications  of  non  omni-directional  antennas. 
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GPS  receivers  with  less  than  four  channels  must  measure  pseudo-range  to  four 
satellites  sequentially.  Since  an  unaided  GPS  receiver  doesn't  measure 
acceleration,  any  relative  acceleration  In  the  period  between  measurements 
is  unaccounted  for  and  creates  errors  in  the  navigated  position  and  velocity. 
(Ref  [2],  pp  12,  13)  states  that  the  position  and  velocity  errors  due  to  this 
effect  are: 

aS  »  1/2  at2 
AV  -  1/2  jt2 

where: 

AS  ■  position  error 
aV  *  velocity  error 

a  ■  vehicle/satell ite  relative  acceleration 
j  ■  time  rate  of  change  of  acceleration 
t  ■  computation  time. 

4.6.4  External  Aids 

In  light  of  the  acquisition  and  tracking  problems  and  errors  for  highly  dynamic 
user  test  vehicles,  most  GPS  applications  for  these  vehicles  are  envisioned  as 
a  complementary  combination  with  other  navigation  systems  (or  other  informa¬ 
tion  sources)  in  order  to  obtain  a  more  accurate  and  reliable  navigation  sol¬ 
ution  than  any  one  system  by  Itself.  This  section  briefly  describes  several 
such  aiding  systems  and  how  they  can  overcome  the  problems  described  previously. 

4. 6. 4.1  Inertial  Navigation  System  (INS)  Integration 

Inertial  Navigation  Systems  (INS)  have  long  been  used  as  stand  alone  guidance 
systems.  Recently  the  technological  gains  made  in  their  accuracies  as  well 
as  miniaturization  and  cost  reduction  has  led  to  extensive  use  in  both  commer¬ 
cial  and  military  aircraft  as  well  as  both  large  and  small  missiles.  The 
errors  In  these  systems  are  well  understood  and  modeled.  When  two  navigation 
systems  such  as  GPS  and  INS  are  integrated  together  many  of  the  weaknesses  of 
both  systems  may  be  overcome. 

A  very  simple  example  of  the  complementary  nature  of  these  two  systems  can 
be  seen  by  recalling  that  one  difficulty  in  rapidly  acquiring  track  in  GPS 
receivers  is  knowing  an  approximate  position  and  velocity  of  the  receiver; 
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the  INS  can  supply  this  data  to  aid  acquisition.  Further,  the  acceleration 
measurements  made  by  the  INS  *  s  accelerometers  can  be  used  to  aid  the  tracking 
loops  In  high  dynamics;  thus,  allowing  the  loop  bandwldths  to  be  narrowed  which 
Improves  the  system  signal  to  noise  ratio  which  In  turn  decreases  the  probabil¬ 
ity  of  skipping  cycles  or  otherwise  losing  track.  The  GPS  can  In  turn  be  used 
to  update  INS  error  estimates  In  an  Integrated  GPS-INS  Kalman  filter  (see 
Section  4.10).  In  addition  the  INS  can  accurately  navigate  by  Itself  during 
reasonable  periods  of  GPS  outage  due  to  the  slowly  varying  nature  of  INS  errors 
even  under  high  dynamics. 

In  addition  to  aiding  acquisition  and  maintenance  of  GPS  track,  the  loop 
tracking  aids  supplied  by  the  INS  can  effectively  eliminate  the  pseudo-range 
measurement  errors  due  to  loop  lags.  In  order  to  optimize  the  tracking  loops' 
performance,  a  process  known  as  adaptive  tracking  may  be  performed  In  which 
the  tracking  loops'  parameters  (e.g.,  bandwidth)  are  changed  dynamically  In 
order  to  adapt  to  changing  environments  (e.g.,  jamming);  the  INS  Inputs  can 
greatly  facilitate  this  adaptive  process.  Finally,  the  attitude  data  derived 
from  the  INS  can  aid  the  control  of  the  beam-pointing  antenna  so  as  to  main¬ 
tain  sight  of  the  GPS  satellites  (Ref  [16],  pp  144-153). 

The  pseudo-range  errors  due  to  high  acceleration  effects  on  the  frequency  of 
the  GPS  crystal  oscillator  can  be  compensated  for,  to  some  degree,  by  the  use 
of  an  accelerometer  (Ref  [71])  thus  an  INS  system  should  be  capable  of  elim¬ 
inating  some  of  these  errors  also. 

Of  course  many  other  measurements  (e.g.,  altitude,  air  speed,  etc.)  may  supple¬ 
ment  the  GPS  when  a  Kalman  filter  Is  used. 

4. 6. 4. 2  A  Priori  Information  Aiding 

In  some  cases  It  Is  practical  to  aid  the  tracking  by  using  a  priori  Informa¬ 
tion  stored  In  the  GPS  processor.  An  example  of  this  Is  a  missile  borne 
receiver  using  a  nominal  acceleration  profile  for  the  pre-determined  trajec¬ 
tory.  References  [46]  and  [47]  describe  a  missile  borne  GPS  receiver  operation 
using  a  priori  aiding,  and  the  "mother-daughter"  concept  discussed  in  Section 
4.9  would  make  use  of  a  similar  aiding. 
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4.7  The  GPS  Receiver  versus  Translator  For 
Mobile  and/or  Stationary  Applications 

4.7.1  Background 

Key  performance  objectives  of  the  GPS  system  which  distinguish  It  from  other 
satellite  and  landbased  navigation  systems  include  the  capability  to  provide 
continuous  and  worldwide  three  dimensional  navigation  coverage  to  users  with 
relatively  high  dynamics.  With  the  advent  of  this  system  It  becomes  practical 
to  measure  the  position  and  velocity  of  high  dynamic  vehicles  such  as  ballis¬ 
tic  missiles,  reentry  bodies  etc.,  with  the  precision  and  accuracy  heretofore 
obtainable  only  through  the  use  of  metric  radar  systems.  However,  one  major 
objection  to  this  system  has  been  the  requirement  for  some  type  of  onboard 
receiver  to  process  and  telemeter  the  position  and  velocity  data  to  a  record¬ 
ing  facility.  This  has  been  accomplished  In  the  past  at  the  WTR  by  the  in¬ 

stallation  of  an  onboard  receiver  (MBRS-MIsslle  Borne  Receiver  Set) (see 
Ref  [46]  and  [47])  Interfaced  to  the  telemetry  system  which  relayed  the  posi¬ 
tion  and  velocity  data  to  ground  recording  sites.  Although  this  system  proved 
to  be  very  accurate.  It  also  proved  to  be  very  expensive  and  subject  to  mal¬ 
function.  The  ETR  currently  uses  a  "bent  pipe"  arrangement  (SATRACK)  de¬ 
signed  to  only  relay  the  received  GPS  signals  to  off-board  data  processing 
sites  where  It  Is  processed  either  In  near  real-time  or  post-flight  to  pro¬ 
vide  the  position  and  velocity  data.  This  system  has  proved  highly  reliable 
and  at  a  significantly  lower  cost  since  all  of  the  processing  equipment  re¬ 
mains  on  the  ground  and  Is  not  lost  at  the  completion  of  the  flight. 

The  WSMC  has  funded  the  Johns  Hopkins  University/Applied  Physics  Laboratory 

( JHU/APL)  to  develop  a  GPS-Sonobuoy  Missile  Impact  Location  System  (GPS-SMILS) 
for  scoring  weapons  systems  reentry  vehicle  impacts  on  a  worldwide  datum. 

The  principal  component  of  this  system  Is  an  expendable  sonobuoy  which  uses  a 
GPS  signal  translator/digitizer  and  relay  system  to  provide  data  to  accurately 
estimate  the  position  of  the  SMILS  array.  In  this  system,  digitized  acoustic 
and  GPS  data  are  multiplexed  and  transmitted  to  an  aircraft  for  recording, 
processing  and  relay  to  the  post  mission  processor  for  further  data  reduction. 
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4. 7. 1.1  Review  of  GPS  C/A  Code  Signal  Transmission  and  Reception 

Before  considering  design  considerations  for  translating,  receiving,  record¬ 
ing  and  processing  GPS  signals,  it  Is  useful  to  review  relevant  signal  char¬ 
acteristics. 

4. 7. 1.2  Signal  Transmission 

The  technique  used  for  GPS  satellite  CA  Code  signal  generation  is  illustrated 
In  Figure  4.7.1.  Note  that  the  transmitted  signal  is  at  1.5  GHz  (nominal) 
and  the  modulation  Is  PCM-PSK  at  a  1.023  MHz  bit  rate  with  50  bps  data  em¬ 
bedded  In  the  PN  Code.  (PSK  modulation  Induces  a  180°  phase  change  of  the 
carrier  each  time  the  PN  Code  changes  state.)  The  RF  main  lobe  signal  spec¬ 
trum  Is  2.048  MHz  wide. 

4. 7. 1.3  Signal  Reception 

A  simplified  block  diagram  of  a  GPS  analog  receiver  channel  Is  shown  in  Fig¬ 
ure  4.7.2.  Received  GPS  signal  levels  may  typically  be  -30  dB  S/N  in  a  2  MHz 
bandwidth.  Each  GPS  satellite  transmits  on  the  same  "carrier"  frequency,  so 
the  received  signals  are  separated  only  by  Doppler  offset.  Therefore,  the  GPS 
signal  spectra  overlap.  In  order  for  a  receiver  channel  to  lock  up  on  one 
signal,  a  code  generator  In  the  receiver  with  the  same  code  as  generated  in 
a  satellite  Is  clocked  into  a  double  balanced  mixer.  For  Initial  lockup, 
the  receivers'  clock  Is  set  somewhat  higher  (or  lower)  than  the  satellite 
clock's  known  rate,  and  each  time  the  code  generator  changes  state  a  180° 
phase  shift  is  Induced  In  each  GPS  signal.  (See  Appendix  A.)  When  the 
receiver's  code  generator  approaches  time  synchronism  with  one  particular 
signal's  code,  the  180°  phase  shifts  which  were  induced  on  that  signal  by 
the  transmitter  are  shifted  another  180°  for  a  total  shift  of  360°,  thus 
restoring  the  carrier.  Code  lock  Is  then  maintained  by  a  tau  dither  loop 
which  Introduces  a  very  small  amount  of  amplitude  modulation  on  the  signal. 

Simultaneously  with  code  lockup,  the  receiver's  down  converter  must  be  swept 
to  place  the  IF  signal  within  the  receiver's  IF  band  pass  filter.  The  IF 
filter  bandwidth  required  to  accommodate  50  bps  data  recovery  is  100  Hz.  3y 
collapsing  the  spectrum  from  2  MHz  to  100  Hz,  the  S/N  is  raised  from  -30  dB 


FIGURE  4.7.1  GPS  SATELLITE  C/A  CODE  SIGNAL  GENERATION  SIMPLIFIED  BLOCK  DIAGRAM 


In  a  2  MHz  bandwidth,  typical,  to  +23  dB  in  a  100  Hz  bandwidth.  After  initial 
lockup,  lock  is  maintained  by  standard  phaselock  techniques. 

4. 7. 1.4  Points  of  Special  Interest 

From  the  proceeding  discussion  of  signal  reception,  several  points  of  specific 
Interest  to  this  note  can  be  made.  First,  consider  a  hypothetical  case  when 
two  of  the  GPS  signals  have  carriers  which  are  offset  in  frequency,  but  with 
Identical  codes  in  time  synchronization.  When  the  receiver  locks  to  one  of 
these  signals,  the  spectrum  of  the  second  signal  is  also  collapsed  as  180° 
shifts  are  Induced  in  both  signals.  This  presents  no  particular  problem  if 
the  carrier  spacing  is  sufficient  to  allow  the  100  Hz  IF  filter  to  pass  one 
signal  and  reject  the  second.  Carrier  spacings  of  a  few  kilohertz  can  easily 
be  acconmodated.  This  feature  allows  GPS  signals  to  be  down  converted  to  a 
very  low  IF  frequency  and  then  upconverted  for  input  to  a  standard  GPS  re¬ 
ceiver. 

A  second  point  of  note  Is  that  the  GPS  receiver  local  oscillator  need  not  be 
derived  from  a  high  stability  reference.  This  Is  because  local  oscillator 
drift  Is  cannon  to  all  receiver  channels  and  therefore  the  restored  carrier 
frequency  differences  are  not  affected.  Similarly,  the  translator's  local 
oscillator  drift  is  cannon  to  all  signals  relayed  through  the  translator. 

A  local  oscillator  stability  of  1  part  in  10®  is  compatible  with  receiver 
lockup  and  ranging  requi rements . 

Finally,  it  should  be  noted  that  interfering  signals  are  spread  during  signal 
reception.  A  CW  signal,  for  example  would  be  PSK  modulated  at  1.023  Mbit  and 
would  have  a  2  MHz  spectrum  as  shown  in  Figure  4.7.1.  This  feature  provides 
excellent  jamnlng  protection.  This  is  especially  important  in  relay  applica¬ 
tions  as  spurs  generated  by  the  relay  transmitter  are  less  troublesome  in 
terms  of  receiver  desensitization. 

4. 7. 1.5  Receiving,  Recording  of  GPS  Signals 

GPS  signals  can  be  received,  recorded  and  played  back  using  conventional  IRIG 
telemetry  pre-detect  equipment  and  techniques.  Referring  to  Figure  4.7.3 


the  L-band  signals  can  be  received  and  up  converted  to  S-band  or  down  con¬ 
verted  to  P-band  or  to  a  receiver's  first  IF.  The  conversion  allows  the 
use  of  conventional  telemetry  receivers  with  IRIG  standard  1.5  MHz  IF's.  . 

The  output  of  these  receivers  is  an  IRIG  standard  pre-detect  signal  centered 
at  900  kHz  which  can  be  patched  to  an  IRIG  standard  telemetry  analaog  recorder 
running  at  120  ips  to  provide  a  2  MHz  frequency  response.  A  2  MHz  reference 
tone  Is  also  recorded  on  the  data  track  in  order  to  reduce  recorder  time 
base  error  during  playback. 

Recording  at  half  speed  can  also  be  accomplished  by  down  converting  the 
signals  to  a  very  low  IF  such  as  30  kHz.  Referring  to  Figure  4.7.3,  both 
In-phase  and  quadrature  signals  are  recorded  on  adjacent  tracks  of  the  same 
record  head.  A  1  MHz  reference  tone  is  utilized. 

4. 7. 1.6  Playback  of  GPS  Signals 

Playback  of  the  signals  for  Input  to  an  analog  receiver  is  accomplished  as 
shown  in  Figure  4.7.4(a).  The  time  base  error  corrector  is  an  off-the-shelf 
IRIG  standard  unit  which  samples  the  signal  in  synchronism  with  the  90° 
crossings  of  the  reference  signal  giving  a  sampling  rate  four  times  higher 
than  the  reference  tone.  The  amplitude  of  the  signal  Is  represented  by  an 
8  bit  word.  These  8  bit  samples  are  then  read  out  at  a  rate  set  by  a  stable 
reference,  and  a  dlgital-to-analog  converter  restores  the  original  signal  with 
sufficient  stability  to  allow  the  receiver  phase  lock  loops  to  acquire  and 
maintain  lock.  There  Is  no  appreciable  Increase  In  noise  due  to  the  8  bit 
samples  taken  at  an  8  MHz  rate,  but  the  total  degradation  of  recording  and 
playback  with  the  suggested  method  has  not  yet  been  measured.  It  Is  expected 
that  there  will  be  a  small  degradation  due  to  residual  recorder  time  base 
error. 

If  a  digital  receiver  is  utilized.  It  Is  not  necessary  to  re-construct  the 
analog  signal.  Rather,  the  sample  rate  can  be  set  to  any  convenient  multi¬ 
digital  processor  as  shown  in  Figure  4.7.4(b).  This  is  convenient  with 
signals  recorded  at  60  ips  as  both  in-phase  and  quadrature  signals  can  be 
simultaneously  processed  at  a  reduced  sample  rate. 
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Design  Considerations 


There  are  two  principal  categories  of  remote  position  determining  systems: 
(1)  the  onboard  receiver  system  is  characterized  by  a  relatively  high  cost 
and  real-time  total  position  processing  capability.  This  system  is  designed 
to  completely  process  the  received  GPS  signals  and  provides  as  an  output  its 
vehicle  position  and  velocity  with  respect  to  time.  (2)  The  second  system 
is  characterized  by  low  cost  and  almost  total  lack  of  GPS  signal  processing 
capability,  it  acts  only  as  signal  relay  equipment.  These  relayed  signals 
are  processed  for  the  vehicle  position  and  velocity  data  at  a  remote  data 
processing  center,  either  near  real-time  or  tape  recorded  for  post  operation 
processing.  A  comparison  of  the  two  catagories  is  shown  in  Table  4.7.1 
(Ref  [16]). 

As  can  be  seen  from  Table  4.7.1,  the  full  GPS  receiver  system  requires  a 
method  df  tracking  loop  aiding,  normally  derived  from  an  onboard  guidance 
or  inertial  system,  facilitate  the  tracking  loops  during  high  dynamic  man- 
uevers  of  the  transporting  vehicle.  In  addition  to  the  simplicity  of  the 
translator/ relay  systems,  aiding  may  not  be  required  in  these  systems  if 
the  data  is  to  be  processed  post-flight.  Real-time  processing  of  relayed 
data  is  complicated  If  timely  acceleration  data  is  not  available  as  an  aid. 

The  onboard  GPS  receiver  may  output  its  position/velocity  solution  to  an 
existing  telemetry  system  or  to  an  independent  transmitter  relay.  In  either 
case,  the  RF  bandwidth  required  to  relay  the  data  is  relatively  narrow,  in 
the  order  of  a  few  kilohertz  and  at  a  relatively  low  power  penalty.  The 
translator/ re lay  system  on  the  other  hand  would  require  an  RF  bandwidth 
equal  to  or  greater  than  the  bandwidth  requirements  of  the  GPS  satellite. 
Thus  the  translator/relay  may  require  In  excess  of  2  MHz  RF  bandwidth  for 
the  C/A  code  relay  and  20  MHz  for  the  P  code  relay.  This  requirement  some¬ 
what  limits  the  relay  data  to  the  C/A  ..ode  and  further  limits  the  number  of 
links  relayed  at  any  one  time  due  to  the  excessive  bandwidth.  A  sequential 
sampling  of  the  multiple  link  would  be  a  possible  solution  but  would  re¬ 
quire  a  scheme  of  coding  unique  to  each  relay. 


TABLE  4.7-1  COMPARISONS  OF  RELAY  AND  RECEIVER  SYSTEMS 


PARAMETER 

RECEIVER  SYSTEM 
(HDUE)* 

RELAY  SYSTEM 

Cost 

high  (*  25K)  (HOUE) 

low  (*  5K) 

Weight 

Obs) 

77 

<5 

Power 

(watts) 

>200 

<50  (Depends  on  relay  trans¬ 
mitter  class  &  power) 

Size 
(cu  ft) 

3.5 

<0.5 

Time  to 
first  fix 
(sec) 

152 

152  (A  rough  position  fix  may 
be  derived  from  less  than  1 
sec  of  cont.  data  through 
computer  analysis  of  the 
relayed  signals) 

Code 

Demodulation 

P,  C/A 

C/A 

Psuedorange 
accuracy  (ic) 

1 .5  meters 

10  meters 

MTBF  (Hrs) 

500 

>2000 

Loop 

tracking 

aid 

Inertial  aids  required 

No  aid  required  for  post  Op 
data  reduction-Aids  required 
for  near  realtime-May  be  der¬ 
ived  from  telemetry,  radar  or 
Doppler  data. 

RF 

bandwidth 

2  kHz  or  less  depending 
on  update  rate 

2  MHz  (C/A  code)  (P  Code  NA) 

Transmitter 

power 

Very  low 

Up  to  20  watts  depending  on 
range  (C/A  code  relay) 

Probability 

of 

acquisition 

0.95  (MBRS) 

>0.95  if  data  reduction  is 
post-fl  ight. 

*HOUE  -  High  Dynamic  User  Equipment 


4.7.3 


Receiver/Translator  Techniques/Designs 


Having  reviewed  the  basic  techniques  of  generating,  receiving,  recording 
and  playback  of  GPS  signals,  it  is  possible  to  readily  understand  the 
operation  of  several  types  of  translators  which  may  be  used  to  relay  GPS 
signals  from  a  moving  vehicle  to  a  receiving  site  and  which  are  compatible 
with  IRIG  standard  receive/record  equipment.  Four  techniques  are  illus¬ 
trated  in  Figure  4.7.5  and  are  discussed  below. 

4.7.3. 1  Onboard  Receiver  Relay 

For  the  onboard  receiver  relay,  a  GPS  receiver  is  placed  on  the  vehicle  and 
the  computed  positons,  velocity,  and  raw  measurement  data  may  be  relayed 
on  a  conventional  telemetry  link.  It  may  be  necessary  to  provide  Doppler 
steering  to  accommodate  vehicle  motion  and  assist  the  onboard  receiver  in 
locking  to  the  satellite  signals.  A  1  kHz  data  rate  is  typically  required 
on  high  dynamic  vehicles.  Tabel  4.7.2  (Ref  [16])  is  a  comparison  of  major 
performance  requirements  for  GPS  systems. 

4. 7. 3. 2  Translator  Systems 

GPS  translator  systems  may  be  grouped  into  two  major  categories,  the  "bent 
pipe"  design  which  performs  no  operations  on  the  GPS  signals  other  than  to 
translate  them  to  a  second  carrier  frequency  then  relay  to  a  processing/ record¬ 
ing  site,  and  those  which  alter  the  signals  by  digitizing,  remodulating,  etc., 
then  relaying.  In  either  system,  it  is  possible  but  usually  not  practical  to 
relay  the  P  code  due  to  the  extremely  wide  bandwidths  involved. 

4. 7.3.2. 1  Bent  Pipe  Linear  Relay 

The  bent  pipe  translator  approach  is  used  on  the  SATRACK  system  associated 
with  the  TRIDENT  program.  The  received  signal  is  converted  to  another  fre¬ 
quency  and  re-radiated.  It  is  necessary  to  radiate  enough  power  to  the  re¬ 
ceiving  site  to  maintain  the  S/N  ratio  achieved  at  the  vehicles  receiver 
input.  It  is  also  necessary  to  retain  amplitude  and  phase  linearity  of  the 
transmitted  signal  which  is  typically  accomplished  with  a  linear  output  amp¬ 
lifier.  A  linear  final  amplifier  introduces  a  5  dB  penalty  in  battery  power 
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FIGURE  4.7.5  GPS  RELAY  TECHNIQUES  SIMPLIFIED  BLOCK  DIAGRAM 


TABLE  4.7.2  GPS  RECEIVER  REQUIREMENT  COMPARISONS 


PARAMETER 

MBRS* 

HOUE** 

MVUE*** 

Code  demodulation 

P, C/A  Code 

P.C/A  Code 

P,C/A  Code 

User  velocity  max. 
(meters/sec) 

7,520 

1 ,100 

25 

User  Acceleration 
max  (meters/sec2) 

100 

30 

6 

User  jerk  max. 
(meters/sec3) 

9 

50 

NA 

Jamming  levels 

J/S  (dB) 

Acquisition 

Track/lock 

Weak  signal 

25 

40 

60 

24  (C/A) 

40  (P) 

47  (P) 

25  (C/A) 

40  (P) 

NA 

Time  to  first  fix 
(seconds) 

120 

152 

240 

Pseudo-range  accuracy 
at  (C/Nq)  (lo) 

1.34m  (25dB-Hz) 
2.4m  (20  dB  Hz) 

1.5m  i30dB-Hz] 

10  m  (30  dB  Hz) 
(position) 

Pseudo-range  rate  (1  a) 
accuracy  (meters/sec) 

.012 

.2 

.3 

Time 

<  1  us 

50  ns 

1  ms 

MTBF  (Hrs) 

Ps  -  0.995 

500 

2,000 

Size  (ft3) 

3.5 

1 .6 

0.75 

Weight  (lbs) 

77 

61 

25 

Power  (watts) 

20B 

200 

45 

MBRS  -  Missile  Borne  Receiver  Set  (Mlnuteman  III) 
HDUE  -  High  Dynamic  User  Equipment 
MVUE  -  Manpack/Vehicular  User  Equipment 


as  compared  to  systems  which  employ  class  C  amplifiers  as  the  final. 

4. 7. 3. 2. 2  Transdigitizer  Relay 

The  C/A  code  transdigitizer  is  being  developed  by  APl/JHU  for  GPS-SMILS.  In 
this  system,  the  received  signals  are  down  converted  to  a  30  kHz,  nominal, 
intermediate  frequency.  Since  the  Doppler  shift  may  be  as  great  as  +  15  kHz, 
and  local  oscillator  uncertainty  as  high  as  1.5  kHz,  the  down  converted  GPS 
"carriers"  may  lie  between  13.5  kHz  and  46.5  kHz.  Following  the  down  conver¬ 
sion,  the  composite  signal  (GPS  signals  embedded  in  noise)  is  hard  limited. 
Unfortunately,  non-linear  operations  such  as  limiting,  drive  the  weaker  GPS 
signals  further  into  the  noise.  The  effect  is  equivalent  to  reducing  the  S/N 
on  the  received  link. 

After  limiting,  the  signal  is  sampled  at  a  2  MHz  rate  providing  a  PCM-NRZ 
wavetrain.  JHU/APL  reports  that  the  combined  limiting  and  sampling  actions 
cause  an  effective  S/N  degradation  on  the  received  link  of  7  dB.  The  sampled 
signal  is  used  to  modulate  a  PSK  transmitter,  the  signal  spectrum  of  the 
transmitted  signal  is  approximately  2  MHz  wide.  QPSK  modulation  of  this 
signal  would  reduce  the  spectrum  to  1  MHz.  The  primary  advantages  of  the 
transdigitizer  approach  are  the  class  C  output  stage  which  utilizes  battery 
power  efficiently  and  the  capability  to  digitally  encrypt  the  relayed  data. 

4. 7. 3. 2. 3  Nonlinear  Relay  Systems 

Nonlinear  relay  systems  are  somewhat  more  complicated  than  the  "Bent  Pipe" 
linear  relay  system  in  that  they  require  some  additional  signal  processing, 
however  power  efficiency  of  these  systems  is  greatly  improved  at  the  expense 
of  some  loss  in  signal -to-noise  ratio  as  compared  with  the  "Bent  Pipe"  design 
but  considerably  better  performance,  than  the  transdigitizer. 

The  received  GPS  signal  is  down  converted  to  the  relay  transmission  frequency, 
limited  and  filtered  to  conserve  transmission  bandwidth,  then  amplified  and 
used  to  directly  drive  a  class  C  final  amplifier.  The  signal-to-noise  ratio 
loss  due  to  the  limiting  action  is  approximately  1  dB.  The  transmission  band¬ 
width  required  is  slightly  greater  than  2  MHz. 
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Bandwidth  of  the  retransmission  link  may  be  conserved  by  down  converting  the 
received  SPS  signals  to  approximately  30  kHz,  limiting  and  filtering  then 
transmitting  at  the  desired  frequency.  The  bandwidth  required  for  this  system 
is  approximately  1.1  MHz  but,  some  loss  (<  3dB)  is  encountered  due  to  the  low 
down  conversion  frequency  and  the  resultant  spectrum  foldover. 

4. 7. 3. 3  Signal -to-Noise  Considerations 

The  translator  serves  the  basic  function  of  receiving  the  GPS  satellite 
signals  and  transmitting  the  signals  to  a  remote  site.  The  transmitted  signal 
"sidebands"  will  "feed  back"  Into  the  receive  channel  as  the  transmit  and  re¬ 
ceive  antennas  and  electronics  are,  of  necessity,  located  in  close  proximity. 
The  relay  trasmltter  emits  "sidebands"  which  are  essentially  distributed 
noise  that  falls  within  the  L-band  receive  spectrum  thus  adding  noise  to  the 
receive  channel.  Therefore,  there  Is  little  advantage  to  using  a  low  noise 
L-band  pre-amp  unless  the  transmitter  induced  noise  at  (L-band)  can  be  reduced 
to  about  the  same  level  as  the  receiver's  pre-amp  noise.  An  L-band  notch  fil¬ 
ter  on  the  transmitter  may  be  required  to  reduce  the  noise  to  an  appropriately 
low  level.  With  regards  to  sinusoidal  spur  interference  from  the  transmitter, 
prudent  design  would  require  that  spurs  be  moved  outside  the  L-band  receive 
bandwidth.  In  any  case,  it  is  important  that  the  L-band  receiver  noise  figure 
be  determined  in  the  operational  configuration  with  the  relay  transmitter  at 
full  power. 

Keeping  the  relay  transmitter  from  interfering  with  the  relay  receiver  favors 
a  design  wherein  the  transmit  frequency  is  significantly  higher  than  the  re¬ 
ceiver  channel.  Small  commercial  units  which  translate  from  L-band  to  S-band 
are  available.  These  units  are  useful  when  the  signal-to-noise  ratio  on  the 
transmitter  link  can  be  maintained  by  high  gain  receiving  antennas  at  the  re¬ 
mote  site.  However,  in  some  applications,  the  receiving  site  may  have  a  low 
gain  antenna,  and  It  may  be  necessary  to  utilize  a  transmit  frequency  which 
is  as  low  as  practical  to  avoid  excessive  space  loss.  In  either  case,  it  is 
usually  advantageous  to  radiate  the  minimum  amount  of  power  to  establish  a 
reliable  link  in  order  to  conserve  the  vehicle's  battery  power. 

When  several  translators  are  used  In  the  same  vicinity,  care  must  be  used  in 
setting  the  relay  link  power  to  ensure  that  sidebands  from  one  translator  do 
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not  Interfere  with  the  mainlobe  signal  of  another  transmitter.  Where  received 
levels  cannot  be  conveniently  matched,  guard  bands  must  be  employed. 

4. 7. 3. 4  Translator  Link  Slgnal-to-Noise  Ratio  Considerations 

The  typical  SPS  satellite  receive  signal  level  is  -160  dBW  from  a  0  dBi 
antenna.  The  received  signal -to-noise  ratio  In  the  translator  RF  amplifier, 
S/NRp,  Is  therefore: 

S/NRF  "  STjJT  or  1n  dB*  "  160  "  101og10  AR 

+  228.6  -  101og1Q  T$  -  62,  dB. 

Where  PR  ■  received  power  from  a  0  dBi  antenna,  -  160  dBW 

AR  ■  L-band  receive  antenna  gain 

K  «  Boltsman's  constant,  -  228.6  dBW/K°  Hz 

T$  «  L-band  receive  noise  temperature  Including  transmitter 
Induced  noise  degradation,  K° 

B  »  GPS  signal  processing  bandwidth,  1.3  MHz. 

The  signal -to-noise  ratio  of  the  transmitted  signal,  S/Nt  is: 

S/Nt.S/NRF  .  -f 

Where  represents  the  signal -to-noise  ratio  degradation  due  to  translator 
"processing". 

For  the  linear  bent  pipe,  T OTog-^ q  — r-  is  0  dB,  while  for  the  limited  bent  pipe, 
It  is  -1  dB  and  for  the  transoigitizer  it  is  -5.5  dB. 

The  relayed  signal  spectra  f r  the  transdigitizer,  bent  pipe  and  limited  bent 
pipe  techniques  can  be  characterized  as  band  limited  noise.  As  the  GPS  signal 


power  Is  very  small  compared  with  the  noise  power,  the  signal -to-noise  ratio 
of  the  relayed  signal,  S/N^  received  at  the  remote  site  is: 

S/Nf  »  S/N^  jr?p“ 
r  n 

Where 

?r  is  the  signal  power  of  the  relayed  signal  at  the  remote  site  receiver 
input  and 

Pn  is  the  noise  power  of  the  remote  site  receiving  system  referenced  to 
the  receiver  input. 

4. 7. 3. 5  Translator  Comparisons 

The  four  translator  techniques  discussed  In  this  note  are  all  compatible  with 
IRIS  standard  receive/record/playback  techniques.  The  on-board  receiver  tech¬ 
nique  requires  the  smallest  data  spectrum,  typically  a  few  kilohertz,  and  is 
therefore  useful  in  an  arena  where  many  translators  will  be  used  simultaneously. 
It  Is  especially  advantageous  when  relay  from  the  vehicle  to  another  vehicle 
(such  as  a  satellite)  via  a  standard  narrowband  conmuni cations  channel  is 
needed.  However,  the  complexity  and  cost  as  well  as  the  requirement  for  on¬ 
board  Doppler  steering  may  well  preclude  widespread  use  of  the  on-board  re¬ 
ceiver  In  many  test  range  applications.  In  addition,  the  on-board  receiver 
battery  and  transmitter  power  penalties  are  typically  10  dB  because  of  the  re¬ 
quirement  to  maintain  a  10  dB  signal -to-noise  ratio  in  a  bandwidth  equal  to  the 
data  bit  rate. 

With  regards  to  the  linear  bent  pipe,  limited  bent  pipe  and  transdigitizer 
techniques,  all  three  should  be  similar  in  terms  of  complexity,  cost  and  spec¬ 
tral  requirements.  However,  the  transdigitizer  suffers  from  severe  transmitter 
power  penalties  as  shown  in  Figure  4.7.6*.  The  limited  bent  pipe  offers  the 
potential  for  a  3  dB  battery  power  savings  if  a  3  dB  signal -to-noise  ratio 


When  both  telemetry  data  and  GPS  signals  are  to  be  simultaneously  radiated  from 
the  same  vehicle,  however,  the  transdigitizer  becomes  a  viable  option. 

See  4. 7. 3. 6 


FIGURE  4.7.4  TRANSMITTER  POWER  AMD  TRANSMITTER  FINAL  BATTERY  PENALTIES  FOR  EQUIVALENT 

SIGNAL -TO -NOISE  PERFORMANCE  AT  THE  RECEIVE  SITE 


degradation  Is  tolerable.  In  general,  however,  when  the  vehicle's  L-band 
receive  antenna  has  low  gain  (perhpas  -10  dBi  or  lower),  it  may  be  necessary 
to  utilize  the  linear  bent  pipe  in  order  to  minimize  the  up-link  degradation 
without  having  to  overpower  the  link.  A  comp-'ison  of  the  important  character¬ 
istics  of  each  relay  system  is  outlined  in  Table  4.7.3. 

4. 7. 3. 6  Telemetry  Data  Relay 

In  many  range  support  applications,  the  vehicle  may  be  designed  to  carry  a 
telemetry  transmitter  and  antenna,  and  it  may  be  planned  to  add  a  GPS  trans¬ 
lator  to  the  vehicle.  However,  such  an  addition  may  be  impractical  because 
of  size,  weight  or  spectrum  constraints.  In  such  cases,  consideration  should 
be  given  to  developing  a  package  of  the  same  size  and  weight  as  the  telemetry 
package  alone  but  which  provides  both  telemetry  transmission  and  GPS  trans¬ 
lation.  Many  schemes  can  be  conceived  to  transmit  both  GPS  and  telemetry  data 
on  a  common  link  using  a  single  transmitter  final.  When  the  telemetry  data 
Is  PCM,  a  convenient  approach  is  a  QPSK  transmission  scheme  wherein  the  trans- 
digitizer  stream  Is  clocked  at  2  MHz  and  the  PCM  data  rate  is  set  to  1  Mbit 
or  a  convenient  sub-multiple  of  1  Mbit.  With  QPSK,  one-half  the  transmitted 
power  Is  in  the  telemetry  stream,  and  one-half  is  in  the  GPS  stream,  and  a 
10  dB  S/N  is  required  by  the  QPSK  receiver  demodulator  (in  a  bandwidth  equal 
to  the  telemetry  data  rate).  All  other  co-channel  transmission  techniques 
considered  to  date  require  the  same  or  greater  data  spectral  occupancy,  as 
the  QPSK  technique,  and  several  other  techniques  can  potentially  be  used  to 
reduce  the  transmitter  power  required  by  nearly  3  dB.  However,  the  possibil¬ 
ity  of  providing  secure  GPS  and  secure  telemetry  transmission  makes  the  QPSK 
approach  an  attractive  option.  * 

4.7.4  Multi  Vehicle  Interference 

In  most  cases  where  there  are  only  a  few  vehicles  (*  10)  whose  positions  must 
be  determined,  the  wide  transmission  bandwidth  required  by  the  translator  sys¬ 
tems  would  pose  few  problems;  however,  in  situations  such  as  military  field 
maneuvers,  where  possibly  hundreds  of  vehicles  would  require  positioning,  the 
problem  of  bandwidth  becomes  more  complex  and  alternate  arrangements  may  be 

★ 

When  links  are  secured,  a  10  dB  S/N  ratio  must  be  maintained  in  a  2  MHz 
bandwidth. 
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TABLE  4.7.3  COMPARISON  OF  RELAY  SYSTEM  CHARACTERISTICS 


PARAMETER 

■■ 

Transdigitizer 

Non  Linear 

Relay 

RF  Bandwidth 

2.1  MHz  (C/A) 

1.1  (C/A) 

2.1  MHZ  (C/A) 

required  (code) 

10.1  MHZ  (P) 

5.1  (P) 

(QPSK  Modulation) 

10.1  MHz  (?) 

Type  of  relay 
Transmitter 

Linear 

Class  C 

Class  C 

Remodulation 

Loss  (dB) 

0 

7 

1  to  3  dB 

Battery  power 
drain 

High 

(Oue  to  class  A 
final ) 

Low 

Low 

Compatible  with 
existing  range 
receive/record 
equipment 

No 

Yes 

P  Code  -  No 
(bandwidth  limit) 

C/A  Code  -  yes 

P  Code  -  No 
(bandwidth  limit) 

Compatible  with 
standard  GPS 
Receiver 

Yes 

Yes 

Yes 

Susceptibility  to 
jamming  on  relay 
link 

Low 

(Same  as  down 

1  Ink) 

High 

Low 

(Same  as  downlink) 

Susceptlbil ity 
to  interception 

High 

Low  IF 

Encrypted 

High 

required.  In  low  vehicle  dynamic  situations,  translator  center  frequencies 
may  be  spaced  as  close  as  1  kHz.  This  close  spacing  causes  the  transmitted 
spectrums  to  overlay,  however,  due  to  the  narrow  bandwidths  (<  100  1  )  employ¬ 
ed  in  the  processing  center  GPS  receivers,  the  overlapping  spectrums  will  not 
interfere.  Since  onboard  receiver  systems  require  only  very  narrow  bandwidths 
to  relay  their  information,  spectral  occupancy  is  not  a  problem. 

4.7.5  Same  Examples  of  Potential  Translator  Applications 

Figure  4.7.7  illustrates  several  potential  translator  applications.  As  shown, 
"Worldwide"  aircraft  testing  can  be  accomplished  using  an  on-board  receiver 
relaying  to  a  communications  satellite. 

Cruise  missile  type  vehicles  can  be  supported  worldwide  using  a  tracking 
aircraft  as  an  intermediate  relay  to  a  communications  satellite. 

Ballistic  vehicles  can  be  supported  by  relaying  to  telemetry  tracking  sites. 

Broad  ocean  scoring  can  be  accomplished  by  relaying  to  the  telemetry  support 
aircraft. 

Note  that  each  application  tends  to  favor  the  use  of  a  particular  kind  or 
kinds  of  GPS  signal  relay  because  of  the  mission  peculiar  vehicle  features 
and/or  mission  requirements. 
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EXAMPLES  OF  POTENTIAL  TRANSLATOR  APPLICATIONS 


4.8 


Access  to  Raw  Oata  Considerations 


This  section  describes  the  type  of  data  that  ma^  be  available  from  a  GPS  re¬ 
ceiver  and  the  potential  effect  on  processing.  As  an  example  the  Missile 
Borne  Receiver  Set  (MBRS)  receiver  flown  on  Mlnuteman  III  output  pseudo-range, 
range  rate  (delta  range),  computed  position,  velocity  and  acceleration,  ephem- 
erls  information  and  receiver  status  data.  These  data  allowed  for  real-time 
evaluation  using  the  on  board  GPS  computed  position  and  velocity  output,  and 
flexible  post-flight  processing  using  pseudo-range,  range  rate  and  ephemeris 
data.  This  receiver  did  not  exercise  the  dual  frequency  ionospheric  refrac¬ 
tion  capability  forcing  data  users  to  apply  an  a  priori  refraction  model  or 
externally  supplied  correction  (from  Aerospace  Corp.)  or  else  neglect  this  term 


In  contrast  the  Missile  Accuracy  Evaluation  (MAE)  receiver,  that  was  to  have 
been  flown  on  the  MX  missile,  was  to  have  had  a  wide  variety  of  raw  data  and 
receiver  performance  data  as  well  as  dual  frequency  output.  The  potential 
user  would  have  been  required  to  obtain  GPS  ephemeris  data  from  a  ground  re¬ 
ceiver  (for  real-time  applications)  or  wait  for  corrected  ephemeris  data  to 
be  supplied  for  post-flight  use.  Use  of  the  MAE  data  in  real-time  would  have 
required  far  more  ground  based  equipment/software  than  MBRS  but,  for  post-test 
applications,  the  MAE  data  would  probably  have  yielded  a  more  accurate  traj¬ 
ectory  reconstruction.  Subsequent  subsections  will  detail  these  issues  further 


4.8.1 


Raw  (R,  R)  Versus  Processed  (P,  V)  Data 


"Raw11  data  (pseudo-range  and  pseudo-range  rate)  is  somewhat  of  a  misnomer  as 
these  data  are  calculated  by  a  GPS  receiver  and  represent  the  basic  data  for 
subsequent  processing.  The  sample  rates  of  these  data  vary  but  .an  be  as  high 
as  50  sps.  In  the  case  of  an  onboard  receiver  (such  as  MBRS)  it  is  required 
that  these  data  be  multiplexed  into  a  down  link  (usually  S-band)  for  acquisi¬ 
tion  and  subsequent  proessing.  Care  must  be  taken  to  ensure  that  all  measure¬ 
ment  samples  are  recoverable.  Under  some  situations  where  the  down  link  con¬ 
tains  other  data,  such  as  INS  parameters,  it  is  possible  that  the  down  link 
frame  rate  is  asynchronous  with  the  raw  data  rate  resulting  in  the  potential 
loss  of  some  data. 


.8- 


The  availability  of  raw  data  allows  the  user  to  apply  processing  algorithms 
tailored  to  his  requirement  and  computer  capability.  These  requirements  may 
be  quite  different  as  to  real-time  versus  post-test. 

Most  military  GPS  receivers  output  position  and  velocity  data,  in  various  co¬ 
ordinate  systems,  for  either  visual  CRT  display  or  input  to  a  navigation  sys¬ 
tem.  Extensive  use  of  microprocessors  are  made  for  Internal  data  processing 
in  these  receivers  with  the  effect  of  limiting  the  parameters  output  to  a 
minimum.  This  category  of  receivers  cover  most  "off  the  shelf"  equipment  with 
the  exceptions  being  the  R&D  receivers  (such  as  MBRS  and  MAE)  and  some  high 
dynamic  receivers.  It  should  be  stressed  that  for  most  Range  applications  a 
receiver  providing  flexible  parameter  output  (and  possibly  aiding  inputs  - 
due  to  the  many  high  dynamic  applications)  will  be  required.  In  low  and  some 
medium  dynamic  field  environments  where  real-time  positioning  is  desired,  off 
the  shelf  receivers  (such  as  the  manpack)  are  suitable.  It  should  also  be 
noted  that  many  test  range  application  receivers  operate  in  a  controlled  en¬ 
vironment  with  technical  support  and  maintenance  thus  reducing  the  need  for 
stringent  reliability  requirements. 

4.8.2  Code  and  Carrier  Loop  Error  Signals 

In  order  to  acquire  a  GPS  signal  it  is  necessary  to  obtain  code  correlation 
with  a  receiver  generated  replica  code  (sometimes  aided  by  an  anti- jamming 
T  code).  Once  code  lock  has  been  achieved,  a  carrier  loop  is  used  to  collapse 
the  spectrum  for  the  50  bps  data  acquisition  block  (1500  bits)  and  to  further 
recover  the  carrier  to  aid  the  code  correlation  function.  The  bandwidth  of  the 
carrier  loop  (and  its  design)  becomes  a  trade  off  between  higher  carrier  SNR 
and  dynamic  tracking  performance  (see  Section  4.6).  It  should  be  noted  that 
the  carrier  loop  will  only  lock  as  long  as  code  correlation  is  maintained  so 
the  net  effect  Is  to  have  a  carrier  aided  code  tracking  loop  and  not  continu¬ 
ous  carrier  tracking  (see  Section  4.6). 

The  common  code  loop  contains  early,  late  and  prompt  (or  punctual)  outputs 
that  determine  the  status  of  the  correlation.  The  input  to  the  code  loop 
is  either  a  VCO  or  0C0  that  drives  the  replica  PN  generator.  It  may  be  desir¬ 
able  to  have  status  indications  of  the  mode  of  this  loop  such  as  C/A  or  P  code 
track;  early,  late,  and  prompt  status  and  the  PN  generator  clocking  rate. 


The  narrow  band  carrier  loop  may  consist  of  a  phase  lock  loop  and  a  frequency 
lock  loop.  The  carrier  loop  control  signals  are  provided  from  a  Costas  error 
generator.  The  error  signals  that  are  provided  from  this  narrow  band  loop 
are  similar  to  that  of  the  wideband  code  loop  but  may  also  indicate  carrier 
cycle  slip.  Other  types  of  narrow  band  carrier  lock  loops  are  the  non  coher¬ 
ent  delay  lock  loop  and  the  t  dither  loop.  The  choice  of  design  is  generally 
dependent  on  the  expected  dynamics  and  sequential /parallel  channel  design. 

All  of  these  signals,  for  both  loops,  can  be  used  for  analysis  purposes  and 
for  possible  aid  in  post-test  data  reduction. 


4.9  Standard  User  Hardware  Characteristics 

Although  there  are  numerous  different  GPS  user  hardware  units  already  avail¬ 
able  and  under  development,  the  consideration  here  will  be  mainly  restricted 
to  those  currently  in  the  full  scale  engineering  development  phase  for  the 
GPS  JPO.  Magnavox  and  Collins  were  contracted  to  perform  the  Phase  II  Full 
Scale  Development.  The  previous  four  years  were  devoted  to  Phase  I,  Concept 
Validation  development  and  testing  and  involved  Texas  Instruments  as  well. 
This  section  will  describe  the  overall  design  concepts  of  the  user  hardware 
being  developed  and  tested  in  Phase  II  as  well  as  a  few  more  specialized  GPS 
units  which  may  be  useful  for  Range  applications. 

The  design  philosophy  has  been  dictated  by  the  concept  of  using  GPS  equip¬ 
ment  on  hundreds  of  different  vehicles  whose  dynamics  vary  from  a  stationary 
user  to  high  velocity  and  acceleration  aircraft  and  missiles.  In  order  to 
produce  this  range  of  units  at  a  reasonable  cost  the  design  goal  was  to  pro¬ 
duce  the  minimal  number  of  unique  hardware  (as  well  as  software)  modules 
which,  through  various  combinations,  could  meet  all  the  performance  and  host 
vehicle  requirements.  The  "System  Segment  Specification  for  the  User  System 
Segment,  NAVSTAR  Global  Positioning  system.  Phase  II,"  (United  States  Air 
Force,  SS-US-200,  31  January  1979)  should  be  consulted  for  details.  Most  of 
what  follows  borrows  heavily  from  Ref  [37]. 

4.9.1  Phase  II  Hardware  Design  Overview 

The  modularization  can  be  broken  into  Line  Replaceable  Units  (LRU)  which  are 
plug-in  compatible  units  which  can  be  used  in  various  combinations  to  meet 
several  requirements  -  Many  more  variations  can  be  formed  (and  cost  savings 
achieved)  at  the  level  of  the  Shop  Replaceable  Unit  (SRU);  for  instance,  the 
circuit  boards  used  In  a  single  channel  receiver  can  be  replicated  to  form 
a  multichannel  receiver.  A  description  of  the  basic  LRU's  and  an  indication 
of  the  type  of  applications  which  would  dictate  their  use  follows. 

4. 9. 1.1  Antenna/Antenna  Electronics 

A  user  obviously  must  have  an  antenna  capable  of  receiving  the  signals  trans¬ 
mitted  from  the  GPS  satellite.  The  simplest  of  these  are  omni-direitional 


antennas  suitable  for  a  user  unlikely  to  encounter  jamming  and  whose  vehicle 
dynamics  combined  with  antenna  placement  insures  no  shading  of  the  antenna; 
the  design  may  vary  but  all  are  referred  to  as  Fixed  Recept  n  Pattern  Antennas 
(FRPA).  For  applications  requiring  more  complexity  there  are  Controlled  Recep¬ 
tion  Pattern  Antennas  (CRPA)  which  have  control  electronic  modules  associated 
with  them.  The  most  prevalent  type  of  CRPA  senses  jamming  energy  arriving 
from  certain  directions  and  adjusts  the  antenna  pattern  so  as  to  create  nulls 
in  those  directions. 

4. 9. 1.2  Receiver/Processors 

The  GPS  receiver/processor  consists  of  hardware  and  software  necessary  to  ac¬ 
quire  and  track  a  GPS  signal,  demodulate  the  navigation  message,  and  compute 
the  navigated  position  and  velocity  of  the  user.  The  receiver  section  might 
be  considered  the  signal  processing  sections  comprising  the  tracking  loops, 
etc.;  while  the  processor  generates  the  appropriate  PN  code,  performs  control 
functions,  has  a  memory  containing  the  almanac,  and  performs  navigation  com¬ 
putations,  etc.. 

a.  The  simplest  receiver  has  only  a  single  channel;  hence,  it  must 
sequentially  track  four  satellites  to  obtain  the  data  necessary  for 
navigation.  Phase  I  tests  revealed  that  such  a  receiver  can  meet  the 
combined  jamming,  dynamics,  accuracy,  and  time-to-first  fix  speci¬ 
fications  for  troops  and  land  vehicles. 

b.  For  medium  dynamics  vehicles  (e.g.,  helicopter)  a  two  channel  re¬ 
ceiver  is  adequate.  This  cycle  functions  much  like  the  one  channel 
but  reduces  measurement  cycle  time. 

c.  High  dynamics  aircraft  and  missiles  (as  well  as  submarines  which 
need  to  have  a  very  short  time  to  first  fix)  requires  a  five  channel 
receiver  and  in  some  applications  aiding  will  be  necessary.  The 
five  channel  receiver  can  be  tracking  five  satellites  simultaneously 
so  that  changing  constellations  (of  four)  need  not  interrupt  the 
navigation  solution. 


4. 9. 1.3  Flexible  Modular  Interface  (FMI) 


These  units,  where  required,  contain  the  hardware  and  software  necessary  to 
integrate  the  GPS  set  into  the  specific  host  vehicle.  Since  these  units  are 
the  most  applications  specific,  the  flexibility  must  be  at  the  SRU  level. 

The  different  developers  seem  to  have  very  flexible  concepts  of  what  the  FMI 
comprises;  for  instance,  some  contain  CPU's  themselves  with  resident  programs 
that  the  processors  may  pull  in.  This  unit  can  serve  as  the  interface  with 
"aiding"  measurements  from  inertial  navigation  systems,  altimeters,  etc.. 

As  indicated  in  Ref  [72]  the  number  and  type  of  output  channels  has  been  left 
up  largely  to  the  discretion  of  the  manufacturer;  this  may  be  the  way  by  which 
Test  Ranges  may  gain  access  to  the  "raw"  GPS  data  which  would  greatly  aid 
post-test  processing. 

4. 9. 1.4  Control  and  Display 

This  unit,  where  needed  or  desirable,  serves  to  display  the  GPS  data,  such 
as  position  and  velocity,  and  to  allow  manual  input  of  data.  This  unit  is 
not  needed  in  applications  where  navigation  control  panels  already  exist. 

4.9.2  Special  Applications  GPS  Hardware 

Although  the  Phase  II  designs  will  cover  most  applications,  there  are  specific 
uses  where  cost  savings  can  be  realized  using  simpler  hardware. 

4. 9. 2.1  Frequency  Translator 

For  vehicles  which  are  expendable  a  simple,  inexpensive  frequency  translator 
is  being  developed.  The  translator,  located  on  the  vehicle  whose  position 
and  velocity  is  to  be  determined,  simply  receives  all  the  GPS  signals  it  can 
and  retransmits  them  on  S-band.  A  surface  (or  an  auxilliary  vehicle)  based 
GPS  receiver  (modified  to  receive  S-band)  decodes  the  signal  and  computes  the 
location  of  the  vehicle.  Thus,  a  four  channel  ground  based  receiver/processor 
can  determine  the  position  and  velocity  of  the  vehicle  with  almost  the  same 
accuracy  as  a  four  channel  receiver/processor  on  the  vehicle  while  saving 
the  electronics  that  would  necessarily  be  destroyed  (see  Section  4.7). 


4. 9. 2. 2  Tactical  Mid-Course  Guidance 


Another,  slightly  more  complicated,  unit  designed  for  expendable  use  (i.e., 
tactical  missiles)  is  based  on  a  so  called  "mother-daughter"  concept.  The 
missile  (the  "daughter")  contains  a  simplified  receiver  which  need  only  ex¬ 
tract  the  pseudo-range  data  which  it  uses  in  conjunction  with  its  inertial 
navigation  system.  The  launching  vehicle  ("mother")  contains  a  complete 
GPS  receiver  and  preselects  the  constellation,  initializes,  aligns,  and 
calibrates  the  daughter  receiver  prior  to  launch.  This  technique  is  for  re¬ 
latively  short  flight  tactical  missiles. 

4.9.3  Configurations  for  Test  and  Training 

Following  Ref  [2]  it  is  assumed  that,  with  some  modifications,  the  hardware 
just  described  could  be  used  for  test  and  training  ranges  also.  Primarily, 

I  some  provision  to  transmit  the  GPS  data  to  a  ground  site  must  be  made;  it 

would  also  seem  desirable  to  transmit  some  of  the  raw  data  (see  Section  4.8). 

It  should  be  stressed  that  this  author  has  no  knowledge  of  whether  the  designs, 
say  on  the  SRU  level  and  in  the  software,  have  any  requirement  to  make  this 
transmission  option  available. 

For  purposes  of  classifying  the  types  of  vehicles  to  be  equipped  with  GPS  units. 
Ref  [2]  defines  the  following  four  categories. 

A.  High  velocity  with  high  acceleration  -  This  is  the  class  of  vehicles 
with  velocities  over  200  kts  and  accelerations  over  5  g;  hence,  it 
Includes  many  missiles  and  the  high  performance  aircraft. 

B.  High  velocity  with  medium  acceleration  -  These  vehicles  have  veloc¬ 
ities  over  200  kts  but  accelerations  below  5  g;  examples  are 
tactical  aircraft  and  helicopters. 

C.  Low  velocity  with  medium  acceleration  -  Vehicles  with  velocities 
below  200  kts  and  acceleration  less  than  5  g  such  as  ships  and 
ground  vehicles  (perhaps  representing  moving  targets  for  testing). 


0.  Low  velocity  with  low  acceleration  -  This  category  is  primarily 
for  ground  units  which  are  basically  stationary  but  move  from  time 
to  time  during  an  exercise. 

Since  these  categories  were  defined  in  particular  for  Navy  Test  Range  pur¬ 
poses,  it  is  quite  possible  they  will  have  to  be  modified  and  extended  for 
a  Triservices  application.  (Translators  and  pseudo-range  receivers  are 
examples.) 

Again  following  Ref  [2],  five  basic  SPS  hardware  configurations  (with  aiding 
where  necessary)  will  now  be  described  to  cover  the  four  test  and  training 
categories.  For  most  of  them  a  down  link  transmitter  and  antenna  for  trans¬ 
mitting  raw  and  processed  data  to  a  ground  (or  airborne)  receiving  site 
would  be  necessary  for  test  purposes. 

Type  I  -  High  performance/airborne  -  This  configuration  consists  of  an  inte¬ 
grated  inertial  guidance  system  and  a  two  frequency  (L^  and  L£),  five- 
channel  GPS  receiver.  This  type  Is  for  use  with  category  A  user  aircraft. 

Type  II  -  High  performance/ expendable  -  For  category  A  vehicles  which  are  not 
recoverable  one  of  the  special  applications  units  may  be  used.  One  example 
would  be  a  ballistic  missile  equipped  with  an  INS  and  a  GPS  translator. 

The  ground  station  would  receive  the  translated  GPS  data  as  well  as  the  INS 
acceleration  data  (so  as  to  maintain  lock  under  high  dynamics  maneuvers). 
Another  example  would  be  the  tactical  missile  case  using  the  "mother-daughter 
configuration.  The  data  that  would  be  needed  to  be  transmitted  to  a  ground 
station  is  not  clear  though. 

Type  III  -  High  speed/airborne  -  This  configuration  is  intended  for  use  with 
category  B  vehicles  and  with  category  C  vehicles  that  require  high  accuracies 
It  consists  of  a  two  frequency,  five  channel  GPS  unaided  unit  with  at  least 
a  10  Hz  update  rate. 

Type  IV  -  Low  speed /airborne  -  This  Is  intended  for  category  C  vehicles  that 
can  tolerate  a  relatively  low  accuracy  and  with  category  0  vehicles.  The 
unit  is  a  two  frequency,  two  channel  unaided  GPS  unit. 


Type  V  -  Manpack  -  This  configuration  of  a  two  frequency,  single  channel  unit 
will  give  accuracies  sufficient  only  for  category  D  users. 

In  addition  there  should  be  a  sixth  category  to  cover  the  possible  use  of 
translators  for  applications  such  as  SPS-SMILS. 

Table  4.9.1  (reproduced  from  Ref  [2])  illustrates  the  accuracies  expected 
from  the  five  configurations. 

4.9.4  Specifications  of  an  Existing  Hardware  Configuration 

For  concreteness,  an  example  of  an  existing  configuration  and  its  specifica¬ 
tions  will  be  described.  This  is  drawn  from  Collins  (Ref  [45]);  similar 
data  from  Magnavox  was  not  at  hand. 

Type  I  configuration/attack  aircraft  -  This  GPS  unit  Is  capable  of  stand  alone 
or  aided  navigation  solutions  for  position,  velocity,  and  time  for  a  high 
performance  aircraft.  The  antenna  system  consists  of  a  bottom  mounted  FRPA 
and  a  top  mounted  CRPA  (with  six  auxiliary  elenents  to  provide  spatially 
adaptive  nulling  for  simultaneous  jamming  sources);  both  are  capable  of  re¬ 
ceiving  both  the  L1  and  L2  frequencies.  An  antenna  electronics  module  re¬ 
ceives  signals  from  both  antennas  simultaneously  and  does  the  null  steering 
process  for  the  CRPA.  The  receiver/processor  module  has  five  channels  and 
consists  of  a  central  processing  unit  which,  along  with  distributed  micro¬ 
processors,  performs  the  signal  processing  and  navigation  tasks.  There  is 
a  twelve  element  Kalman  filter  for  the  integrated  processing  of  GPS  measure¬ 
ments  and  host  vehicle  aiding  sensors  (inertial  navigation,  Doppler  radar, 
altimeter,  leading/roll /pitch  sensor,  and  true  airspeed  sensor).  Fix  up¬ 
dates  are  once  a  second  and  can  be  propagated  up  to  20  times  a  second.  The 
FMI  subassembly  Is  connected  to  the  receiver/processor  unit  and  serves  as  the 
interface  with  the  aircraft's  other  related  systems.  It  contains  a  second 
CPU  for  weapons  delivery  computations  using  the  navigation  solution.  There 
Is  an  optional  data  loader  subsystem  for  bulk  loading  host  vehicle  and  mis¬ 
sion  unique  data  and  an  optional  control /display  unit. 

The  GPS  unit's  accuracy  In  a  stand  alone  mode  is  quoted  as:  Position  -  10  m; 
Velocity  -  .1  m/s;  and  time  -  50  ns.  These  are  for  dynamics  up  to: 
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Velocity 
Acceleration 
Jerk 
Yaw 
Pitch 
Roll 

When  significant  jamming  Is  present  these  accuracies  can  be  maintained  only 
with  aiding. 

This  unit  requires  4.5  minutes  from  power  on  to  first  accurate  solution. 

Its  Tlme-To-First-FIx  (TTFF)  after  initialization  is  1.5  minutes. 

Examples  of  configurations  suitable  for  the  lesser  dynamic  applications  can 
be  found  In  References  [39]  through  [44], 

4.9.5  GPS  Receiver  and  Translator  Sources 

In  addition  to  the  JPO  contracted  companies  (Magnavox  and  Collins),  the 
following  are  possible  sources  of  User  Equipment  for  range  applications. 

a.  Texas  Instruments  -  In  addition  to  TI's  work  developing  a  sequential, 
time  multiplexed  single  hardware  channel  receiver  for  commercial  use; 
they  are  developing  a  “mother -daughter"  configuration  for  the 
Tomahawk  missile. 

b.  Stanford  Telecommunications  Inc.  -  STI  has  a  C/A-L^  receiver  design 
under  development  for  commercial  use. 

c.  Interstate  Electronics  Corp.  -  IEC  has  designs  for  an  "all  digital" 
full  C/A -P,  Lj,  l2  receiver  as  well  as  a  C/A  -  L1  receiver. 

d.  RCA  -  Designs  for  a  SATRACK  translator. 

e.  Cubic  Corp.  -  Designs  for  a  SATRACK  translator. 

f.  Applied  Physics  Lab  -  Designing  a  translator  for  GPS-SMILS. 


1200  m/s 
118  m/s/s 
180  m/s/s/s 

+  1  rad/s,  +  3  rad/s/s 
+  1  rad/s,  +  6  rad/s/s 
+  5.5  rad/s,  +  17.5  rad/s/s 
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Filtering  and  Smoothing 


The  problem  of  processing  a  set  of  measuremer.ts  to  produce  "optimal"  estim¬ 
ates  of  selected  random  variables  has  been  of  some  interest  since  at  least 
the  early  nineteenth  century  (Ref  [65]).  In  recent  years  the  widely  success¬ 
ful  application  of  linear  estimation  algorithms  has  been  most  notable  (Ref  [66], 
[67],  [68]). 

Linear  estimation  algorithms  have  generally  been  either  of  the  "batch"  or 
the  "recursive"  type.  In  the  former,  all  measurements  are  operated  upon  at 
once,  while  in  the  latter,  measurements  are  processed  a  few  at  a  time,  gener¬ 
ally  in  time  sequence.  The  Bayesian  Weighted  Least  Squares  (ULS)  method  is 
the  most  general  batch  processing  algorithm;  and  the  Kalman  Filter/Smoother 
(KFS)  method  is  the  most  general  recursive  algorithm. 

While  there  exist  many  similarities  and  much  commonality  of  application  be¬ 
tween  the  WLS  and  KFS  methods,  the  KFS  method  actually  is  applicable  to  a 
larger  class  of  problems  than  the  WLS  method.  For  this  reason,  the  KFS  method 
is  considered  exclusively  in  the  sequel. 

The  problem  of  processing  measurements  to  obtain  optimal  estimates  of  selected 
variables  of  Interest  Is  a  fundamental  one.  In  the  general  case,  the  problem 
is  formulated  by  expressing  the  variable  of  Interest  together  with  augmented 
auxiliary  variables,  collectively  called  state  variables,  as  dependent  vari¬ 
ables  in  a  set  of  stochastic  differential  equations  in  which  time  is  the  in¬ 
dependent  variable.  The  measurements  are  then  expressed  as  stochastic  func¬ 
tions  of  the  state  variables  and  time. 

A  brief  mathematical  treatment  In  Appendix  A  demonstrates  the  application  of 
the  KFS  method  to  the  estimation  problem  formulated  above. 

In  order  to  apply  the  KFS  algorithm,  it  is  necessary  to  formulate  a  state 
variable  representation  for  the  quantities  of  interest.  In  GPS  applications 
these  variables  generally  Include  position  and  velocity.  Moreover,  measure¬ 
ment  errors  either  in  pseudo-range  or  from  auxiliary  sources  such  as  iner¬ 
tial  sensors,  which  are  not  sequentially  uncorrelated  (i.e.,  non-"white")  must 


be  modeled  in  terms  of  augmented  Markovian  state  variables  to  achieve  optim¬ 
ality.  In  fact,  in  post-test  applications,  in  which  refined  ephemeris  and 
clock  data  for  each  SV  are  available,  it  may  yet  be  worthwhile  representing 
errors  in  these  quantities  as  state  variables,  although  real istically,  there 
is  little  improvement  expected  in  the  estimate  of  these  quantities.  This  is 
an  example  of  a  common  principle  —  augmentation  of  states  of  little  or  no  in¬ 
terest  to  improve  the  estimates  of  those  states  which  are  of  primary  interest. 

In  the  error  budgets  presented  in  Section  2.0,  errors  were  categorized  as  either 
slowly  or  rapidly  varying  with  time.  The  rapidly  varying  errors  are  treated 
as  white  noise  in  the  KFS.  The  slowly  varying  errors  must  optimally  be 
modeled  as  biases  or  Markovian  random  processes. 

Applications  of  KFS  to  SPS  are  widely  reported  in  the  literature.  A  few  of 
these  will  be  briefly  discussed  below. 

A  variety  of  applications  can  be  readily  envisioned  in  which  vehicle  naviga¬ 
tion  is  performed  using  a  combination  of  INS  and  GPS.  Two  examples  will 
serve  to  illustrate  most  of  the  points  of  interest;  namely,  a  real-time  cruise 
navigation  example  and  a  post-test  ballistic  missile  guidance  accuracy  analy¬ 
sis  example. 

An  INS  designed  for  long  term  real-time  application  will  benefit  considerably 
from  GPS  (Ref  [16],  pp  144-153).  These  systems  normally  are  used  to  provide 
vehicle  position  and  velocity  as  well  as  vehicle  attitude  ( i . e. ,  heading, 
pitch  and  roll)  relative  to  a  locally  level  frame  of  reference.  Without  some 
form  of  velocity  reference  these  systems  are  subject  to  Schuler  oscillations 
in  the  horizontal  plane;  and  without  an  independent  altitude  (geodetic  height) 
indicator,  they  are  unstable  in  the  vertical  channel.  Moreover,  because  of 
instrument  noise,  the  position  and  velocity  errors  slowly  diverge,  resulting 
in  the  requirement  for  frequent  position  reset. 

GPS  is  an  inherently  more  accurate  navigation  system  than  the  INS,  but  the 
INS  provides  attitude  sensing  and  is  more  responsive  to  maneuvers  than  GPS. 

When  measurements  from  the  two  systems  are  combined  in  a  real-time  Kalman 
Filter,  very  accurate  position,  velocity  and  attitude  can  result.  Moreover, 
the  INS  can  aid  GPS  tracking  during  high  dynamic  conditions,  and  the  INS  can 
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maintain  high  accuracy  for  brief  intervals  when  GPS  is  not  available. 

The  Inertial  Guidance  Systems  (IGS)  which  are  deployed  on  ballistic  missiles 
are  very  similar  in  some  respects  to  the  inertial  navigation  systems  pre¬ 
viously  discussed.  The  significant  difference  is  in  the  duration  of  the 
navigation  interval,  which  for  ballistic  missiles  is  sufficiently  short  that 
damping  of  Schuler  oscillations  and  the  vertical  channel  are  unnecessary. 

During  ballistic  missile  tests,  inertial  guidance  data  is  collected  by  tele¬ 
metry  systems  for  post-test  analysis  purposes.  If  GPS  data  is  collected  and 
telemetered  during  the  missile  flight,  from  either  a  receiver  or  translator 
onboard  the  missile,  this  data  can  be  used  as  a  basis  for  guidance  system 
analysis  in  much  the  same  fashion  as  ground  based  sensors  have  been  used  in 
the  past  (Ref  [69]). 

By  combining  telemetered  GPS  and  IG  data  In  a  KFS,  it  is  possible  to  estim¬ 
ate  initial  alignment,  accelerometer  and  gyro  instrument  error  parameters  for 
the  IGS.  Estimates  of  these  parameters  can  be  used  to  isolate  hardware  and 
software  problems  which  can  then  be  corrected  in  future  flights.  Furthermore 
the  KFS  procedure  produces  an  estimate  of  the  missile  trajectory  and  provides 
the  basis  for  an  accuracy  analysis  of  the  weapon  system. 

Perhaps  the  most  widespread  application  of  GPS  is  that  of  fixed-point  surveys. 

In  such  applications,  one  has  the  luxury  of  collecting  data  over  extended 
periods  of  time.  Data  reduction  via  KFS  (or  WIS)  can  be  employed  to  achieve 
absolute  accuracies  which  are  potentially  in  the  sub-meter  range  (Ref  [6],  p  85) 
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floppier  Preprocessing 


This  subsection  deals  with  Doppler  signal  processing  options  which  can  be  ex¬ 
ercised  prior  to  inputting  GPS  measurements  into  the  post-test  estimator. 

This  section  is  based  on  analysis  detailed  in  Ref  [18]  which  should  be  con¬ 
sulted  for  more  details. 


The  measurements  of  interest  are  the  code  and  carrier  loop  outputs  of  pseudo¬ 
range  and  delta-range.  The  pseudo-range  measurement  is  given  by 


R  -  R  ♦  R$  ♦  Rt  ♦  Rj  ♦  Rm  *  Rn 


where  R  is  the  true  range  from  the  SV  to  the  receiver,  R^  is  a  systematic 
term  which  Includes  SV  and  user  clock  errors,  Ry  is  tropospheric  delay,  Rj  is 
ionospheric  delay,  R^  is  multipath  noise,  and  R^  is  receiver  noise  and  quant¬ 
ization  e»*ror.  In  similar  fashion,  the  delta-range  measurement  is  given  by 


Cr]t  -  [RL  ♦  cy!  ♦  ty,  -  *  [rB]t  +  [rn:‘ 


where  upper  case  letters  denote  the  same  quantities  as  before,  and  lower  case 
letters  denote  quantities  which  apply  only  to  the  delta-range  measurements . 


From  a  purely  theoretical  viewpoint,  the  pseudo-range  and  delta-range  measure¬ 
ments  can  be  input  directly  to  an  optimal  estimator,  provided  stochastic  state 
variable  representations  can  be  formulated  for  each  term  (except  white  noise) 
in  the  measurement  equations.  State  variable  representation  for  R^  is 
straightforward.  It  is  probable  that  state  variable  representation  of  Ry  is 
possible  also,  but  may  not  be  worth  the  additional  estimator  complexity.  It 
is  questionable  whether  adequate  state  variable  representation  of  Rj  is 
possible.  In  any  case,  if  state  variable  representations  for  R^ ,  Ry,  and  Rj 
are  Introduced,  the  measurements  can  be  input  directly  to  the  estimator. 

In  the  case  of  delta-range,  an  additional  state  must  be  introduced  for  R(t  ), 
but  this  is  easily  done,  since  R(t  )  is  a  perfectly  correlated  function  of 
SV  and  receiver  position  states.  Thus  the  a  priori  estimate  of  R(t  )  and  the 
a  priori  covariance  of  R(t  )  are  simple  functions  of  the  a  priori  estimates 
and  covariances  of  SV  and  receiver  position  states. 
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In  practice,  Ry  and  Rj  are  not  typically  modeled  in  terms  of  state  variables 
Instead  Ry  is  generally  estimated  from  local  measurements  of  atmospheric 
conditions  (viz,  temperature,  pressure,  relative  humidity),  and  Rj  is  estim¬ 
ated  from  dual  frequency  ranging  measurements. 

The  ionospheric  delay  terms  Rn  and  at  the  respective  frequencies  f^  and 
f2  are  related  by 


where  the  approximation  improves  with  Increasing  frequencies  and  is  excel¬ 
lent  at  L-band.  The  above  relation  is  the  basis  for  dual  frequency  ionos¬ 
pheric  delay  compensation. 

When  continuous  measurements  (i.e.,  no  dropouts  or  cycle  slips)  exist  simul¬ 
taneously  on  and  l_2  for  both  pseudo-range  and  pseudo-Doppler,  a  very 
accurate  means  of  Ionospheric  delay  compensation  is  possible.  For  this  pur¬ 
pose  the  following  measurements  are  available. 

"l  ■  R  *  RS  *  "t  +  RI1  *  "m  *  *N1 

*2  •  R  +  Rs  ♦  +  Rjj  +  Rj,  +  Rn2 

0  0  0  0 

"tRn^t  +  ^t  + 

0  0 

[?2]'  -  W*  *  D4]*  *  [Rt]‘ 

0  0  0  0 

-CRI2]‘  *  CrM]ts  *  [rN2]J 

0  0  i 

It  now  follows  that 

R2  -  Rl  *  RI2  -  Rn  +  Rn2  -  RN1 


^Nl^t 

o 


and 
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F- 


Pi 


[r2]J  -  *  ^12^  +  ^rN2^t  '  ^rNl 

0  0  0  0  0  0 

Now  if  [tQ,  t^j  is  the  interval  over  which  the  above  measurements  are 
able  and  if  T  •  t^  ■  t  ,  then  using  the  notation 

tf 

/ 


<X>  2 


X(t)dt 


it  follows  that 


<R2  -  R-j  +  C»*2]t  -  Cr]]t  > 
o  o 


3  RI2^V  ■  RI1^V  +  <RN2  +  ^rN2^t  ‘  RN1  ■  ^Nl^t  > 

0  0 


Thus,  with 


K  * 


~T~T 

fVrz 


we  take  as  an  estimate  of  R^Ctg): 


RIl(to}  5  K  <  R2  ’  R1  +  ^r2^t  “  ^rl^t  > 

o  o 

3  RI1^V  +  K  <  RN2  +  ^rN2^t  ‘  RN1  "  ^rNl^t  51 

0  Q 

It  can  also  be  shown  thai  an  estimate  of  [Rj.j]£  is  given  by: 

[Rn^t  *  K  *  l-rl^t  "  ^r2^t  } 

0  0  0 

•CRn]‘  *K([rN,]J  -CrN2]J) 

0  0  0 

Combination  of  the  last  two  expressions  yields  an  estimate  for  R^(t), 


avail  - 


namely 
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Rn  ^  =  Rn (t)  +  :<  <  rn2  +  pn2  "  rni  *  rNi> 

+  K  [rN1  (t)  -  rfj2  (t)] 

Now,  Rj-j(t)  can  be  used  to  compensate  R.j  and  for  ionospheric  delay, 

and  an  estimate  of  R(tQ)  can  be  obtained  by  a  Doppler  smoothing  technique. 

Thus 

R  (V  3  <Ri  '  ^ri^t  "  2  Rn  + 

0 

■  »<'„>  *  W  *  ty'o* 

*  *  RH  *  SN1  '  CrM^t„  ■  ^rN1^t  * 

0  0 

+  K  <  RN1  -  [rN1]t  -  Rn2  +  ^rN2^t  * 

0  0 

■  »<t0)  *  as('„>  -  rtU0) 

+  <Rm  +  (1+K)  R^  -  KRn2> 

-  <[rM]t  +  O+K)  CrNl]t  -  <  CrN23t  >  . 

0  0  0 

Also  an  estimate  of  [R]£  can  be  obtained  by  compensating  [r^*  with 

-  t  0  0 

L *n]J  •  Thus 

"0 

[R]J  3  *  [Rsl*  +  CRt]J  +  0MJt  +  CrN1 3 t 

ooooo  0 

+  K  {  ^rNl^t  "  ^rN2^t  ' 

0  0 
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Finally,  combination  of  the  last  two  equations  yields 


R(t)  »  R(t)  +  R$(t)  +  RT(t)  +  <Rm  ♦  (UK)  Rn1  -  KRn2> 

-  <r„  *  (1*K)  r„,  -  KrN2> 

+  rM(t)  +  (1+K)  rN1(t)  -  KrN2(t) 

The  last  equation  provides  a  pseudo-range  measurement  which  has  been  fully 
compensated  for  ionospheric  delay.  Except  for  the  systematic  error  term  R^ 
and  the  tropospheric  error  term  Ry,  the  only  errors  are  a  bias  (due  to  init¬ 
ialization  @  t  )  and  ultra-low  phase  noise  terms.  The  bias  term 
o 

RB  =  <RM  +  RN1  *  KRN2> 

-  <rM  +  (UK)  rm  -  KrN2> 

is  approximately 

RB  *  <RM  +  ^1+K^  RN1  ' 

since  the  RMS  values  of  R^,  RN1 ,  RN2  are  greater  than  the  corresponding  values 
for  rM,  rN1 ,  rN2  by  at  least  a  factor  of  100. 

Furthermore,  the  RMS  value  of  Rg  is  simply 

aB  *  a/~T"  (ffRM  +  (UK)2  ffRMl  +  ^  jRN2)V2 


where  tr  is  the  correlation  time  of  the  code  receiver  (t^  *  1/BR  *  1  sec). 
If  T  >>  Tp,  then 


( UK)  ‘ 


K2  a2  )1/2 

K  °RN2 
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The  last  expression  for  R(t)  represents  the  most  practical  combination  of 
pseudo-range  and  delta-range  measurements  with  dual  frequency  ionospheric 
delay  compensation.  It  is  this  value  of  R(t)  which,  after  compensation  for 
tropospheric  delay,  should  be  input  to  the  optimal  estimator  to  obtain  maxi¬ 
mum  navigation  accuracy. 

The  pseudo-range  measurement  developed  above  contains  unavoidable  systematic 
error,  R$,  tropospheric  delay  compensation  error,  5Rj,  an  initialization  bias 
with  RMS  value  proportional  to  1  /J~^~  and  an  ultra-low  noise  error  which  is 
serially  uncorrelated  at  a  10  Hz  rate  (BQ  *  1/tg,  where  tQ  *  0.1  sec  is  the 
Doppler  receiver  correlation  time). 

In  contrast  to  the  above,  there  exists  a  similar  processing  technique  in  which 
Doppler  smoothing  is  employed  over  discrete  intervals  of  specified  length,  but 
only  the  initial  condition  on  each  interval  is  input  to  the  estimator  algorithm 
In  this  case,  the  expression  for  R(t^),  where  t.  denotes  the  initial  time  of 
each  discrete  interval,  is  approximated  by 

R(tt)  *  R(tf)  +  Rs(tj)  +  RT(tj) 

*  0+K)  Rfl'j  “  i  * 


where  now 


1 


■i+1 


<x>i 


ti+l  -  ti  J  t. 


X(t)dt 


and  phase  noise  terms  have  been  neglected.  With  at  =  ,t.+1  -  t . ,  it  follows 
that  the  RMS  noise  error  in  the  measurement  sequence  R(t . ) ,  i  -  0,  1 
is  given  by 


i  •  •  •  1 


/  tR  r  2  Mx i^2  „  2  .  2  T  1/2 

I  At  <■  ffQM  +  ”+*)  +  * 


RN2" 


In  this  method  the  data  rate  is  1/at,  and  the  noise  level  is  determined  by 
the  code  noise.  Thus  in  comparison  to  the  previous  method,  this  method  is 
deficient  both  with  respect  to  data  rate  and  noise  level.  Quantitatively 
this  deficiency  can  be  expressed  by  a  product  of  ratios.  Thus 


In  reality,  this  latter  method  is  equivalent  to  processing  only  the  code  meas- 
surements  in  the  presence  of  a  smooth  trajectory  constraint.  Hence,  the  latter 
method  does  not  make  very  effective  use  of  the  high  precision  and  high  data 
rate  inherent  in  the  pseudo-Ooppl er  range  measurements. 


5.0 


CONCLUSIONS  ANO  RECOMMENDATIONS 


This  document  orients  the  Range  analyst  toward  use  of  GPS  technology  for  meet¬ 
ing  Test  and  Training  Range  TSPI  requirements.  The  highly  accurate  GPS  naviga¬ 
tion  capability  has  potential  benefits  of  improved  accuracy  and  lower  costs 
for  Range  instrumentation.  The  Range  analyst  is  aided  in  evaluating  accuracy 
by  the  definition  of  error  budgets  and  accuracy  measures.  Development  of  GPS 
compatible  configurations  is  aided  by  discussions  of  system  and  data  processing 
applications  considerations.  Use  of  GPS  technology  in  Range  instrumentation 
is  a  manageable  problem  involving  the  development  of  configurations  meeting 
Test  and  Training  Ranges  requirements. 

In  order  to  expedite  availability  of  this  document  at  this  time,  it  is  not 
possible  to  develop  each  of  the  topics  to  the  In-depth  detail  it  might  deserve. 
Additional  work  should  supplement  other  projects  such  as  the  GPS  evaluation 
tests  at  Yuma,  the  Tri services  GPS  Ranges  Applications  study,  and  development 
of  the  GPS-SMILS.  It  is  assumed  that  the  Triservices  Study  will  provide  ge¬ 
neric  GPS  compatible  range  instrumentation  configurations  and  requirements  for 
receivers,  translators  and  pseudo-satellites.  The  configurations  should  in¬ 
clude  data  flow/processing  diagrams. 

This  document  addresses  technical  aspects  of  evaluating  use  of  GPS  technology 
in  range  instrumentation  but  not  cost  considerations.  The  Triservices  study 
will  develop  cost  Information  for  selected  generic  ranges  which  can  be  used 
by  the  Range  analyst  to  develop  cost  comparisons  for  specific  situations. 

Recommendations  for  follow-on  work  are: 

1)  Utilize  the  Yuma  Inverted  GPS  facility  to  determine  near/far 
effect  constraints. 

2)  Develop  lag  error  models  and  calibration  techniques  for 
dynamic  targets  carrying  various  receiver  and  translator 
configurations. 

3)  Determine  applicability  of  the  translocation  approach  to 
ionospheric  delay  correction  for  various  target  flight 

paths,  accuracy  requirements,  and  time  and  distance  separation. 


4)  Initiate  a  continuing  effort  to  monitor  GPS  error  budgets  and 
refine  accuracy  estimation  techniques  for  real-time  and  post-test 
modes  using  relative  and  absolute  techniques. 

5)  Investigate  use  of  operational  GPS  receivers  for  meeting  range 
requirements:  real-time  result  accuracy  and  timeliness, 
incorporation  of  a  recorder  or  data  link  in  the  vehicle,  and 
availability. 

6)  Identify  and  categorize  test  range  raw  data  output  requirements 
for  real-time  and  post- test  processing. 

7)  Develop  GPS  receivers,  translators  and  pseudo-satellites  for 
use  in  Range  instrumentation. 

8)  Design  and  perform  experiments  to  allow  determination  of  fast 
and  slow  error  sources  and  magnitudes  for  refinement  of  GPS 
data  processing  error  correction. 

9)  Investigate  the  use  of  C/A  code  vs  P-code  for  meeting  Range 
TSPI  requirements, 

10)  Evaluate  the  uses  of  C/A  and  P-code  on  combinations  of  L.|  and  l_2 
links  and  L1  and  l3  links  for  estimating  ionospheric  delay. 

11)  Develop  optimal  estimation  algorithms  for  GPS  applications. 


BIBLIOGRAPHY 


1.  Erb,  Hinely,  Seiders  and  White,  Assessment  of  Enhanced  TSPI  Capability 
Requirements  at  Navy  T&E  and  Fleet  Training  Ranges,  Technical  Report 
7205-81-TR-1 15,  July  t 981 ,  SRI  International,  Menlo  Park,  Calif. 

2.  J.  N.  Fredericksen,  Applicability  of  NAVSTAR  GPS  to  Test  and  Training, 

May  1980,  The  MITRE  Corp.,  McClean,  Virginia.  — 

3.  Principles  of  Error  Theory  and  Cartographic  Applications,  ACIC 
Technical  Report  No.  96,  February  1962,  United  States  Air  Force, 

St.  Louis,  Missouri. 

4.  SS-GPS-3008,  System  Specification  for  the  NAVSTAR  Global  Positioning 
System,  Code  Ident  07863,  3  March  1980. 

5.  R.  B.  Pickett,  F.  L.  Matthews  and  G.  0.  Trimble,  Use  of  GPS  and 
Telemetry  Ooppler  in  Trajectory  Measurement  Applications,  Technical 
Note  AS300-N-81-40,  August  l9fil.  Federal  Electric  Corporation-ITT, 
Vandenberg  AF3,  California. 

6.  Paul  D.  Perreault,  Civilian  Receivers  Navigate  by  Satellite,  January  1981, 
Microwave  Systems  News,  Pgs  61-93,  Stanford  Telecommunications,  Inc. 

7.  Wayne  F.  Brady  and  Paul  S.  Jorgenson,  Worldwide  Coverage  of  the  Phase  II 
NAVSTAR  Satellite  Constellation.  The  Aerospace  Corp.,  El  Segundo, 
California. 

8.  A.  J.  MacMillan,  NAVSTAR  Global  Positioning  System  Flight  Test  Program 
Overview,  The  Aerospace  £orp.,  El  Segundo,  California. 

9.  S.  P.  Teasley,  W.  M.  Hoover,  and  C.  R.  Johnson,  Differential  GPS 
Navigation,  Texas  Instruments,  Inc. 

10.  Phil  Ward,  Advanced  NAVSTAR  GPS  Multiplex  Receiver,  Texas  Instruments, 

Inc.,  Lewisville,  Texas.  - 

11.  Kenneth  Putkovitch,  USNO  GPS  and  TRANSIT  Programs,  U.  S.  Naval  Observa¬ 
tory,  Washington,  0 .  ?TI 

12.  S.  A.  Book,  W.  F.  Brady,  and  Mazaika,  The  Nonuniform  GPS  Constellation, 
Mazaika,  The  Aerospace  Corp.,  El  Segundo,  California. 

13.  Paul  S.  Jorgenson,  Combined  Pseudo  Range  and  Ooppler  Positioning  for  the 
Stationary  NAVSTAR  User,  The  Aerospace  Corp.  ,  El  Segundo,  CA 

14.  W.  C.  Melton,  Global  Positioning  System  Measures  Time  of  Arrival, 

Stanford  Telecommunications,  Inc.,  Microwave  Systems  News, 

January  1982,  Vol .  12,  No.  1,  pgs  136-142. 

15.  Phase  I  NAVSTAR/GPS  Major  Field  Test  Objective  Report,  Ephemeris  and 
Space  Vehicle  Clock  Prediction  Accuracy,  4  May  1979,  NAVStAR/6Ps  Joint 
Program  Office,  SAMS0,  Los  Angeles,  California. 


16.  Vol.  25,  No.  2,  Global  Positioning  System,  Reprint  of  Summer  1978 
Institute  of  Navigation  Journal,  1980.  ™ 

17.  Caroline  F.  Lerory,  The  Impact  of  GRS  80  on  DMA  Products,  Defense 
Mapping  Agency,  Washington,  D.  C. 

18.  Or.  R.  A.  Brooks,  Ionospheric  Refraction  Compensation  in  GPS  Applications, 
Technical  Note  AS300-N-32-07,  January  1982,  Systems  Performance  Analysis 
Department,  Federal  Electric  Corp.-ITT,  Vandenberg  AFB,  California. 

19.  P.  8.  Levine  and  H.  L.  Jones,  GPS  Range  Instrumentation  Equipment  Issues 

and  Applications,  30  January  1982,  The  Analytic  Sciences  Corp. ,  ” 

Reading,  Massachusetts . 

20.  P.  S.  Jorgenson,  Ionospheric  Measurements  from  NAVSTAR  Satellites, 

December  1978,  The  Aerospace  Corp.,  El  Segundo,  California. 

21.  A.  J.  Mallinckrodt,  Group  and  Phase  Phenomena  in  an  Inhomogeneous 
Ionosphere,  Santa  Ana,  California. 

22.  Phase  I  NAVSTAR/GPS  Major  Field  Test  Objective  Report,  Tropospheric 
Correction,  4  May  1979.  Navstar/GPS  Joint  Program  bffice, 

Los  Angeles,  California. 

23.  Phase  I  NAVSTAR/GPS  Major  Field  Test  Objective  Report,  Ionospheric 
Retraction,  4  May~T577,  NAvstAR/GPS  Joint  Program  Office,  Los  Angeles, 
California. 

24.  V.  L.  Pisacane  and  M.  M.  Feen,  Ionospheric  Effects  on  Trans ionospheric 
Measurements  of  Range  and  Range  Difference  TG  1267,  October  1974, 

Johns  Hopkins  University,  Applied  Physics  Laboratory,  Silver  Spring, 
Maryland. 

25.  J.  A.  Klobuchar,  M.  A.  and  0.  A.  Pearson,  A  Preliminary  Evaluation  of 
the  Two-Frequency  Ionospheric  Correction  for  the  NAvSTAR-GlobaT 
Positioning  System,  Air Force  Geophysics  Laboratory,  Hanscom  A?B,  The 
Aerospace  Corp.,  Los  Angeles,  Calif. 

26.  J.  A.  Klobuchar,  Ionospheric  Time  Delay  Corrections  for  Advanced 
Satellite  Ranging  Systems,  Air  Force  Geophysics  Laboratory,  Hanscom  AF3,  MA. 

27.  J.  A.  Klobuchar  and  J.  M.  Johnson,  Correlation  Distance  of  Mean 
Daytime  El ectron  Content,  AFGL-TR-77-0185,  Air  force  Surveys  in 
Geophysics,  No.  3^3,  22  August  1977,  Air  Force  Geophysics  Laboratory, 

Hanscom  AFB,  MA. 

28.  John  A.  Klobuchar,  A  First-Order,  Worldwide,  Ionospheric,  Time-Qelav 
Algorithm,  AFCRL-TR-75-0502,  Air  Force  Surveys  in  Geoohysics,  Mo.  324, 
25~September  1975,  Air  Force  Geoohysics  Laboratory,  Hanscom  AFB,  MA. 

29.  E.  G.  31akcwe11 ,  J.  S.  Cline,  E.  A.  Erb,  and  J.  R.  Olmstead,  GPS-SMILS 
Concept  Definition,  WSMC  TR  30-1,  March  1980,  SRI  International,  Menlo 
Park,  California. 


30.  John  H.  Painter,  Designing  Pseudorandom  Coded  Ranging  Systems,  Motorola, 
Inc.,  Scottsdale,  Arizona. 


31.  S.  C.  Gupta  and  0.  H.  Painter,  Correlation  Analysis  of  Linear. 

Processed  Pseudo  Random  Sequences,  December  1^65,  Motorola.  Inc. , 
Scottsdale,  AZ.  ~ 

32.  Rodney  J.  Boehm,  Simulation  of  a  Kalman  Filter  for  GPS.  June  1979, 

Texas  A&M  University,  College  Station,  Texas. 

33.  Steven  B.  Hyde,  A  NonLinear  Technique  to  Compensate  for  Range  Bias 
Error  in  a  low  Cost  GPS  Set,  TCSL  Memo  ?  &009,  August  1980.  Texas 
A&M  University,  College  Station,  Texas. 

34.  Rodney  J.  Boehm,  An  Adaptive  Gain  Technique  for  the  GPS  Extended 
Kalman  Filter.  TCSL  Memo  7910,  November  1979,  Texas  A&M  University, 

College  Station,  Texas. 

35.  M.  D.  Eggers  and  0.  M.  Ballock,  Computer  Simulation  of  Ionospheric 
Wave,  Propagation  for  Detection  of  Range  Error  for  Satellite  Navigation 
Systems,  TCsL  Research  Memorandum  Bl-Oo,  May  19,  1981,  Texas  A&M 
University,  College  Station,  Texas 

36.  Col.  Donald  W.  Hederson,  Joseph  A.  Strada,  Lt.  Com.,  NAVSTAR  Field 
Test  Results.  MX-TM-3316-80A,  July  20,  1980,  USAF  Space  Division,  El 
Segundo,  CA. 

37.  L.  J.  Jacobson,  and  V.  Calbi,  Engineering  Development  of  NAVSTAR 

GPS  User  Equipment,  April  1981,  Magnavox  Advanced  Products  and  Systems 
Company,  'Torrance,  California 

38.  Richard  Enossen,  Low-Cost  GPS  Navigation  Receiver  for  General  Aviation, 
February  1980,  Magnavox  Government  and  Industrial  Electronics  Company, 
Torrance,  California 

39.  Technical  Data  Sheet,  Collins  GPS  User  Equipment  Manoack/Vehicular 
Application,  Rockwell  International,  Collins  Government  Avionics  Division. 


40.  Technical  Data  Sheet,  Collins  GPS  User  Equipment  Surface  Ship  Application, 
Rockwell  International,  Collins  Government  Avionics  Division. 

41.  Technical  Data  Sheet,  Collins  GPS  User  Equipment  Submarine  Application, 
Rockwell  International,  Collins  Government  Avionics  Division.  ~~ 

42.  Technical  Data  Sheet,  Collins  GPS  User  Equipment  Rotary  Wing  Aircraft 
Aopl ication,  Rockwell  International,  Collins  Government  Avionics  Division. 

43.  Technical  Data  Sheet,  Collins  GPS  User  Equipment  Transport/Patrol  Aircraft 
Application,  Rockwell  International,  Collins  Government  Avionics  Division. 


*6.  G.  Trimble,  R.  Weigel  and  S.  Cresswell,  Analysis  of  the  Ballistic 
Receiver  Evaluation  Wafer  flight  Test.  OP  7434  (PVM-19),  June  1981 
Technical  Report  Number  AS3Q0-T-31-27,  Systems  Performance  Analysis 
Department,  Federal  Electric  Corporation,  ITT,  Vandenberg  AF3,  CA. 

47.  G.  Trimble,  R.  Weigel  and  S.  Cresswell ,  Analysis  of  the  Ballistic 
Receiver  Evaluation  Wafer  Flight  Test,  0P1767  (PVM-18)  technical  Report 
AS300-T-80-33,  November  1980,  Systems  Performance  Analysis  Department, 
Federal  Electric  Corporation,  ITT,  Vandenberg  AFB,  California 

48.  Multi-Ob.lect  Tracking,  Range  Instrumentation  Symposium,  3-4  March  1981, 
Volume  1  -  Unclassified 

49.  GPS-RTP,  ICD-4,  MBRS  Functional  Characteristics,  Interface  Control, 

June  3  1978  plus  Revisions  SAMSO  AFSC,  Norton  AFB,  CA. 

50.  Space  Vehicle  Nav  Subsystem  and  NTS  PRN  Navigation  Assembly/User 
System  Segment  and  Monitor  Station.  August  T97j.  Rockwell  international 
Co rp..  Space  Oivision,  Downey,  CA. 

51.  R.  L.  Barkley,  Jr.,  and  R.  T.  Herzog,  Global  Positioning  System  Mobile 
Missile  Tracking  Platform  Feasibility  Study,  301 96-6002-TU-00,  30  June 
1977,  TRW  Defense  and  Space  Systems  Group,  Redondo  Beach,  CA. 

52.  Nicholas  Stilwell,  Edward  H.  Martin,  and  Jack  Moses,  GPS  Study  for  Range 
Instrumentation,  Final  Report  for  Space  and  Missile  Test  Center  (XREA), 

30  March  1 §77 , Magnavox  Government  and  Industrial  Electronics  Co., 
Advanced  Products  Division,  Torrance,  CA. 

53.  R.  1.  Bogusch,  F.  W.  Guigliano,  D.  L.  Knepp,  and  A.  H.  Michelet, 

Frequency  Selective  Propagation  Effects  on  Spread-Spectrum  Receiver 
Tracking,  Mission  (Research  (Dorp.,  Santa  Barbara.  CA:  Proceedings  of 
the  IEEE,  Vol  69,  No.  7,  July  1981,  pgs  787-796. 

54.  W.  R.  Fried,  A  comparative  Performance  Analysis  of  Modern  Ground-Based, 
Air-Based,  and  Satellite-Based  Radio  Navigation  Systems,  Hughes  Aircraft 
Co.,  Fullterton,  CA;  Journal  of  the  Institute  of  Navigation,  Vol  24,  No.l, 
Spring  1977,  pgs  48-58. 

55.  Technical  Digest,  Johns  Hopkins  API,  Navy  Navigation  Satellite  System 
(Transit)  Issue;  January-March  1981,  Vol  2,  No.  1. 

56.  Jack  K.  Holmes,  Coherent  Spread  Spectrum  Systems.  John  Wiley  &  Sons, 

New  York,  N.  Y.  1982. 

57.  James  J.  Spilkor,  Digital  Communications  by  Satellite.  Prentice  Hall, 
Englewood  Cliffs,  NJ  19??. 

58.  John  0.  Cardall  and  Richard  S.  Crossen,  Civil  Application  of  Differ¬ 
ential  GPS,  Oct.  1981,  Magnavox  Advanced  Products  and  Systems  Company, 
Torrance,  CA. 


59.  Jacques  Beser  and  Bradford  W.  Parkinson,  The  Application  of  NAVSTAR 
Differential  GPS  in  the  Civilian  Community,  Intermetrics ,  Inc., 
Huntington  Beach,  CA. 


in  Crystal  Oscillators 


,  Jet  Propulsion  Laboratory  Report 


73.  P.  S.  Jorgenson,  "Normalized  Accuracy  of  the  NAVSTAR/Global  Positionin 
The  Aer 
ry  1978 


A.  LINEAR  ESTIMATION 

Consider  a  linear  stochastic  system  described  by 

0)  xjt,  .  «,x,  *  U(,  i  •  <3,  1 . . 

(2)  yi  3  H.x.  +  v.,  i  *  0,  1 . N  . 

In  these  equations  x  is  the  state  vector,  and  y  is  the  measurement  vector, 
tu. ,  i  3  0,  1,  ....  N-1}  is  a  sequence  of  random  vectors  called  the  state 
noise  process,  and  {v,.,  i  »  0,  1 ,  . . . ,  N}  is  a  similar  sequence  called  the 
measurement  noise  process.  These  processes  are  assumed  to  be  zero  mean, 
sequentially  uncorrelated,  and  mutually  uncorrelated  with  each  other  and  xQ. 
Mathematically  these  assumptions  are  expressed  as  follows: 

Mij’  1,13  “  °’  1 . N"1  * 

Ri6ij;  j  3  °’  1 . N 


3  0,  E(x0v])  3  0  ; 


Q  and  R  are  the  state  noise  and  measurement  noise  covariance  matrices,  re¬ 
spectively.  To  complete  the  system  description,  it  is  assumed  that  the 
a  priori  mean  and  covariance  of  xQ,  denoted  by  x’  and  P‘,  respectively,  are 
also  specified. 


E(u.)  3  0,  E^uT)  3 
E(vi)  3  0,  E(v1 vj)  * 
and 

E(u.vT>  *  0,  E(xqu^} 

where 


The  estimation  procedures  to  be  considered  in  the  remainder  of  this  appendix 
will  apply  to  the  linear  system  given  by  (1)  and  (2). 


A.  1 


Optimal  Linear  Estimation 


For  the  system  (1),  (2),  let  x(i|j)  denote  a  function  which  is  in  the  form 

of  a  constant  plus  a  linear  function  of  the  measurement  set  y  ,  ....  y.  and 

**  *  J 

which  has  the  property  E[x(i|j)]  *  Z(x.).  Such  a  function  is  said  to  be  a 

linear  unbiased  estimate  of  x,  given  y  ,  . . . ,  y . . 

'  ^  J 

Let  x(i|j)  be  a  linear  unbiased  estimate  of  x..  given  yQ,  ....  y  j ,  and  suppose 
x  ( i  |  j )  has  the  property* 

ECllidlJ)  -  X,||2]  <  EClIxdlJ)  -  xf||2] 

for  all  linear  unbiased  estimates  x(i|j)  of  xi  given  yQ . y^..  Then 

x(i|j)  is  said  to  be  an  optimal  linear  estimate  of  x^  given  yQ,  ....  y^.. 

It  can  be  shown  that  an  optimal  linear  estimate  of  xi  given  yQ,  ....  y. 
always  exists  and  is  unique.  The  notation  x(i|j)  will  be  used  exclusively 
to  denote  the  optimal  estimate  defined  above,  and  the  notation  P(i|j)  will 

A 

be  used  to  denote  the  error  covariance  of  x(i|j)  defined  by 
P ( i I J )  S  E[(x(i|j)  -  x^)  (x(i | j )  -  x1 )T3  . 

P(i|j)  is  also  called  the  state  covariance  of  x^  given  yQ,  y^. 

A. 2  Kalman  Estimation 

A  recursive  procedure  for  realization  of  the  optimal  linear  estimator  for 
the  system  (1),  (2)  has  been  developed  by  Kalman  [3],  [4].  The  procedure 
consists  of  two  stages.  The  first  stage  employs  a  filter  algorithm,  while 
the  second  uses  a  smoother  algorithm**. 


*  The  notation  | |zj 


1  T 

=  z  z  denotes  the  ordinary  Euclidean  norm  of  z. 


**The  terminology  employed  here  is  due  to  N.  Wiener  and  has  been  adopted  by  R. 
Kalman.  An  estimator  which  estimates  x,  given  measurements  with  indices  up 
to  and  including  j  is  called  a  fi 1  ter  ‘if  i  ■  j;  it  is  called  a  predictor 
if  i  >  j,  and  it  is  called  a  smoother  if  i  <  j. 


m 


I*. 

Apolication  of  the  Kalman  filter  to  the  system  (1),  (2)  yields  x(i|i) 
and  P(iji)  =  P.  by  recursion.  The  algorithm  is 


(3) 

Ki  “ 

a;h>,.p:h(  *  Rjj 

(A) 

X1  * 

-  v; 

(5) 

Pi  * 

p;  -  (  ■ 

(6) 

>> 

xi+1 

*  «i*i  • 

(7) 

p;+i 

*  *  V  1 

-1 


The  smoother  algorithm  is  initiated  when  the  filter  stage  is  complete.  The 
smoother  uses  the  filter  outputs  in  a  recursive  process,  which  runs  in  reverse 

a 

order  to  the  filter  recursion,  to  compute  x(i|N)  and  P(i|N),  i  *  0,  1,  ....  N. 
The  algorithm  for  smoothing  is  given  by 


M  w  *  view'  • 

(9)  x(i|N)  *  x.  +  Aj[x(i+1 |N)  -  xT+1]  , 

■  (10)  P(1|N)  *  Pf  ♦  A^PO+llN)  -  Pj+1]AT,  i  *  N-l . 0. 
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GLOSSARY 


ACRONYMS 

SECTION 

BET 

= 

Best  Estimate  of  Trajectory 

4.81 

bps 

= 

Bits  Per  Second 

4.8.2 

C/A 

= 

Coarse  or  Clear  Acquisition  Code,  2  MHz 

2.0 

CEP 

= 

Circular  Error  Probable 

3.0 

CPU 

s 

Computer  Processor  Units 

4. 9. 1.3 

CRPA 

= 

Controlled  Reception  Pattern  Antenna 

4. 9. 1.1 

CRT 

= 

Cathode  Ray  Tube 

4.8.1 

dB 

= 

Decibels 

2.0 

OBI 

s 

Decibels 

2.0 

DBW 

= 

Decibels  Referred  to  1  w 

2.0 

DCO 

3 

Digitally  Controlled  Oscillator 

4.8.2 

DOP 

» 

Dilution  of  Precision 

3.0 

ETR 

3 

Eastern  Test  Range 

4.7.1 

FMI 

= 

Flexible  Modular  Interface 

4. 9. 1.3 

FOM 

= 

Figure  of  Merit 

3.0 

FRPA 

3 

Fixed  Reception  Pattern  Antenna 

4. 9. 1.1 

9 

= 

Unit  of  Gravity 

4.9.3 

GPS 

= 

Global  Positioning  System 

2.0 

GPS-SMILS 

3 

GPS-Sonobuoy  Missile  Impact  Location  System 

4.7.1 

HDUE 

3 

High  Dynamic  User  Equipment 

4.7.2 

Hz 

3 

Hertz,  Cycles  per  Second 

2.0 

I  CBM 

= 

Intercontinental  Ballistic  Missile 

4.0 

IF 

3 

Intermediate  Frequency 

4. 7. 3. 2. 3 

IG 

3 

Inertial  Guidance 

4.10 

IGS 

3 

Inertial  Guidance  System 

4.9.3 

INS 

= 

Inertial  Navigation  System 

4.9.3 

JHU-APL 

3 

Johns  Hopkins  University-Applied  Physics  Lab 

4.7.1 

JPO 

3 

Joint  Program  Office 

1.1 

KFS 

3 

Kalman  Filter/Smoother 

4.10 

kHz 

3 

KiloHertz 

4. 7. 3. 2. 2 
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u'l . 

ACRONYMS 

SECTION 

h’L2 

3 

Designations  for  GPS  SV  Signals 

2.0 

(J 

LEP 

3 

Linear  Error  Probable 

3.0 

•*  ’• 

LRU 

3 

Line  Replaceable  Units 

4.9.1 

MAE 

3 

Missile  Accuracy  Evaluation 

4.8 

MBRS 

S 

Missile  Borne  Receiver  Set 

4.7.1 

MHZ 

3 

Mega  Hertz 

2.0 

MX 

3 

Missile  X 

4.8 

P 

3 

Precision  Code,  20  MHz 

2.0 

PCM-NRZ 

3 

Pulse  Code  Mod  -  Non  Return  To  Zero 

4. 7. 3. 2. 2 

PC(R) 

3 

Circular  Error  Probability 

3.0 

Pe(K;3) 

3 

Elliptical  Error  Probability  in  3  Dimensions 

3.0 

Pl(R) 

3 

Linear  Error  Probability 

3.0 

PM 

3 

Phase  Modulation 

4. 7. 3. 2. 3 

1 

PN 

3 

Pseudo-Noise 

4.8.2 

PSK 

3 

Phase  Shift  Keyed 

4. 7. 3. 2. 2 

PS(R) 

3 

Spherical  Error  Probability 

3.0 

qpsk 

3 

Quad  Phase  Shift  Keyed 

4. 7. 3. 2. 2 

a 

R&O 

3 

Research  and  Development 

4.8.1 

RF 

3 

Radio  Frequency 

2.0 

RH 

3 

Right  Hand  (Polarization) 

2.0 

RMS 

3 

Root  Mean  Square 

4.2.1 

a 

SEP 

3 

Spherical  Error  Probable 

3.0 

jC 

3 

Circular  Standard  Error 

3.0 

* 

aS 

3 

Spherical  Standard  Error 

3.0 

SNR,  S/N 

3 

Signal  to  Noise  Ratio,  dB 

4. 7. 3. 2. 2 

-  •* 

SPS 

3 

Samples  Per  Second 

4.8 

m 

SRU 

3 

Shop  Replaceable  Units 

4.9.1 

;. 

SV 

3 

GPS  Satellite  "Space  Vehicle" 

2.0 

TTFF 

3 

Time  to  First  Fix 

4.9.4 

UERE 

3 

User  Equivalent  Range  Error 

2.0 

-i 

r— . 

VCO 

3 

Voltage  Controlled  Oscillator 

4.8.2 

;  -• 

MLS 

3 

Weighted  Least  Squares 

4.10 

WSMC 

3 

Western  Space  and  Missile  Center 

PREFACE 

r-. 

MTR 

3 

Western  Test  Range 

4.7.1 

b  ■■ 
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